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Abstract 



Quantum wire networks are novel artificial nano-objects that represent a two dimen- 
sional (2D) grid formed by superimposed crossing arrays of parallel conducting quantum 
wires, molecular chains or metallic single-wall carbon nanotubes. Similar structures arise 
naturally as crossed striped phases of doped transition metaloxides. The mechanical 
flexibility of the networks, the possibility of excitation of some of their constituents (a 
nanotube or a single wire) by external electric field, and the existence of bistable confor- 
mations in some of them (molecular chain) make the networks one of the most attractive 
architectures for designing molecular-electronic circuits for computational application. 

Since such networks have the geometry of crossbars, we call them "quantum crossbars" 
(QCB). Spectral properties of QCB cannot be treated in terms of purely ID or 2D electron 
liquid theory. A constituent element of QCB (quantum wire or nanotube) possesses the 
Luttinger liquid (LL) fixed point. A single array of parallel quantum wires is still a 
LL-like system qualified as a sliding LL phase provided only the density- density or/and 
current-current interaction between adjacent wires is taken into account. Two crossing 
arrays (QCB) coupled only by capacitive interaction in the crosses have similar low-energy, 
long-wave properties characterized as a crossed sliding LL phase. QCB with electrostatic 
interaction in the crosses possess a cross-sliding Luttinger liquid (CSLL) zero energy fixed 
point. 

In this Thesis we develop a theory of interacting Bose excitations (plasmons) in a 
superlattice formed by m crossed arrays of quantum wires. The subject of the theory is 
the spectrum of excitations and response functions in the 2D Brillouin zone so it goes far 
beyond the problem of stability of the CSLL fixed point. 

In the first part we analyze spectrum of boson fields and two-point correlators in double 
(m = 2) square QCB [double tilted and triple (m = 3) QCB are considered in Appendices 
O and El respectively] . We show that the standard bosonization procedure is valid, and 
the system behaves as a cross-sliding Luttinger liquid in the infrared limit, but the high 
frequency spectral and correlation characteristics have either ID or 2D nature depending 
on the direction of the wave vector in the 2D elementary cell of reciprocal lattice. As a 
result, the crossover from ID to 2D regime may be experimentally observed. It manifests 
itself as appearance of additional peaks of optical absorption, non-zero transverse space 
correlators and periodic energy transfer between arrays ("Rabi oscillations"). 

In the second part, the effectiveness of infrared (IR) spectroscopy is studied. IR 
spectroscopy can be used as an important and effective tool for probing QCB at finite 



frequencies far from the LL fixed point. Plasmon excitations in QCB may be involved in 
resonance diffraction of incident electromagnetic waves and in optical absorption in the 
IR part of the spectrum. The plasmon velocity is much smaller than the light velocity. 
Therefore, an infrared radiation incident on an isolated array, cannot excite plasmons at 
all. However in QCB geometry, each array serves as a diffraction lattice for its partner, 
giving rise to Umklapp processes of reciprocal super-lattice vectors. As a result, excitation 
of plasmons in the center of the Brillouin zone (BZ) occurs. 

To excite QCB plasmons with non-zero wave vectors, an additional diffraction lattice 
(DL) coplanar with the QCB can be used. Here the diffraction field contains space 
harmonics with wave vectors perpendicular to the DL that enable one to eliminate the 
wave vector mismatch and to scan plasmon spectrum within the BZ. In the general case, 
one can observe single absorption lines forming two equidistant series. However, in case 
where the wave vector of the diffraction field is oriented along some resonance directions, 
additional absorption lines appear. As a result, an equidistant series of split doublets can 
be observed in the main resonance direction (BZ diagonal). This is the central concept 
of dimensional crossover mentioned above with direction serving as a control parameter. 
In higher resonance directions, absorption lines form an alternating series of singlets and 
split doublets demonstrating new type of dimensional crossover. The latter occurs in a 
given direction with frequency as a control parameter. 

In the third part, dielectric properties of QCB interacting with semiconductor sub- 
strate are studied. It is shown that a capacitive contact between QCB and a semicon- 
ductor substrate does not destroy the Luttinger liquid character of the long wave QCB 
excitations. However, the dielectric losses of a substrate surface are drastically modified 
due to diffraction processes on the QCB superlattice. QCB-substrate interaction results 
in additional Landau damping regions of the substrate plasmons. Their existence, form 
and the spectral density of dielectric losses are sensitive to the QCB lattice constant and 
the direction of the wave vector of the substrate plasmon. Thus, the dielectric losses in 
the QCB-substrate system serve as a feasible tool for studying QCB spectral properties. 

In the fourth part we formulate the principles of ultraviolet (UV) spectroscopy, i.e., 
Raman-like scattering. UV scattering on QCB is an effective tool for probing QCB spec- 
tral properties, leading to excitation of QCB plasmon(s). Experimentally, such a process 
corresponds to sharp peaks in the frequency dependence of the differential scattering cross 
section. The peak frequency strongly depends on the direction of the scattered light. As a 
result, ID — > 2D crossover can be observed in the scattering spectrum. It manifests itself 
as a splitting of single lines into multiplets (mostly doublets). The splitting magnitude 
increases with interaction in the QCB crosses, while the peak amplitudes decrease with 



ii 



electron-electron interaction within a QCB constituent. 

The following novel results were obtained in the course of this research: 

• It is shown that the bosonization procedure may be applied to the Hamiltonian of 2D 
quantum networks. QCB plasmons may have either ID or 2D character depending 
on the direction of the wave vector. The crossover from ID to 2D regime may be 
experimentally observed. Indeed, due to inter-wire interaction, unperturbed states, 
propagating along the two arrays are always mixed, and transverse components of 
correlation functions do not vanish. For quasi-momenta near the resonant line of the 
BZ, such mixing is strong, and the transverse correlators possess specific dynamical 
properties. One of the main effects is the possibility of a periodic energy transfer 
between the two arrays of wires. 

• The principles of spectroscopic studies of the excitation spectrum of quantum cross- 
bars are established, which possesses unique property of dimensional crossover. The 
plasmon excitations in QCB may be involved in resonance diffraction of incident 
light and in optical absorption in the IR part of the spectrum. One can observe 
ID 2D crossover behavior of QCB by scanning an incident angle. The crossover 
manifests itself in the appearance of a set of absorption doublets instead of the set 
of single lines. At special directions, one can observe new type of crossover where 
doublets replace the single lines with changing frequency at a fixed direction of a 
wave vector. 

• It is shown that a capacitive contact between QCB and semiconductor substrate 
does not destroy the LL character of the long wave excitations. However, quite 
unexpectedly the interaction between the surface plasmons and plasmon-like exci- 
tations of QCB essentially influences the dielectric properties of a substrate. First, 
combined resonances manifest themselves in a complicated absorption spectra. Sec- 
ond, the QCB may be treated as the diffraction grid for a substrate surface, and an 
Umklapp diffraction processes radically change the plasmon dielectric losses. So the 
surface plasmons are more fragile against interaction with superlattice of quantum 
wires than the LL plasmons against interaction with 2D electron gas in a substrate. 

• The principles of inelastic UV Raman spectroscopy of QCB are formulated. An 
effective Hamiltonian for QCB-light interaction is expressed via the same boson 
fields as the Hamiltonian of the QCB themselves. One can observe ID — > 2D 
crossover of QCB by scanning an scattered angle. The crossover manifests itself in 
the appearance of multiplets (mostly doublets) instead of single lines. 

iii 
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Chapter 1 

Background and Objectives 



The behavior of electrons in arrays of one-dimensional (ID) quantum wires was recog- 
nized as a challenging problem soon after the consistent theory of elementary excitations 
and correlations in a Luttinger liquid (LL) of interacting electrons in one dimension was 
formulated (see pQ for a review). In contrast to the Fermi liquid (FL) theory j2j EJ, one 
dimensional electron liquids exhibits the following properties: 

• There is no elementary fermionic quasi-particles, the generic excitations are bosonic 
fluctuations. 

• Charge and spin quasi-particles are spatially separated and move with different 
velocities (charge-spin separation). 

• The correlations between these excitations show up as an interaction-dependent 
non-universal power law. 

There are several possible regimes in which a ID electron liquid can exist j3]. First, there 
is an insulating regime, where charge and spin excitations are gapped. Second, there is a 
conducting (Tomonaga-Luttinger) regime where the charge sector is gapless. In this case 
the spin sector is either gapped (Luther-Emery regime) or gapless. 

One of the fascinating challenges is a search for LL features in higher dimensions [3] . 
Although the Fermi liquid state seems to be rather robust for D > 1, a possible way 
to retain some ID excitation modes in 2D and even 3D systems is to consider highly 
anisotropic objects, in which the electron motion is spatially confined in the major part 
of real space (e.g., it is confined to separate linear regions by potential relief). One may 
hope that in this case, weak enough interaction does not violate the generic long-wave 
properties of the LL state. 

Recent achievements in material science and technology have led to fabrication of an 
unprecedented variety of artificial structures that possess properties never encountered in 
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"natural" quantum objects. One of the most exciting developments in this field is fabri- 
cation of 2D networks by means of self-assembling, etching, lithography and imprinting 
techniques [HIE!- Another development is the construction of 2D molecular electronic cir- 
cuits [8] where the network is formed by chemically assembled molecular chains. Arrays 
of interacting quantum wires may be formed in organic materials and in striped phases of 
doped transition metal oxides. Artificially fabricated structures with controllable config- 
urations of arrays and variable interactions are available now (see, e.g., Refs. [H] ITU] ITT]). 
Such networks have the geometry of crossbars, and bistable conformations of molecular 
chains may be used as logical elements [T2] . Especially remarkable is a recent experimen- 
tal proposal to fabricate 2D periodic grids from single-wall carbon nanotubes (SWCNT) 
suspended above a dielectric substrate [H]. The possibility of excitation of a SWCNT by 
external electric field together with its mechanical flexibility makes such a grid formed by 
nanotubes an excellent candidate for an element of random access memory for molecular 
computing. 

From a theoretical point of view, such double 2D grid, i.e., two superimposed crossing 
arrays of parallel conducting quantum wires [TH] HH [THj or nanotubes [IB] , represents a 
unique nano-object - quantum crossbars (QCB). Its spectral properties cannot be treated 
in terms of purely ID or 2D electron liquid theory. A constituent element of QCB 
(quantum wire or nanotube) possesses the Luttinger liquid (LL)-like spectrum ^3 EE] • 
The inter-wire interaction may transform the LL state existing in isolated quantum wires 
into various phases of 2D quantum liquid. The most drastic transformation is caused by 
inter-wire tunneling in arrays of quantum wires with intra-wire Coulomb repulsion. The 
tunneling constant rescales towards higher values for strong intra-wire interaction, and 
the electrons in an array transform into 2D Fermi liquid fHl ED]- The reason for this 
instability is the orthogonality catastrophe, i.e., the infrared divergence in the low-energy 
excitation spectrum that accompanies the inter-wire hopping processes. 

Unlike inter-wire tunneling, density-density or current-current inter-wire interaction 
does not modify the low-energy behavior of quantum arrays under certain conditions. In 
particular, it was shown recently fU\ [221 EH1 EE] that "vertical" interaction which depends 
only on the distance between the wires, imparts the properties of a sliding phase to 2D 
array of ID quantum wires. Such LL structure can be interpreted as a quantum analog 
of classical sliding phases of coupled XY chains [21] . Recently, it was found [23 1201 that 
a hierarchy of quantum Hall states emerges in sliding phases when a quantizing magnetic 
field is applied to an array. 

Similar low-energy, long-wave properties are characteristic of QCB as well. Its phase 
diagram inherits some properties of sliding phases in case when the wires and arrays are 



4 



coupled only by capacitive interaction (T^j. When the inter-array electron tunneling is 
possible, say, in crosses, dimensional crossover from LL to 2D FL occurs [23 EE]- If 
tunneling is suppressed and the two arrays are coupled only by electrostatic interaction 
in the crosses, the system possesses the LL zero energy fixed point [23]. 

The physics of dimensional crossover is quite well studied, e.g., in thin semiconductor 
or superconductor films where the film thickness serves as a control parameter that governs 
the crossover (see e.g. Ref. |2H1ES]). It occurs in strongly anisotropic systems like quasi- 
one-dimensional organic conductors (SHI or layered metals [SU E21 EH1 OH EHj In the 
latter cases, temperature serves as a control parameter and crossover manifests itself in 
inter-layer transport. In metals, the layers appear "isolated" at high temperature, but 
become connected at low temperatures to manifest 3D conducting properties. 

The most promising type of artificial structures where dimensional crossover is ex- 
pected is a periodic 2D system of m crossing arrays of parallel quantum wires or carbon 
nanotubes. We call it "quantum crossbars" (QCB). Square grids of this type consisting 
of two arrays were considered in various physical contexts in Refs. [T3] ITi] EH [2Z] • In 
Refs. |lf)[ 127] the fragility of the LL state against inter-wire tunnelling in the crossing 
areas of QCB was studied. It was found that a new periodicity imposed by the inter-wire 
hopping term results in the appearance of a low-energy cutoff A; ~ Hv/a where v is the 
Fermi velocity and a is the period of the quantum grid. Below this energy, the system is 
"frozen" in its lowest one-electron state. As a result, the LL state remains robust against 
orthogonality catastrophe, and the Fermi surface conserves its ID character in the cor- 
responding parts of the 2D Brillouin zone (BZ). This cutoff energy tends to zero at the 
points where the one-electron energies for two perpendicular arrays and €k 2 become 
degenerate. As a result, a dimensional crossover from ID to 2D Fermi surface (or from 
LL to FL behavior) arises around the points ep 1 = ep 2 . 

Unlike inter-wire tunneling, the density- density or current-current inter-wire interac- 
tion does not modify the low-energy behavior of quantum arrays under certain conditions. 
In particular, it was shown recently fH\ [221 EHj that "vertical" interaction which depends 
only on the distance between the wires, imparts the properties of a sliding Luttinger 
liquid phase to 2D array of ID quantum wires. Such LL structure can be interpreted 
as a quantum analog of classical sliding phases of coupled XY chains [23]. Recently, it 
was found j2E| that a hierarchy of quantum Hall states emerges in sliding phases when a 
quantizing magnetic field is applied to an array. Similar low-energy, long-wave properties 
are characteristic of QCB as well. Its phase diagram inherits some properties of sliding 
phases in case when the wires and arrays are coupled only by capacitive interaction [16 . 
If tunneling is suppressed and the two arrays are coupled only by electrostatic interaction 
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in the crosses, the system possesses a cross-sliding Luttinger liquid (CSLL) zero energy 
fixed point. 

In this Thesis we develop a theory of interacting Bose excitations (plasmons) in a 
superlattice formed by crossed interacting arrays of quantum wires. This theory goes far 
beyond the problem of stability of the CSLL fixed point. We do not confine ourselves 
with the studying the conditions under which the LL behavior is preserved in spite of 
inter-wire interaction. We consider situations where the dimensional crossover from ID 
to 2D occurs. It turns out that the standard bosonization procedure is valid in a 2D 
reciprocal space under certain conditions. The QCB behaves as a sliding Luttinger liquid 
in the infrared limit, and exhibits a rich Bose-type excitation spectrum (plasmon modes) 
arising at finite energies in 2D BZ. We derive the Hamiltonian of the QCB, analyze the 
spectrum of boson fields away from the LL fixed point and compute two-point correlation 
functions in QCB with short range inter-array capacitive interaction. We study new type 
of dimensional crossover, i.e., a geometrical crossover where the quasimomentum serves 
as a control parameter, and the excitations in a system of quantum arrays demonstrate 
either ID or 2D behavior in different parts of reciprocal space. A rather pronounced 
manifestation of this kind of dimensional crossover is related to the QCB response to 
an external ac electromagnetic field. We formulate the principles of spectroscopy for the 
QCB. We consider an infrared (IR) absorption spectroscopy of the QCB and an ultraviolet 
(UV) scattering on the QCB, and study the main characteristics of IR absorption spectra 
and UV scattering observables. 

The structure of the Thesis is as follows. In the second Chapter the progress in 
the theory of interacting fermions in low-dimensional systems (such as quantum wires, 
metallic carbon nanotubes, array of quantum wires, and QCB) exhibiting LL-like behavior 
is briefly reviewed. The bosonization procedure is introduced for a simple model of ID 
spinless interacting electrons. The LL theory is applied for describing the low-energy 
behavior of interacting electrons in real systems such as quasi one-dimensional quantum 
wires and single-wall carbon nanotubes. The existence of sliding LL phase is established 
for an array of weakly coupled parallel quantum wires. This analysis is extended to a 
system of two crossed arrays of ID quantum wires (QCB) with a capacitive inter-wire 
coupling. Such a system exhibits a crossed-sliding LL phase. We also consider QCB with 
virtual wire-to-wire electron tunneling, and find the necessary condition under which the 
one-electron tunneling is suppressed and the cross-sliding LL phase is stable. 

In the third Chapter the spectrum of boson fields and two-point correlation functions 
are analyzed in a double square QCB. We show that the standard bosonization procedure 
is valid, and that the system behaves as a sliding Luttinger liquid in the infrared limit, 
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but the high frequency spectral and correlation characteristics have either ID or 2D 
nature depending on the direction of the wave vector in the 2D elementary cell of the 
reciprocal lattice. As a result, the crossover from ID to 2D regime may be experimentally 
observed. It manifests itself as appearance of additional peaks of optical absorption, 
non-zero transverse space correlators and periodic energy transfer between arrays ( "Rabi 
oscillations" ) . 

In the fourth Chapter the effectiveness of infrared spectroscopy is studied. It is shown 
that plasmon excitations in the QCB may be involved in resonance diffraction of incident 
electromagnetic waves and in optical absorption in the IR part of the spectrum. The 
plasmon velocity is much smaller than the light velocity. Therefore, an infrared radiation 
incident on an isolated array, cannot excite plasmons at all. However in QCB geometry, 
each array serves as a diffraction lattice for its partner, giving rise to Umklapp processes 
of reciprocal super-lattices vectors. As a result, plasmons may be excited in the BZ 
center. To excite QCB plasmons with non-zero wave vectors, an additional diffraction 
lattice (DL) coplanar with the QCB can be used. Here the diffraction field contains space 
harmonics with wave vectors perpendicular to the DL that enable one to eliminate the 
wave vector mismatch and to scan plasmon spectrum within the BZ. In the general case, 
one can observe single absorption lines forming two equidistant series. However, in case 
where the wave vector of the diffraction field is oriented along some resonance directions, 
additional absorption lines appear. As a result, an equidistant series of split doublets can 
be observed in the main resonance direction (BZ diagonal). This is the central concept 
of dimensional crossover mentioned above with direction serving as a control parameter. 
In higher resonance directions, absorption lines form an alternating series of singlets and 
split doublets demonstrating new type of dimensional crossover. The latter occurs in a 
given direction with a frequency as a control parameter. 

The fifth Chapter is devoted to the study of dielectric properties of a semiconductor 
substrate with the imposed 2D QCB. We demonstrate that a capacitive contact between 
the QCB and semiconductor substrate does not destroy the Luttinger liquid character 
of the long wave QCB excitations. However, dielectric losses of a substrate surface are 
drastically modified due to diffraction processes on the QCB superlattice. QCB-substrate 
interaction results in additional Landau damping regions of the substrate plasmons. Their 
existence, form and the density of dielectric losses are strongly sensitive to the QCB 
lattice constant and the direction of the wave vector of the substrate plasmon. Thus, 
dielectric losses in the QCB-substrate system serve as a good tool for studying QCB 
spectral properties. 
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In the sixth Chapter, the principles of UV spectroscopy for QCB are formulated and 
the main characteristics of scattering spectra are described. We study inelastic scattering 
of an incident photon leading to the creation of a QCB plasmon. Experimentally, such a 
process corresponds to sharp peaks in the frequency dependence of the differential scatter- 
ing cross section. We show that the peak frequency strongly depends on the direction of 
the scattered light. As a result, the ID —>■ 2D crossover can be observed in the scattering 
spectrum. It manifests itself as a splitting of single lines into multiplets (mostly doublets). 

All technical details are contained in Appendices |XJ IEJ andEl Double tilted and triple 
QCB are considered in Appendices O and |Dl respectively. 

This work was partially presented by posters and lectures in scientific conferences and 
schools (see List of Presentations). The first part of the results was published in Refs. 1-4 
(see List of Publications). The second and third parts were published in Refs. 5-8. The 
fourth part was published in Refs. 9, 10. 

The author is grateful to V. Liubin, M. Klebanov, and Y. Imry for discussions of 
the effectiveness of the infrared absorption and ultraviolet scattering in probing spectral 
properties of QCB. 
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Chapter 2 

Prom Quantum Wires to Quantum 
Crossbars 

2.1 Introduction 

In this Chapter, a brief review of electron properties of low-dimensional systems is pre- 
sented. In Section 12.21 Luttinger-liquid (LL) theory is introduced. In Sections 12.31 and 
12.41 the LL theory is applied to describe low-energy properties of an electron liquid in a 
quasi one dimensional quantum wire and a carbon nanotube. In Section 12 .5[ an array 
of quantum wires or carbon nanotubes with local density-density and/or current-current 
inter-wire interactions is considered which result in generalized Luttinger-liquid theory. 
Similar Luttinger-liquid behavior of electron liquid in crossed arrays is considered in Sec- 
tions E12 and O 

2.2 Luttinger-Liquid Theory 

Following conventional Luttinger liquid (LL) theory [HHl EH EH1 E2] , we consider in this 
Section a simple model: a ID conductor containing spinless right- and left-moving elec- 
trons (spin sector is assumed to be gapped). In the Tomonaga-Luttinger model, the 
free-electron dispersion is assumed to be linear e^k = ±hvpk around two Fermi points 
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±fcp, and a local electron-electron interaction is parameterized by the dimensionless cou- 
pling constants g 2 and g^. The model Hamiltonian Htl — H^ in + H int is 

L/2 



-L/2 
L/2 

H int = 7rhv F J dx[g i (p 2 L (x) + p 2 R (x)) +2g 2 p L (x)p R (x)] . (2.2) 

-L/2 

Here iPl{x) (iJ)r(x)) is the field operator for left-moving (right-moving) fermions satisfying 
the anti-commutation relations {ip a {x),ipl l ,(x')} = 5 aa t5(x — x') {a, a 1 = L,R); p a (x) = 
ipl i (x)ip a (x) are the density operators for left- and right-movers. The total number of 
left (right) moving electrons is a good quantum number. Therefore all excitations are 
electron-hole-like and hence have bosonic character. 

It is convenient to write the Hamiltonian in terms of bosonic fields. The electron 
density Pl/r(x) can be expressed in terms of derivative fields d x <fL/R( x )' 

1 1 

: Pl(x) := — d x <p L {x), : p R (x) := -—d x <p R (x). (2.3) 

Here : Pl/r '■= Pl,r — (0\pl,r\0) denotes the fermion-normal-ordering with respect to the 
Fermi sea |0) defined as following [SB], to normal-order a function of operators of creation 
and annihilation of fermions, operators of creation of fermions above the Fermi level and 
operators of annihilation of fermions below the Fermi level are to be moved to the left of 
all other operators (namely operators of creation of fermions below the Fermi level and 
operators of annihilation of fermions above Fermi level). The fields ^l,r{x) satisfy the 
following commutation relations 



[<Pl/r(x), ipl/r{x')] = =R7rsign(x - x'), [^l(x), ^r(x')] = 0. (2.4) 
It is convenient to define 

= ~7t= Ivl(x) - <f R (x)) , (f)(x) = [<pl(x) + <Pr(x)) , (2.5) 

V47T V47T 



where 9{x) is a density variable and 4>(x) is the conjugate phase variable jHH]- Then one 
obtains 

L/2 

H TL = y J dx^\x)+g{d x 6{x))^j. (2.6) 

-L/2 

The Hamiltonian ()2.6)1 describes a set of harmonic oscillators, where 9{x) and ir(x) = 
d x (j)(x) satisfy the commutation relations of conventional canonically conjugate operators 



12 



of a "coordinate" and a "momentum": [6(x),tt(x')] = i8{x — x'). The renormalized Fermi 
velocity v and the LL parameter g are given by the equations 



The dimensionless parameter g is a measure of the strength of the electron-electron in- 
teractions. It plays a central role in the LL theory. The noninteracting value of g (i.e., 
for g 2 — 0) is 1, and for repulsive interactions (g 2 > 0) g is less than 1. 

2.3 Quasi One Dimensional Quantum Wire 

In this Section, the LL theory is applied to a conductor slab of length L, width R and 
thickness ro (L^-Ro^^o) containing free spinless electrons (spin sector is still assumed 
to be gapped). In experimentally realizable setups, such a structure can be created by 
cleverly gating 2D electron gas in GaAs inversion layers jlOJ EU H21 EE] , and doped helical 
polyacetilene nanofibres jH]. The position of an electron is described by a 2D vector 
r = {x, y), where the x-axis is taken along the wire direction (0 < x < L), and y is taken 
along transverse direction (0 < y < Ro). The 2D momentum of an electron is p = (p, k) 
and its dispersion is E(p) = h 2 p 2 /(2m e ) (see Fig. 12.1)1 . Here p (k) is the wave number 
in the direction of the wire axes (in the transverse direction), m e is an effective electron 
mass. Up to scales \p\ < po = n/Ro and E < E = f?p\j (2m e ), all the excitations are 
one dimensional. In the subsequent discussion we assume that pf < Po and Ep < E . 
Then we have "left"- and "right" -moving quasi-particles with energies near the Fermi 
level and momenta near pp (—Pf) for right-moving (left-moving) fermions. We introduce 
the momentum index k = p — pp (k = p + pp) for right-moving (left-moving) fermions. 
It is seen that — pf < k < oo (— oo < k < pp) for right-moving (left-moving) fermions. 

Following Ref. [HE], we extend the range of k to be unbounded by introducing addi- 
tional unphysical "positron states" at the bottom of the Fermi sea. Next, we factor out 
the rapidly fluctuating e ±tppX phase factors and express the physical fermion field *&(x) 
in terms of two fields iPl/r{x) that vary slowly on the scale of 1/pf'- 



The energies near the Fermi level can be written as Ek = Ep + fwpk [E^ = Ep — hvpk) 
for right-moving (left-moving) fermions, where Vp = hpp/m e is Fermi velocity. Then in 
this approximation the kinetic energy Hamiltonian has the form of the Hamiltonian (|2.1j) . 
Electron-electron interaction is a Coulomb interaction screened in the long-wave limit 
Indeed, quantum wires are not pure ID objects and screening arise due to their 




(2.7) 



= e ipFX tfj R (x) + e 



y> L (x). 



(2.8) 
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finite transverse size Rq which is the characteristic screening length |46j. It is described 
by the following Hamiltonian 

dxdy f dx'dy' 



int. 



Rq 



Rq 



-U(t - r')p(x,y)p(x',y'). 



(2.9) 



Here 



e 2 C 



U(t) 



x 

R 



[2.10) 



where tq is the stripe thickness, the screening function (introduced phenomeno logi- 
cally) is of order unity for |£| < 1 and vanishes for |£| > 1. In the long- wavelength limit, 
the electron density operator p(x,y) = p(x) = ty'(x)ty(x) can be written as 



p(x) = p L (x) + p R {x) + i)\{x)il) R {x)e lipFX + i) R (x)i) L (x)e 



-2ippx 



(2.11) 



where p a (x) = ip^x^^x) (a = L,R) are density operators for left- and right-moving 
fermions. 
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Figure 2.1: The low energy part of the free-electron spectrum. Here po = tt/Rq and 
E = h 2 pQ/(2m e ). k is measured from ^op {—Pf) for right-moving (left-moving) fermions. 
Tilted lines describe the linear approximation Ep ± hvpk for right/left moving electrons. 



Then one can write the interaction ()2.9|) in the form (J2.2j) . where 

1 



1 fdxdyUir) 2e 2 
TV J Rq Flvf nvp 



_ j ^y^M. (1 _ cos (2p F x)) « ^— (p F R Q ) 2 
71 J Rq Kvf 3Hvf 



With Eqs. (J2~3J) and IjESj) . the Hamiltonian H = H kin + H mt acquires the form (|2.fi|) . 
where renormalized Fermi velocity v and the dimensionless interaction parameter g are 
given by Eqs. (|2.7jl . In the GaAs slab with a density of one electron per 10 nm and the 
width Rq ~ 1 nm, m e = 0.068mo (mo is the free electron mass), Vp ~ 10 7 cm/sec and 
then g ~ 0.97. 
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2.4 Luttinger Liquid Behavior in Single- Wall Carbon 
Nanotubes (SWCNT) 

Nanotubes are tubular nanoscale objects which can be thought as a graphite sheet wrapped 
into a cylinder ^3 HHj- The arrangement of carbon atoms on the tube surface is deter- 
mined by the integer indices 0<m<n of the wrapping superlattice vector T = nai + m&2 
|491 1571] , where &i and a 2 are the primitive Bravais translation vectors of the honeycomb 
lattice (see Fig. 12.2)1 . The first Brillouin zone of the honeycomb lattice is a hexagon 
(Fig. EHJ), and there are two independent Fermi points denoted by K and — K, with two 
linearly dispersing bands around each of these points. The necessary condition of metal- 
licity of SWCNT is 2n + m = 31 for an integer / jlH] • If this condition is not fulfilled, 
the nanotube exhibits the band gap AE ~ 2hvp/(3Ro) ~ 1 eV |4U) 151). where Vp is the 
Fermi velocity and Rq is the nanotube radius. Even if the necessary condition is fulfilled, 
the rearrangement of local bonds due to the curvature of the cylinder can introduce a 
gap, AE ~ 10 meV, which implies narrow-gap semiconducting behavior. For very small 
tube diameter (1 nm or less), due to the strong curvature-induced hybridization of o and 
7r orbitals, this effect can be quite pronounced In the cases of "armchair" (n = m) 
and "zigzag" (m = 0) nanotubes, however, the formation of a secondary gap is prevented 
by the high symmetry, and therefore armchair and zigzag nanotubes are metallic 

In this Section, a single metallic SWCNT is under consideration. The effective low- 
energy description of SWCNT is derived. Coulomb interaction between electrons induce 
a breakdown of Fermi liquid theory. It is shown that the bosonization procedure (similar 
to the procedure derived in Section f2.1)l is valid. As a result, interacting electrons in a 
metallic SWCNT exhibit Luttinger liquid behavior. 

The electronic properties of carbon nanotubes are due to special band-structure of the 
7r-electrons in graphite [54J. The band structure exhibits two Fermi points kK (k = ±) 
with a right- and left-moving (a = R/L) branch around each Fermi point. These branches 
are highly linear with Fermi velocity v p ~ 8 x 10 7 cm/s. The R- and L- movers arise as 
linear combinations of the r = A,B sublattice states reflecting the two C atoms in the 
basis of the honeycomb lattice. The dispersion relation is linear for energy scale E < D, 
with the bandwidth cutoff scale D « hvp/Ro for tube radius _R . For typical SWCNT, D 
is of the order 1 eV. The large overall energy scale together with the structural stability 
of SWCNTs explain their unique potential for revealing Luttinger liquid (LL) physics. 
The fermionic quasi-particles in the vicinity of the Fermi level of graphite are described 
by the 2D massless Dirac Hamiltonian This result can also be derived in terms of 
k • p theory [5T1] • 
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(a) 

Figure 2.2: Honeycomb lattice of 2D carbon sheet and the coordinate system. Hexagonal 
sublattices A and B are labelled by • and o, respectively. Two primitive translation 
vectors are a.i and a 2 , where |ai| = |a 2 | = a = 2.5 A j^. The vector directed from an 
A site to a nearest neighbor B site is d^. A nanotube is specified by a chiral vector T 
corresponding to the circumference of the nanotube whereas the x-axes is oriented along 
the nanotube axes. 

Wrapping the graphite sheet onto a cylinder then leads to the generic band-structure 
of a metallic SWCNT shown in Fig. 12.41 Quantization of transverse (azimuthal) motion 
now allows for a contribution oc exp(imy/R ) to the wave function. Here the x-axis 
is taken along the tube direction, the circumferential variable is < y < 2irR (y is 
really the azimuthal angle multiplied by the nanotube radius Ro). However, excitation 
of angular momentum states other than m = costs a huge energy of order D « 1 eV. 
In an effective low-energy theory, we may thus omit all transport bands except m = 
(assuming that the SWCNT is not excessively doped). Evidently, the nanotube forms 
a ID quantum wire with only two transport bands intersecting the Fermi energy. This 
strict one- dimensionality is fulfilled up to a remarkably high energy scales (eV) here, in 
contrast to conventional ID conductors. The Hamiltonian of kinetic energy is: 

H kin = -ihv F ^ dx (4>\ Ka (x)d x 4> BKa (x) + i) ] BK(T {x)d x i) AKC7 (x)^ , (2.12) 

where ip rKC r{x) is the "smooth" field operator of fermions in the sublattice r = A, B in 
the vicinity of the Fermi point kK (k = ±) with spin a =T,i |51j . Switching from the 
sublattice (r = A, B) description to the right- and left-movers (a = R,L), 

i)R/L K o{x) = -j= (^A K a(x) =F t/WW) , 
implies two copies (k = ±) of massless ID Dirac Hamiltonians (similar to the Hamiltonian 
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(b) 

Figure 2.3: First BZ of 2D carbon sheet. Here o and • are two vertices of the BZ 
corresponding to the vectors K and — K respectively, |K| = 47r/(y / 3ao)- Wave vectors of 
low-energy quasi-particles lie at the axes k in the vicinity of the points ±K. 

(J2HJ) for each spin direction, 

recr 

The electron-electron interaction is the screened Coulomb interaction between charge 
fluctuations 46^ The re-distribution of a charge induced by "external" charge can be 
described by the envelope function (introduced phenomenologically) C(0> £ = x /Ro, 
= C( — 0) C(0) ~ 1- This function is of order of unity for |£| ~ 1 and vanishes outside 
this region. Thus similarly to Eq. (j2.1(J|) . the interaction is introduced as 



U(v) 



( x 

\R 



^ + 4^ sin 2 QjA+rl 



1 1/2 ' 



(2.13) 



where r# denotes the average distance between a 2p z electron and the nucleus, i.e., the 
"thickness" of the graphite sheet. In the long- wavelength limit, electron-electron interac- 
tions are then described by the second-quantized Hamiltonian [3T] 

Hl = \ 5l/ ^ dx '^}( x_ ^)^i»i<T( K )^ l c 2 a'(^)V'T'i e 3^(a ; ')V'r«4<T(a:), (2-14) 

tt'cjcj 1 {re,:} 

with the ID interaction potentials 

V { %(X - X') = j ei K((re 1 -re 4 ) r+ (re 2 -rea)O f/(r . _ ^ + ^ (2 . 15) 

These potentials depend only on x — x' and on the ID fermion quantum numbers. For 
interactions involving different sublattices t ^ t' for r and r' in Eq. ([2.15J1 . one needs 
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Figure 2.4: The low-energy part of the band structure of a metallic nanotube in the 
vicinity of the vertices K and K' of the BZ (see Fig. 12.3)1 . Lines 1,2 (3,4) correspond to 
excitations with orbital quantum number m = (m = ±1). The bandwidth cutoff scale 
is u = Vp/R . The Fermi level E F (and momentum q F = Ep/vp) can be tuned by an 
external gate. 

to take into account the shift vector d pp i (d^B = ~ d.BA, I^ab] — 1-44 A jl^l) between 
sublattices (see Fig. 12. 2j) . 

To simplify the resulting ID interaction ()2.14j) . we now exploit momentum conserva- 
tion, assuming Ep ^ (see Fig. I2.4j) so that Umklapp electron-electron scattering can be 
ignored, and only the processes conserving the number of electrons for each channel tk 
will be considered, i.e., K\ + k 2 = k 3 + k 4 . We then have "forward scattering" processes 
and "exchange interaction", where Ki = and K2 = k^. In addition, "backscattering" 
processes may be important, where K\ = —k 2 = k% = —K4. 

In the next step, we employ the fact that the potentials V{ Ki }(x — x') are screened 
with the radius of screening being of order of the nanotube radius Rq, whereas the field 
operators ^'s are slowly varying on this distance scale. As a result, one can approximate 
the short-range potentials V{ Ki y(x — x') by delta-like potentials. Then, switching to the 
right- and left-mover representation, the Hamiltonian ()2.14|) can be written in the form 



Here Hf describes the "forward scattering" processes, H x is the exchange interaction, and 




Hi 



(2.16) 
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H;, corresponds to "backscattering" , 

Pouut'{x) = ^1pl KCT {x)lp aK '*{x), Pan{x) = Pa KK (x), p{x) = ^p aK (x), 

(7 OLK 

S aKK >(x) =^2lpl Ka (x)T (7a <iJ aK , a/ (x), 

T aa i is the vector of Pauli matrices. The effective coupling constants 70, 71, and 75 are 
given by the equations: 

72„ „ r 



e 2 7o 



where d^P is the x-component of the vector d pp i, q F = Ep/ihvp)- For carbon nanotube, 
d AB = lM ^ qp w 1/(3^), R « 4 A [18J, AT = 4^/(00^), a = 2.5 A [4QJ, then 
70 = 1.3, 71 « 0.09, 7 fc « 0.08, i.e., 7^76 < 70. 

The electron density p QK(J can be expressed in terms of the derivative fields d x ip aKcr 
similar to (|2.3|) 

: PL K a(x) := ^-d x (p LKa (x) , : p RK(T {x) := -^-d x (p Bjia {x) . 

Here : p aKtJ := Pq, K(J — (0|p areCT |0) denotes fermion normal ordering with respect to the 
Fermi sea |0) j3H]. The fields (f aKa (x) satisfy the following commutation relations 

[<fL/R Ka(x),Lp L/R K ' a '{x')] = ^ilt5 KK > 5 acT 'Sign(x - X'), [<PLK*(x), ( PR K 'a'{%')) = 0. 

It is natural to introduce the standard linear combinations 6\ u (x) and their dual fields 
<Pxv( x ) subject to the algebra 

[6\u{%), <t>\'v'{x')] = --5 AA '5 w 'sign(x - x') . (2.20) 

The bosonic density fields 0\ v {x) for the total [y = g) and relative [y = u) charge (A = c) 
and spin (A = s) channels are constructed as 

0c,g/u = -j— f= - <Pr+~\ ± <Pl-t =F <Pr-i + - <PR.+[ ± fL~l =F PR-[] ) 

0s,g/u = -7^1= [<Pl+T - <PR+\ ± ^l-T =F <PR-1 - <Pl+i + Vr+i =F (Pl-i ± (pR-i] ■ 
4V 71 " 

Their dual phase fields (fi\ u are defined similarly 
1 



} c,g/u 



40F 
1 

4v^ 



[<Pl+T + V^R+T ± ¥>l-t ± <Pfl-t + V^L+I + <PR+| ± ¥>l-.i ± <Pr-i] 



^s, g /u = ~rj= [<Pl+t + <PR+T ± ¥>l-T ± ^K-T _ V?l+| ~ <Pfl+i =F <Pl-| =F <Pk-J 



19 



The Hamiltonian Hf can be written purely in terms of charge bosonic field operators 

O c ,g/u(x) and tt C;9/u (x) = d x (f) Cig/u (x), 



dx { (2 7o - |) {d x d cg {x)f - | [*i(aO + 7r c 2 u (x) - (5^ cu (x)) 2 ] 



The Hamiltonian H x leads to nonlinearities in the 9 CU charge field and the S)9 / U spin fields. 
The four channels are obtained by combining charge and spin degrees of freedom as well 
as symmetric and antisymmetric linear combinations of the two Fermi points, K = ±. 
The bosonized expression for H x reads 



H, r 



47i 
L 2 

+ cos ( VAtt 6 cu {x) ) cos ( VAtt 9 sg (x) ) — cos ( VAtt 9 sa (x) ) cos ( VAtt 8 su (x 



/^:[cos(^Mx))oos(^Mx)) + 



(2.21) 



'sg\ 



Here : . . . : denotes boson-normal-ordering defined as follows: to boson-normal-order a 
function of operators of creation and annihilation of bosons, all creation operators are to 
be moved to the left of all annihilation operators. 

Similar to H x , the backscattering Hamiltonian leads to nonlinearities in the 8 CU and 
9 SU fields. The bosonized expression for Hi, takes the form [HH] 



H b = jj J dx : [cos (VA^6 cu {xj) cos (V4^6 m {xj) + (2.22) 
+ cos ( VAtt 9 cu (x) ) cos ( VAtt <p su (x) \ + cos ( VAtt 9 su (x) ) cos ( VAtt 4> su (x) 



Writing the non-interacting Hamiltonian H (j2.12j) in terms of bosonic field operators, 
one obtains 



H = ^ H Xv + H x + H b , H Xu = I 



dx 



9\uirl v {x) + — {d x 6 Xl/ (x)f 
9\u 



(2.23) 



where A = c, s; v — g,u] H x is given by Eq. (j2.21|) . H b is given by Eq. ()2.22|) . 



VF 

9\u 



Here 



9cg = 9 



9i 



l-9i 



l + 9o + 9i 



1/2 



0.25, « 0,, fl /u « 1. 



7i e 



9o 



4 7o e 2 



Clearly, the charged (eg) mode propagates with significantly higher velocity than the three 
neutral modes. There is a further renormalization of the values v cu and v St9 / u , however, 
this effect is very small and can be neglected. Renormalization group analysis [ST] shows 
that the contribution H x is marginally irrelevant, whereas the backscattering part H b is 
marginally relevant. 
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There are several possible regimes in which a nanotube can exist 0]. First, there is 
an insulating regime with the density at half filling (qp = 0), where all excitations are 
gapped. Second, there are conducting states which can be realized by applying various 
external fields. These fields may close some gaps or even all of them, provided their 
magnitudes exceed certain critical values. For example, by varying the chemical potential 
(or changing qp), one can close all the gaps. This leads to a transition into a metallic 
regime. 

2.5 Sliding Luttinger Liquid Phase 

The simplest 2D ensemble of ID nanoobjects is an array of parallel quantum wires or 
nanotubes. The inter-wire interaction may transform the LL state existing in isolated 
quantum wires into various phases of 2D quantum liquid. However, the density- density 
or/and current-current inter- wire interactions do not modify the low-energy behavior of 
quantum arrays under certain conditions. In particular, it was shown recently [22] that 
"vertical" interaction which depends only on the distance between the wires, imparts the 
properties of a sliding phase to 2D array of ID quantum wires. 

In this Section, a 2D array of coupled ID quantum wires is considered and the question 
of existence of a stable 2D phase that retains some of the properties of ID Luttinger liquid 
is addressed. Following Ref. [22] this phase will be called as sliding Luttinger liquid (SLL). 
An anisotropic 2D system composed of parallel chains with spinless Luttinger liquids (LL) 
in each chain (LL in the spin-gapped phase) is considered. It will be shown that the 
long- wavelength density- density and/or current- current interactions between neighboring 
Luttinger liquids are marginal operators which result in the sliding Luttinger-liquid phase. 

Let us consider an array with N chains, each labelled by an integer n<i = 1, 2, . . . , N. 
The conventional LL regime in a single ID quantum wire is characterized by bosonic 
fields describing charge modes (LL in the spin-gapped phase). It is assumed that all wires 
of the array are identical. They have the same length L, Fermi velocity v and Luttinger 
liquid parameter g. The period of the array is a. The axis x\ is chosen along the array 
direction, whereas the x 2 axis is perpendicular to the array. The excitation motion in 
QCB is one-dimensional in major part of the 2D plane. The anisotropy in real space 
imposes restrictions on the possible values of the coordinate x 2 . It should be an integer 
multiple of the array period a, so that the vector r = (xi,n 2 a) characterizes the point 
with the ID coordinate x\ lying at the ri2-th wire of the array. The low-energy Luttinger 
liquid of each wire with spinless interacting fermions is described by the Hamiltonian 
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(|2.6|) . The Hamiltonian of the array without inter- wire interaction reads 

L/2 

H °i=Y^ J ^i[9^ 2 (^n 2 a) + ^(d xl 9(x 1 ,n 2 a)) 2 Y (2.24) 

™ 2 -L/2 

where 7r are the conventional canonically conjugate boson fields. 

The interactions between the chains correspond to couplings between the long wave- 
length components of the densities p(xi,n 2 a) and of the currents j(xi,n 2 a) [23J: 

L/2 

H in t = — ^2 / dxi Wj(n 2 - n 2 )j(xi,n 2 a)j(x 1 ,n 2 a) + (2.25) 

+W p (n 2 - n' 2 )p(x 1 ,n 2 a)p(x 1 , n' 2 a) 

With Eqs. (j2.3|) and ()2.5|) . the electron density p = Pl + Pr an d current j = Pl — Pr 
can be expressed in terms of d Xl 8 and tt respectively. Then the bosonized form of the 
Hamiltonian Hi = H® + for the interacting liquids has the form 

L/2 

Hi = — / dxi{Kj(n 2 - n' 2 )ir(xi ,n 2 a) n(xi,n' 2 a) + 

n 2 n' 2 _ L/2 

+K p (n 2 -n' 2 )(d xi 9(xi,n 2 a))(d xi e(x 1 ,n 2 a))}, 
where the coupling matrices are 

Kj(n 2 - n' 2 ) = {g5 mn > 2 + Wj(n 2 - n 2 )) , K p (n 2 - n' 2 ) = (^^l + W p (n 2 - n' 2 )^j . 

The Hamiltonian Hi describes coupled harmonic oscillators and can be diagonalized. By 
introducing Fourier transforms in the direction perpendicular to the wires, 

9(xi,n 2 a) = -^=J2 el92n2a9 ^i), 7i(xi,n 2 a) = -±= £ e^ir^), 



N — ' ** VN 

<?2 <?2 



(\q 2 \ < Q/2, Q = 2n/a, 6^ = 0^ q2 , 7rJ 2 = n^ q2 ) the Hamiltonian Hi can be rewritten in 
the form similar to the Hamiltonian ()2.24|) : 



L/2 



#i = E^T / + (2-26) 

92 -L/2 

where the Luttinger liquid parameters g q2 and the velocities v q2 are defined as 



v q2 = vJKj(q 2 )K p (q 2 ), g q2 = J K J/p (q 2 ) = £ K J/p (n 2 )c^ 



KM) 
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The Hamiltonian ()2.25|) is invariant under the transformations 4>{x\, n 2 a) — > 4>{x\, n 2 a) + 
C n2 and 9(xi,n 2 a) — > 6 , (xi,n 2 a) +-D n2 , where 7r(xi,n 2 a) = d xi (f)(xi, n 2 a), C n2 and L> n2 
are constants on each wire. The corresponding phase is called as a SLL one [22J. In 
this phase, the total numbers of left (right) moving electrons on each chain are good 
quantum numbers and expectation values (ip a (xi, n 2 a; t)ijr al {x'i, n 2 a; 0)) (a, a' = L, R) for 
n 2 ^ n' 2 are necessarily zero in this phase. This corresponds to a perfect charge insulator 
in the transverse direction. Density correlations in the transverse direction are, however, 
nontrivial. For short ranged density-density and current-current interactions, they decay 
exponentially with increasing distance between the wires. The low energy modes are sets 
of ID density oscillations (sound) propagating along each wire of the array with a wave 
number k\ and a phase shift q 2 a between adjacent wires. The dispersion of the modes is 
linear with respect to ID wave number ki, E(ki,q 2 ) = v q2 \ki\. These modes can, for in- 
stance, transport heat perpendicular to the chains although the system is a perfect charge 
insulator in this direction. 

2.6 Cross-Sliding Luttinger Liquid Phase 

Next, a square network of ID wires formed by coupling two perpendicular arrays of 
chains is considered. In experimentally realizable setups jH] these are cross-structures of 
suspended single-wall carbon nanotubes placed in two parallel planes separated by an 
inter-plane distance d. However, some generic properties of QCB may be described in 
assumption that QCB is a genuine 2D system. The system consists of two periodically 
crossed arrays of ID quantum wires. It is assumed that all wires are identical. They have 
the same length L, Fermi velocity v and Luttinger parameter g. A coordinate system 
is chosen so that the axes Xj and the corresponding basic unit vectors e,,- are oriented 
along the j-th array (j = 1,2 is the array index). The period of crossbars is a, and 
the basic vectors are a, = aej (Fig. 12.5)1 . The interaction between the excitations in 
different wires includes both interaction between wires from the same array (intra-array 
interaction) and wires from different arrays (inter-array interaction) . The former is given 
by Eq. (|2.25jl . the latter is assumed to be concentrated around the crossing points with 
coordinates niai + n 2 a 2 = {n\a, n 2 a). The integers rij enumerate the wires within the 
j-th array. Following Refs. [2U I2H] h will be shown that it exhibits a new crossed sliding 
Luttinger liquid (CSLL) phase. 

The LL of the first array is described by the Hamiltonian Hi (I2.2(i|l . The Hamiltonian 
H 2 of the second array is obtained from Hi after permutation 1 <-» 2 in the arguments of 
the fields. The density-density interactions between electrons on intersecting wires gives 
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Figure 2.5: 2D crossbars formed by two interacting arrays of parallel quantum wires. 
Here ei,e 2 are the unit vectors of the superlattice, a\ = a 2 =a is the superlattice period 
(the case ai^a>2 is considered in Appendix EJ) and d is the vertical inter-array distance. 

rise to a term in the Hamiltonian in the form 

L/2 

H12 = / dxidx2V(xi — nia, X2 — n2a)pi(xi,n 2 a)p2(nia, X2), (2.27) 

mn 2 _ L/2 

where V(x\ — nia, X2 — n 2 a) is a short-range inter-array interaction. Using the definitions 
Pi(xi,ri2a) = d Xl 9\(xi, n 2 a)/ v /7r and p2(^ia,x 2 ) = d X2 82(nia, x^) / y/rr, one obtains the 
bosonized form of the Hamiltonian H12 [2*3*] : 

L/2 

H\2 = — yj / dx\dx2V(xi — riia,X2 — n,2a)d Xl 9i(xi,n2a)d X2 92(nia,X2).(2.28) 

71 nin2 -L/2 

Let us introduce Fourier transforms according to jHEj 

6 l {x u n 2 a) = -^^e i(9l+miQ)a:i+4 ' ?2n2 ^ lq+mi , (2.29) 

v N L 

qm 2 

where wave vector q belongs to the first Brillouin zone, {qx^l < Q/2, m i, 2 — m i,2Q e i,2 are 
reciprocal superlattice vectors ""HI U23] (Q = 2n /a is reciprocal superlattice constant and 
mi^ are integers). Then removing degrees of freedom with wavelengthes smaller than 2a 
(that is, considering only the modes with my = 0), one obtains the total Hamiltonian 

H = H 1 + H 2 + H 12 M 

- - N - ) 



H = 2^ I — 2 — lq q — 2 — 2q q J 2flLW^' ^ 2 ' 31 



q «jq 
where z, j = 1,2 denote the array number, 



^ = QiQj 



f^lJ- Sjl + ^^-5 j2 ) *y + ^(q)(l - 5 tJ ) 
\v g q2 J v g qi J 



L/2 

V(q) = — — [ dx 1 dx 2 V(x 1 ,X2)e iqiXl+iq2X2 . 
ixhva J 



-L/2 
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The Hamiltonian ()2.31|) exhibits a CSLL phase. Renormalization group analysis [2*3] 
shows that additional interactions between the two arrays, such as the Josephson and 
charge- density- wave couplings, are irrelevant and the CSLL phase is stable. 

It should be noted that the Hamiltonian (j2.31j) does not include interactions of the 
long- wavelength plasmons with quasi-particles whose wavelengthes smaller than 2a, i.e., 
plasmons in higher energy bands. Being weak, these interactions renormalize the plasmon 
velocity (see Eq. (|B.1J) ). When the interaction strength increases, the lowest QCB modes 
soften and their frequencies vanish in a whole BZ at a certain critical of the interaction 
strength (see Section l3~"""p . 

2.7 Quantum Crossbars with Virtual Wire-to- Wire 
Tunneling 

To finalize the substantiation of the CSLL family, consider the condition when the tun- 
neling of electrons in quantum crossbars is suppressed and the Luttinger-liquid-like phase 
is stable. For the case when the charge degrees of freedom are quenched, we derive the 
effective spin Hamiltonian of QCB. 

Let us consider two non-parallel metallic nanotubes. One of them belongs to the 
first array and has the number n 2 , and another one belongs to the second array and 
has the number n\. The intersecting point has the coordinates (ra 1 a,n 2 a). The elec- 
tron hopping between the wires gives rise to a term in the Hamiltonian of the form: 
S CT tl^ia( n i a i n 2 a)^/ 2a (riia, n 2 a) + h.c], where t is an effective hopping constant, and 
a =|, | is a spin index. Next we introduce "slowly varying" field operators ipL a and ipRa 
(similar to (|2.8j) ). Assuming that the Fermi vector kp is not commensurate with the 
reciprocal superlattice vector Q = 2n/a, we can write k F = m F Q + q , with m F being 
integer and |go| < Q/2. Then we represent the hopping between the two arrays as 



where q F = qo and qn = —go- Here and below we assume that t is real and positive. 

The energy cost of the electron wire-to-wire tunneling is the energy 2Eq necessary to 
charge both wires, 



where L is the nanotube length and Rq is the nanotube radius. The tunneling is suppressed 
if t Ec- The tunneling constant t can be estimated from the transport experiment 




i){ aa (nxa, n 2 a)ilj2a'a(n 1 a, n 2 a) + h.c. 



(2.32) 




(2.33) 



25 



through crossed nanotubes. The elastic and van der-Waals interactions between crossed 
nanotubes determine two equilibrium positions [0] with inter-wire distances d = 1 nm and 
2 nm. It is shown that for d = 2 nm, the resistance is R ~ 10 10 fl On the other hand, 
resistance can be evaluated from the well known Landauer formula 

1 „ 2e 2 z/ t 2 L 

where uq is the density of states in the quantum wires. For real nanotubes, the length of 
a ballistic transport is L exp ~ lpm and E c ~ 20 meV [THJ- Taking v = 8 • 10 7 cm/sec, 
Rq = 0.4 nm, and a = 20 nm, we have t/a ~ 1 peV, i.e., t <C E'c. Then, single-particle 
hopping between nanotubes is suppressed, and only backward scattering interaction can 
be relevant. However, if kp is not commensurate with the reciprocal lattice constant Q, 
the back scattering is also suppressed. 

In this case the effective interaction Hamiltonian can be obtained from the initial one 
()2.32j) by means of elimination the states from adjacent charge sectors by using the second 
order perturbation theory with respect to t/(aEc) <C 1: 

t 2 

Hint = tttt - ^ {pi(n 1 a,n2a)p2(n 1 a,n 2 a) + 2Si(nia, n 2 a) ■ S 2 (n 1 a,n 2 a)}. (2.35) 

nvni 

Here pj is the electron density operator and Sj is the spin operator, 



aaa' 



where r is the vector of Pauli matrices. 

The first term under the sum in the right-hand side of Eq. ()2.35|) describes the potential 
interaction between charge fluctuations and simply renormalizes the coupling strength 
V(xi,x 2 ) of the Coulomb inter-wire interaction (|2.27jl . The second term under the sum 
in the right-hand side of Eq. (|2.35jl is spin-spin inter- wire interaction. 

The Hamiltonian ()2.35j) shows that the virtual wire to wire hopping results in only a 
slight renormalization of the strength of the capacitive inter-wire interaction. However, 
the virtual hoping results in the effective spin-spin interaction. 

Short summary and future outlook: It was demonstrated in this Chapter that 
the LL fixed point is conserved in a quasi 2D periodic network (quantum crossbars). In 
the following chapters we concentrate on another aspect of the problem: we consider the 
properties of QCB at finite wave vectors q and frequencies u in the whole 2D Brillouin 
zone, study the spectrum of Bose excitations u(qi,q 2 ), consider situations where dimen- 
sional crossover ID <— > 2D occurs, formulate the principles of infrared and ultraviolet 
spectroscopy for QCB and study the basic characteristics of the spectra. 
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Chapter 3 



Plasmon Excitations and One to 
Two Dimensional Crossover in 
Quantum Crossbars 

In this Chapter we analyze the spectrum of a specific nano-object - square double quantum 
crossbars (QCB). We show that the standard bosonization procedure is valid and the 
system behaves as cross-sliding Luttinger liquid in the infrared limit. However plasmon 
excitations in QCB demonstrate both ID and 2D behavior depending on the direction of 
the plasmon wave vector. We discuss several crossover effects such as appearance of non- 
zero transverse space correlators and the periodic energy transfer between arrays ("Rabi 
oscillations"). The spectra of a tilted double and triple QCB are considered in Appendices 
Oand|D] 

3.1 Introduction 

A double 2D grid, i.e., two superimposed crossing arrays of parallel conducting quantum 
wires or nanotubes, represents a specific nano-object which is called quantum crossbars 
(QCB). Its spectral properties cannot be treated in terms of purely ID or 2D electron 
liquid theory. A constituent element of QCB (quantum wire or nanotube) possesses 
the LL-like spectrum ^3 EE] • A single array of parallel quantum wires is still a LL-like 
system qualified as a sliding phase 16 provided only the electrostatic interaction between 
adjacent wires is taken into account. If tunnelling is suppressed and the two arrays are 
coupled only by electrostatic interaction in the crosses, the system possesses the LL zero 
energy fixed point. 



27 



In the present Chapter we concentrate on another aspect of the problem of capac- 
itively interacting arrays of quantum wires. Instead of studying the conditions under 
which the LL behavior is preserved in spite of inter-wire interaction, we show that a rich 
Bose-type excitation spectrum (plasmon modes) arises at finite energies in the 2D Bril- 
louin zone (BZ) [SEl EZj and consider situations where the dimensional crossover from 
ID to 2D occurs [SHJ EH|- We start our studies of QCB with a double square QCB, 
namely, m = 2. In the first two Sections 13.21 and 13.31 we introduce basic notations and 
construct the Hamiltonian of the QCB. The main approximations are discussed in Sec- 
tion 13.41 Here we substantiate the method used (separable interaction approximation) 
and show that interaction between arrays in QCB is weak. The energy spectra for square 
QCB are described in detail in Section 13.51 Various correlation functions and related 
experimentally observable quantities (optical absorption, space correlators) are discussed 
in the last Section l3~o1 We predict here the effect of peculiar "Rabi oscillations" - periodic 
energy transfer from one of the QCB array to another. All technical details are placed in 
Appendices |X] and El 

Double tilted QCB is considered in Appendix O Triple QCB (m = 3) formed by three 
arrays lying in parallel planes are studied in Appendix [Dj Such hexagonal grids may be 
useful for three-terminal nanoelectronic devices jSJ. The plasmon spectra of triple QCB 
possess some specific features in comparison with double QCB. We introduce the main 
notions and construct the Hamiltonian of symmetric triple QCB ( section lD.1 J) , analyze the 
peculiarities of the frequency spectrum ( section \D.2\i . and illustrate them by description 
of triple Rabi oscillations - periodic energy transfer between all three arrays (part |D~3|) . 
The results are summarized in the Conclusions. 

3.2 Basic Notions 

Double square QCB is a 2D periodic grid, which is formed by two periodically crossed 
arrays of ID quantum wires. In experimentally realizable setups [Hj these are cross- 
structures of suspended single-wall carbon nanotubes placed in two parallel planes sep- 
arated by an inter-plane distance d. However, some generic properties of QCB may be 
described in the assumption that QCB is a genuine 2D system. We assume that all wires 
of QCB are identical. They have the same length L, Fermi velocity v and Luttinger pa- 
rameter g. The periods of a crossbars is a. The arrays are oriented along the unit vectors 
ei t 2, and corresponding basic vectors are a, = ae,- (Fig. 12. 5jl . 

The interaction between the excitations in different wires is assumed to be concentrated 
around the crossing points with coordinates nis.\ + n 2 a 2 = (nia, 712a). The integers rij 
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enumerate the wires within the j-th array. Such an interaction imposes a super-periodicity 
on the energy spectrum of initially one dimensional quantum wires, and the eigenstates 
of this superlattice are characterized by a 2D quasimomentum q = qxe\ + q^&2 = (<?i, (fe)- 
The corresponding basic vectors of the reciprocal superlattice have the form m = mi+m2, 
where m 12 = mi,2<5 e i,2, Q = 27r/a and are integers. 

However, the crossbars kinematics differs from that of a standard 2D periodic sys- 
tem. In conventional 2D systems, forbidden states in the inverse space arise due to Bragg 
diffraction in a 2D periodic potential, whereas the whole plane is allowed for wave propa- 
gation in real space, at least when the periodic potential is weak enough. A set of Bragg 
lines correspond to different reciprocal lattice vectors m. A Brillouin zone is bounded 
by the Bragg lines. It coincides with a Wigner-Seitz cell of reciprocal lattice (Fig. 13.1k ). 
In contrast, most of the real space is forbidden for electron and plasmon propagation in 
sharply anisotropic QCB. The Bragg conditions for the wave vectors of a quasi-particle 
in an array are modulated by a periodic potential created by another array, unlike those 
in conventional 2D plane. There are two sets of Bragg lines corresponding to recipro- 
cal lattice vectors mi and ni2. These conditions are essentially one-dimensional. The 
corresponding BZ is a Wigner-Seitz cell of a reciprocal lattice shown in Fig. 13.1b . 



q 2 /Q 




1 1.5 qi/Q 



(a) (b) 

Figure 3.1: The Bragg lines in conventional 2D square lattice (a) and that of 2D square 
crossbars (b). In both panels, solid lines are the Bragg lines, whereas points correspond 
to reciprocal superlattice vectors. 

Indeed, the excitation motion in QCB is one-dimensional in major part of the 2D plane. 
The anisotropy in real space imposes restrictions on the possible values of 2D coordinates 
x\, X2 (r = x\e\ + X2e 2 ). At least one of them, e.g., X2 (xi) should be an integer multiple 
of the corresponding array period a, so that the vector r = (xi,ri2a) (r = (nia, X2)) 
characterizes the point with the ID coordinate x\ (a^) lying at the 712-th (rii-th) wire of the 
first (second) array. As a result, one cannot resort to the standard basis of 2D plane waves 
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when constructing the eigenstate with a given wave vector k. Even in non-interacting 
arrays of quantum wires (empty superlattice) the 2D basis is formed as a superposition 
of two sets of ID waves. The first of them is a set of ID excitations propagating along 
each wire of the first array characterized by a unit vector fciei with a phase shift ak 2 
between adjacent wires. The second set is the similar manifold of excitations propagating 
along the wires of the second array with the wave vector k 2 e 2 and the phase shift ak\. The 
dispersion law of these excitations has the form u°(k) = uj{k\)-\-uj{k 2 )- The states of equal 
energy obtained by means of this procedure form straight lines in the 2D reciprocal space. 
For example, the Fermi surface of QCB developed from the points ±kpi,2 f° r individual 
quantum wire consists of two sets of lines |fci )2 | = &fi,2- Respectively, the Fermi sea is not 
a circle with radius &p like in the case of free 2D gas, but a cross in the k plane bounded 
by these four lines [THj (see Fig. 13. 2j) . Finally, the Bragg conditions read 

uj(ki) - uo(ki + niiQ) + uo(k 2 ) - uo(k 2 + m 2 Q) = 0. 

and the lines \k\\ = Q/2 and \k 2 \ = Q/2 satisfying these conditions, form a 2D BZ of 
double QCB. 
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Figure 3.2: Fermi surface of 2D metallic quantum bars in the absence of charge transfer 
between wires. 

Due to the inter-array interaction, the plasmons in QCB (see Figs. 13.31 IC.ll below) 
exist in a 2D Brillouin zone. These excitations are characterized by the quasimomentum 
Q = {QiiQ.2)- However, in case of weak interaction, the 2D waves are constructed from 
the ID plane waves in accordance with the above procedure. Two sets of ID plane waves 
form an appropriate basis for the description of elementary excitations in QCB in close 
analogy with the nearly free electron approximation in conventional crystalline lattices. 
It is easily foreknown that a weak inter-array interaction does not completely destroy 
the above quasimomentum classification of eigenstates, and the 2D reconstruction of the 
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spectrum may be described in terms of wave mixing similarly to the standard Bragg 
diffraction in a weak periodic potential. Moreover, the classification of eigenstates of an 
empty superlattice may be effectively employed for the classification of energy bands in a 
real QCB superlattice where the super-periodicity is imposed by interaction. 

Complete kinematics of an empty super-chain (wave functions, dispersion laws, rela- 
tions between quasiparticle second quantization operators) is developed in Appendix 1X1 
In terms of these ID Bloch functions (see Eqs. (jA.ljl . (|A.2J) of Appendix El we construct 
the 2D basis of Bloch functions for an empty superlattice 

^ss'q(r) = 4> sqi (Xl)4>siq 2 (x2)- (3.1) 

Here s,s' = 1,2,..., are the band numbers, and the 2D quasimomentum q = (91,92) 
belongs to the first BZ, \qj\ < Q/2. The corresponding eigenfrequencies 

LO ss ,(q) = Uu(qi) + ^'(92)- 

Here ujj S (qj) = ui s (qj), and u} s (qj) are formed by two ID acoustic branches propagating 
along the j-th array and belonging to the band s. Their explicit form is defined in Ap- 
pendix 1X1 (Eq. (|A.2|) ) . Alternatively, each mode with quasimomentum q in the energy band 
s (reduced BZ description) propagating along the j-th array can be described by the wave 
vector q+nij (extended BZ description), where mj = rrijQej andm 3 - = (— l) s sign(9j)[s/2]. 
We will use both these bases in the next Section when constructing the excitation spectrum 
of QCB within the reduced band scheme and in the next Chapters where we formulate 
the spectroscopic principles of QCB within the extended band scheme. 

3.3 Hamiltonian 

When turning to the description of interaction in a QCB, one should refer to a real 
geometry of crossbars, and recollect the important fact that the equilibrium distance 
between two arrays is finite and large enough to suppress direct electron tunnelling Pj. 
We neglect also the elastic and van der-Waals components of the interaction between 
real nanotubes, because these interactions are not involved in the formation of collective 
excitations in QCB. Then the full Hamiltonian of the QCB is 

Hqcb = H 1 + H 2 + #12, (3.2) 

where the Hamiltonian H% (|2.24|) describes the ID boson field in the first array, the 
Hamiltonian H2 is obtained from Hi after permutation 1 <-> 2 in the arguments of the 
fields. 
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The inter-wire interaction includes both interactions between wires from the same 
array (intra-array interaction) and wires from different arrays (inter-array interaction). 
The latter is the Hamiltonian 1)2.28)1 which results from a contact capacitive coupling in 
the crosses of the bars. Physically, the short-range inter-array interaction V (x%— nia, n^a— 
X2) represents the screened Coulomb interaction between charge fluctuations around the 
crossing point (rija, n^a). We assume that the crossed nanotubes are suspended in an 
unpolarized medium, and screening arises due to intra-wire interaction. The nanotube 
diameter is the only physical parameter which determines the screening length r$ in a 
tube (see e.g. Ref.[iBj). We describe the re-distribution of a charge in tube j induced by 
the interaction with tube i by the envelope function (introduced phenomenologically) 

aa c- ao = a-o, c(o)~i. (3.3) 

^0 

This function is of order unity for |£| ~ 1 and vanishes outside this region. Thus the 
on-cross interaction is introduced as 



r 12 



2 V r o ?"o 

W1 4-1 

d 2 

where r i2 = r £iei — r ^ 2 G2- It is seen from these equations that $(d,^2) vanishes for 
|^i,2| > 1 and satisfies the condition $(0,0) ~ 1. The effective coupling strength is 

2e 2 

V = (3.5) 

In terms of boson field operators 9j, the inter-array interaction is written in the form 
similar to (J22BJ) 

Hint = Vo ^ J dxidx 2 ^( Xl ^ ni<2 , n2a ^ — 1 d Xl 6 1 (x 1 ,n 2 a)d X2 2 {n l a,X2). (3.6) 



711,712 



r r 



As for the inter-wire interaction within each array, one can neglect it for a couple of 
reasons. First, the inter- wire distance within the same array is much larger than the inter- 
array distance. Second, this interaction is irrelevant in the long- wave limit [23] . Thus Eq. 
(13. 6 j) is the full interaction Hamiltonian. 

In the quasimomentum representation (|3.1|) the full Hamiltonian (|3.2|) acquires the 
form, 

Hv 2 h 2 

Hqcb = ~2{j XMsq 7 ^ + 2^ X^' ss '^q^' Vt i' ( 3 - 7 ) 

3=1 s,q jj'=l ss'q 

where 8j SCl and 7Tj Sq are the Fourier components of the boson fields 8j and ttj. 
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The matrix elements for inter-array coupling are given by: 

fyj'as'q = u s (qj)uj s ,(qf) [6jf8 ss/ + </>$ ss ' q (1 - 8^)] , <p = 9 . (3.8) 

Here 0J s {qj) = v([s/2]Q + (— l) s_1 |gj|) are eigenfrequencies of the "unperturbed" ID mode 
(see Eq. (IA.2J) of Appendix |XJ), pertaining to an array j, band s and quasimomentum 
q = qjej. The coefficients 



-l) s + s 'sign( gi g 2 ) | ^i^2$(6,6)e- Jro(9l€l+ ^ 2) < gi (r 6)^ (?2 (ro6), (3.9) 



($ss'q = ^s'sq) are proportional to the dimensionless Fourier component of the interaction 
strengths, where the Bloch amplitudes u sqj (xj) are given by Eq. (jA.2|) of Appendix lAl 

The Hamiltonian ()3.7|) describes a system of coupled harmonic oscillators, which can be 
exactly diagonalized with the help of a certain canonical linear transformation (note that it 
is already diagonal with respect to the quasimomentum q). The diagonalization procedure 
is, nevertheless, rather cumbersome due to the mixing of states belonging to different 
bands and arrays. However, it will be shown below that provided d ^> r , a separable 
potential approximation is applicable, that significantly simplifies the calculations. 

3.4 Approximations 

As it has already been mentioned, we consider the rarefied QCB with short range capac- 
itive interaction. In the case of QCB formed by nanotubes, this is a screened Coulomb 
interaction at a distance of the order of the nanotube radius Rq [46 , therefore r ~ R Q . 
The minimal radius of a single-wall carbon nanotube is about Ro = 0.35 -j- OAnm (see 
Ref. |54|). The inter-tube vertical distance d in artificially fabricated nanotube networks 
is estimated as d ~ 2nm (see Ref. [S]). Therefore the ratio r^/d 2 ~ 0.04 is really small 
and the dimensionless interaction $(^1,^2) f)3.4|) in the main approximation is separable 

$(a,6)~$o(6,6) = C(£i)C(£2). (3.io) 

To diagonalize the Hamiltonian ()3.7|) . one should solve the system of equations of 
motion for the field operators. The generalized coordinates 9 satisfy the equations 



[u 2 s {qi) - lo 2 ] #i sq + VsCs(qi)v s (qi)— (ftO^V (<?2)#2 S 'q = 0, s = 1,2, . . . , (3.11) 

a 

s' 

and the similar equations obtained by permutation 1 <-> 2. Here 

C.(?) = (-irsign(g) J ^C(0e ,ro9 ^ 9 (ro0, (3-12) 
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the Bloch amplitudes u sq (r ^) are denned by Eq. (jA.2|) of Appendix^ and 



a 



, =(*P) ■ (3-13) 
r J \ hv J 



Due to separability of the interaction, the equations of motion (j3.11j) can be solved exactly 
The corresponding eigenfrequencies are determined by the characteristic equation 

1 

e" 

where 



FJ^)FJ^) = - (3.14) 



FM)= r^a^MM (3 . 15 ) 

The function F q . {u 2 ) has a set of poles at u 2 = w 2 (qj) (s = 1,2,3,...). For squared 
frequency smaller than all squared initial eigenfrequencies u 2 (qj), i.e., within the interval 
[0, u) 2 (qj)], this is a positive and monotonically increasing function. Its minimal value F 
on the interval is reached at uj 2 = 0, and it does not depend on quasimomentum qj 

F v( Q ) = ~ EC 2 fe) = / deC a (0 = Fo (3.16) 

a s J 

(here Eqs. fl3.15j) and ()3.12j) are used). If the parameter e is smaller than its critical 
value e c = 1/Fq, then all solutions u 2 of the characteristic equation are positive. When 
e increases, the lowest QCB mode softens and its frequency vanishes in the entire BZ at 
e = e c . For exponential charge density distribution = exp(— 1£|), one obtains e c ~ 1. 
In our model, the dimensionless interaction e in Eq. (j3.13j) can be written as 

For nanotube QCB, the first factor within parentheses is about 0.35. The second one 
which is nothing but the corresponding QCB "fine structure" constant, can be estimated 
as 0.9 (we used the values of g = 1/3 and v = 8 x 10 7 cm/sec, see Ref. [IB])- Therefore e 
approximately equals 0.1, so this parameter is really small. Thus the considered system 
is stable, its spectrum is described by Eqs. (|3.14|) . (|3.15J) with a small parameter e. 

The general Eq. ()3.14j) reduces in the infrared limit q, u — > to an equation describing 
the spectrum of two coupled sliding phases, i.e., 1 : 1 arrays in accordance with classifi- 
cation scheme, offered in Ref. [23j . Equation (3.13) of this paper is the long wave limit 
of our equations (jB.2|) and (jB.13|) derived in Appendix El Therefore, the general analysis 
of stability of the LL fixed point is applicable in our approach. 
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3.5 Spectrum 



Due to the weakness of the interaction, the systematics of unperturbed levels and states is 
grossly conserved, at least in the low energy region corresponding to the first few energy 
bands. This means that perturbed eigenstates could be described by the same quantum 
numbers (array number, band number and quasimomentum) as the unperturbed ones. 
Such a description fails in two specific regions of reciprocal space. The first of them is the 
vicinity of the high-symmetry lines qj = nQj/2 with n integer. Indeed, as it follows from 
the equations of motion (j3.11|) . around these lines the inter-band mixing is significant. 
These lines with n = ±1 include BZ boundaries. The second region is the vicinity of the 
lines where the resonance conditions are fulfilled 

u 2 s {q l )=u J 2 (q 2 ). (3.18) 

Here inter-array mixing within the same energy band s is significant. In what follows we 
will pay attention first of all to these two regions because in the rest of the BZ the initial 
systematics of the energy spectrum can be successfully used. 

Equations (j3.11|) . ()3.14|) . describing the wave functions and the dispersion laws are 
analyzed in Appendix [5] We describe below some of these dispersion curves for a square 
QCB based on this analysis (the case of tilted QCB is described in Appendix IU|) . 

We start with the simplest case of square QCB formed by identical wires. This means 
that all parameters (wire length, space period, Fermi velocity, LL parameter, screening 
radius) are the same for both arrays. The corresponding BZ is also a square (see Fig. 
13.3)1 . The resonant lines are the diagonals of BZ. 
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Figure 3.3: Two dimensional BZ of square QCB. 

In the major part of the BZ, for quasimomenta q lying far from the diagonals, each 
eigenstate mostly conserves its initial systematics, i.e. belongs to a given array, and mostly 
depends on a given quasimomentum component. Corresponding dispersion laws remain 
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linear being slightly modified near the BZ boundaries only. The main change is therefore 
the renormalization of the plasmon velocity. 



Figure 3.4: The energy spectrum of QCB (solid lines) and noninteracting arrays (dashed 
lines) for quasimomenta at the lines AT, VXi, X\W ', and WT of BZ. 

i. General case: The point q lies away from both the high symmetry lines 
and the resonant lines. This case is illustrated in the left part of Fig. 13.41 where we 
display dispersion curves corresponding to quasimomenta belonging to a generic TA line 
in the BZ. In what follows we use (j, s) notations for the unperturbed boson propagating 
along the j-th array in the s-th band. Then the lowest two curves in this part of Fig. 13.41 
are, in fact, the slightly renormalized dispersion of (2, 1) and (1, 1) bosons. The middle 
curves describe (1,2) and (2,2) bosons, and the upper curves are the dispersion of (2,3) 
and (1,3) bosons. It is seen that the dispersion remains linear along the whole line TA 
except at the nearest vicinity of the BZ boundary (point A in Fig. 13.41 and point T for 
highest energy bands). 

ii. Inter-band resonance in one of the arrays: the point q lies on a high 
symmetry line of only one array. This case is illustrated by the lines TXi and X{W 
in Fig. 13.31 Dispersion curves corresponding to quasi momenta lying at the line q 2 = 0, 
< qi < Q/2 (line TX 1 in Fig. and the BZ boundary q 1 = Q/2, < q 2 < Q/2 
(line XiW in Fig. I3.3|) are displayed in the central parts in Fig. I3.4J) . The characteristic 
feature of these lines is the intra-band degeneracy in one of two arrays. Indeed, in zero 
approximation, two modes (2, s), s = 2,3, propagating along the second array with a 
quasimomentum lying in the line FXi are degenerate with unperturbed frequency u = 1. 
The interaction lifts the degeneracy. This interaction appears to be repulsive. As a 
result, the lowest of the two middle curves in Fig. 13.41 corresponds to (2, it) boson, and 
the upper among them describes (2,g) boson. Here the indices g,u denote a boson 
parity with respect to the transposition of the band numbers. Note that the (2, g) boson 
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exactly conserves its unperturbed frequency 00 = 1. The latter fact is related to the square 
symmetry of the QCB. 

Other curves correspond to almost non perturbed bosons of the first array. The lowest 
two curves describe the dispersion of the (1,1) and (1,2) waves. Plasmons in the third 
band (1, 3) are described by the uppermost curve in the figure. Their dispersion laws are 
nearly linear, and deviations from linearity are observed only near the boundary of the 
BZ (point X% in Fig. 13.3)1 and near the point V for the highest bands. 

Similarly, in zero approximation, two modes (l,s), s = 1,2 (3,4), propagating along 
the first array with a quasimomentum lying in the line XiW are degenerate with un- 
perturbed frequency uo = 0.5 (1.5). The interaction lifts the degeneracy. As a result, 
the lowest of two middle (highest) curves in Fig. 13.41 corresponds to (2, u) boson, and 
the upper of them describes (2, g) boson. As in the previous case, (2, g) boson exactly 
conserves its unperturbed frequency to = 0.5 (1.5). 

Other curves correspond to almost non perturbed bosons of the second array. The 
lowest curve describes the dispersion of the (2, 1) wave. The two middle curves describe 
the dispersion of (2, 2) and (2, 3) plasmons. Their dispersion laws are nearly linear, and 
deviations from linearity are observed only near the corner of the BZ (point W in Fig. 
13. 3|) as well as in the vicinity of the point X\ of the BZ. 

iii Inter-band resonance in both arrays: the point X\ (X2) is a crossing 
point of the two high symmetry lines away from all resonant lines. This case is 
illustrated by the points X± and X 2 in Fig. 13.31 Consider for example point Xl. Here q\ = 
Q/2, Q 2 = 0. In zero approximation, two modes (l,s), s — 1,2 (s = 3,4), propagating 
along the first array are degenerate with unperturbed frequency uo = 0.5 (u = 1.5). The 
lower (higher) two lines correspond to even (1, g) and odd (l,u) superpositions of the 1-st 
array states of the first and second (third and fourth) bands. Similarly, two modes (2, s) 
with s = 2, 3 are degenerate in zero approximation with unperturbed frequency u = 1. 
Therefore the middle two lines describe the same superpositions of the 2-d array states 
from the second and third bands. 

iv. Inter-array resonance: The point q lies only on one of the resonant lines 
away from the high symmetry lines. This case is illustrated by the diagonal TW 
of BZ. Consider now dispersion relations of modes with quasi-momenta on the diagonal 
TW of BZ and start with q not too close to the BZ corner W (q± = q 2 = Q/2) and the 
r point. This diagonal is actually one of the resonance lines. Two modes in the first 
band corresponding to different arrays are strongly mixed. They mostly have a definite 
j-parity with respect to transposition of array numbers j = 1,2. Interaction between 
these modes appears to be attractive (repulsive) for q\q 2 > (qiq 2 < 0). Therefore the 
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odd modes (u,s), at the BZ diagonal TW s = 1,2, correspond to lower frequencies and 
the even modes (g, s) correspond to higher ones. The corresponding dispersion curves are 
displayed in the right part of Fig. 1.141 

v. Inter-array and inter-band resonance: The point q lies at the cross of 
two resonant lines. There are two points, the BZ corner W and the BZ center F in 
Fig. 13.31 At the BZ corner qi = q 2 = Q/2 (point W in Fig. 13. 3|) all four initial modes 
j = 1,2 and s = 1, 2 (s = 3, 4) are degenerate in the lowest approximation. This four-fold 
degeneracy results from the square symmetry of BZ (the resonant lines are diagonals of the 
BZ). Weak inter-array interaction partially lifts the degeneracy. However, the split modes 
have a definite s-parity with respect to transposition of band numbers s = 1, 2 (s = 3, 4). 
The highest frequency corresponds mostly to (g,u) boson, symmetric with respect to 
transposition of array numbers, but antisymmetric with respect to the transposition of 
band numbers. The lower curve describes a (u, u) boson with both j-parity and s-parity 
odd. The two middle modes with even band parity, (g, g) and (u,g) bosons, remain 
degenerate and their frequencies conserve the unperturbed value u = 0.5 (u = 1.5). This 
also results from the square symmetry of QCB f)3.4j) . 

Similar behavior is observed in the BZ center T. All four initial modes j = 1,2 and s = 
2, 3 are degenerate in the lowest approximation. Weak inter-array interaction partially lifts 
the degeneracy. The highest frequency corresponds mostly to (g, u) boson,the lower curve 
describes a (u, u) boson, and the two middle bosons, (g,g) and (u,g), remain degenerate 
and their frequencies conserve the unperturbed value io = 1. 
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Figure 3.5: Lines of equal frequency of the lowest mode for QCB (solid lines) and for 
noninteracting arrays (dashed lines). The lines 1,2,3 correspond to the frequencies u)\ = 
0.1, u 2 = 0.25, u 3 = 0.4. 

All these results show that the quantum states of the 2D QCB conserve the quasi 
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ID character of the Luttinger-like liquid in major part of momentum space, and that 
2D effects can be successfully calculated within the framework of perturbation theory. 
However, bosons with quasimomenta close to the resonant line (diagonal OC) of the BZ 
are strongly mixed bare ID bosons. These excitations are essentially two-dimensional, 
and therefore the lines of equal energy in this part of the BZ are modified by the 2D 
interaction (see Fig. 13. 5|) . It is clearly seen that deviations from linearity occur only in a 
small part of the BZ. The crossover from LL to FL behavior around isolated points of the 
BZ due to a single-particle hybridization (tunnelling) for Fermi excitations was noticed in 
Refs. |2Z1 , where a mesh of horizontal and vertical stripes in superconducting cuprates 
was studied. 



3.6 Correlations and Observables 

The structure of the energy spectrum analyzed above predetermines optical and transport 
properties of QCB. We consider here three types of correlation functions manifesting 
dimensional crossover in QCB. 

3.6.1 Optical Absorption 

We start with ac conductivity 

(Jjy(q, u) = o-'y (q, u) + i<rL(q, u). 



The real part cL/(q, oj) determines an optical absorption. The spectral properties of ac 
conductivity are given by a current-current correlator 



a Jj/ ( q ,u) = - / rfte^([j jlq (t),4 lq (0) 



(3.19) 



Here J Jsq = y/2vgiTj S q is a current operator for the j-th array (we restrict ourselves to the 
first band, for the sake of simplicity). 

The current-current correlator for non-interacting wires is reduced to the conventional 
LL expression [I], 



•W*),«4lq(°) 



-2ivguiq sin(u; lq £)<5 



j j 



with metallic-like peak 



cr' -,(q,o; > 0) = nvg5(uj - Wi q )5,y. 



(3.20) 
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For QCB this correlator is calculated in Appendix EI Its analysis leads to the following 
results. The longitudinal absorption 

<r' n (q,u) oc (1 - 0i q )5(w - w lq ) + 0? q 5(w - u 2q ) 

contains well pronounced peak on the modified first array frequency and weak peak at 
the second array frequency (the parameter 0i q = v^Ci(9i)CifeVi(<?i Vife) is small). 
The modified frequencies u>i q and u) 2q coincide with the eigenfrequencies cj +lq and c<j_ 2q 
respectively, if cj lq > cj 2q . In the opposite case the signs +, — should be changed to the 
opposite ones. 

The transverse absorption component contains two weak peaks 

(q, k>) OC 01q [S(U - Ulq) + 6(u - £ 2 q)] • 

At the resonant line, the results are drastically modified. Both longitudinal and trans- 
verse components of the optical absorption contain two well pronounced peaks corre- 
sponding to slightly split modified frequencies 

a' u (q,uj) oc - [<J(o; -o)i q ) + 5{u - a) 2q )] . 
3.6.2 Space Perturbation 

One of the main effects specific for QCB is the appearance of non-zero transverse momentum- 
momentum correlation function. In space-time coordinates (x, t) its representation reads, 

Gi 2 (x;t) = i ([vr 1 (xi,0;t),vr 2 (0,x 2 ;0)]) , x = (x 1 ,x 2 ). 

This function describes the momentum response at the point (0,x 2 ) of the second array 




Figure 3.6: The transverse correlation function Gi 2 (xi, x 2 ; t) for tq = 1 and vt = 10. 
for time t caused by an initial (t = 0) perturbation localized in coordinate space at the 
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point (xi,0) of the first array. Standard calculations similar to those described above, 
lead to the following expression, 



oo 

,2 



n ( A V ° r o f x a u (i M i ■ (i \ ■ (i ,vk 2 sm(vk 2 t) -vhsm^krf) 

Gi 2 (x;t) = — — / dk 1 dk2<Pi{ki)^2{h)kiK2sm(k 1 x 1 )sin(k2X2) 
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where (f>j(k) is the form-factor (j3.12j) written in the extended BZ. This correlator is shown 
in Fig. 13.61 It is mostly localized at the line determined by the obvious kinematic condition 
+ |^2| — vt ("horizon of events"). The time t in the r.h.s. is the total time of plasmon 
propagation from the starting point (xi,0) to the final point (0,^2) or vice versa, along 
any of the shortest ways compatible with a restricted geometry of the 2D grid. The 
finiteness of the interaction radius slightly spreads this peak and modifies its profile. 

3.6.3 Rabi Oscillations 

Further manifestation of the 2D character of QCB system is related to the possibility 
of periodic energy transfer between the two arrays. Consider an initial perturbation 
which excites a plane wave with amplitude #0 within the first array in the system of 
non-interacting arrays, 

Oi(xi,n 2 a;t) = 8 sin(giXi + q 2 n 2 a - v \qi\t). 

If the wave vector q, satisfying the condition |q| << Q/2, is not close to the resonant 
line of the BZ, weak inter-array interaction <p (|3.8jl slightly changes the 9\ component and 
leads to the appearance of a small 82 ~ component. But for q lying on the resonant 
line (v\q±\ = v\q2\ = w q ), both components within the main approximation have the same 
order of magnitude, 

9i(xi,n 2 a; t) = 9 cos ( ^0 iq u; q t J sin(giXi + q 2 n 2 a - Vqt), 



6 2 (n 1 a, x 2 ; t) 



9 sin Q</»i q o; q tJ cos(ginia + g 2 ^2 - w q t). 



This corresponds to 2D propagation of a plane wave with wave vector q, modulated by a 
"slow" frequency ~ <puj. As a result, beating arises due to periodic energy transfer from 
one array to another during a long period T ~ [cf)uo)~ l (see Fig. 13. 7J) . These peculiar 
"Rabi oscillations" may be considered as one of the fingerprints of the physics exposed in 
QCB systems. Similar periodic energy transfer between three arrays can be observed in 
triple QCB (see Section ID. 31 of Appendix |D|) . 
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Figure 3.7: Periodic energy exchange between arrays (Rabi oscillations). 

3.7 Conclusions 

We discussed in this Chapter the kinematics and dynamics of plasmon spectrum in QCB. 
These nanostructures may be fabricated from single- wall carbon nanotubes (HI EJ- On 
the one hand, QCB is promised to become an important component of future molecular 
electronics [13 OS- On the other hand, the spectrum of elementary excitations (plasmons) 
in these grids possesses the features of both ID and 2D electron liquid. As is shown in 
Refs. j2H EBj and confirmed in the present study, the energy spectrum of QCB preserves 
the characteristic properties of LL at |q|, uj — > 0,. At finite q, u the density and momentum 
waves in QCB may have either ID or 2D character depending on the direction of the 
wave vector. Due to inter-array interaction, unperturbed states, propagating along arrays 
are always mixed, and transverse components of correlation functions do not vanish. 
For quasi- momentum lying on the resonant lines of the BZ, such mixing is strong and 
transverse correlators have the same order of magnitude as the longitudinal ones. Periodic 
energy transfer between arrays ( "Rabi oscillations" ) is predicted. 

The crossover from ID to 2D regime may be experimentally observed. One of the 
experimental manifestations, i.e. the crossover from isotropic to anisotropic (spatially 
nonuniform) conductivity was pointed out in Ref. j2H]. The current may be inserted 
in QCB at a point on an array j and extracted from another array i at a distance r. 
Then a temperature dependent length scale l(T) arises, so that for r ^> I the resistance is 
dominated by small q and therefore, the current is isotropic. In the opposite limit r < / 
the dependence of the current on the points of injection/extraction may be detected. At 
T = the length I becomes infinite, and current can only be carried along the wires. 
These effects are in fact manifestations of the LL behavior of the QCB in the infrared 
limit. 
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To observe the crossover at finite {co>,q}, one should find a way of exciting the cor- 
responding plasmon modes. Then, scanning the u; ((71,(72) surfaces, one may in principle 
detect the crossover from quasi ID to 2D behavior in accordance with the properties 
of the energy spectra presented in Section 13.51 and Appendices O and |DJ Plasmons in 
QCB may be excited either by means of microwave resonators or by means of interaction 
with surface plasmons. In the latter case one should prepare the grid on a corresponding 
semiconductor substrate and measure, e.g., the plasmon loss spectra. The theory of these 
plasmon losses will be presented in Chapter 
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Chapter 4 



Infrared Spectroscopy of Quantum 
Crossbars 

4.1 Introduction 

In this Chapter we consider various possibilities of direct observation of plasmon spectra 
at high frequencies and wave vectors by the methods of infrared (IR) spectroscopy. The 
QCB plasmons can be treated as a set of dipoles distributed within QCB constituents. 
In a single wire, the density of the dipole momenta is proportional to the LL boson field 
6{x) (x is the coordinate along the wire). Few sets of coupled ID dipoles form a unique 
system which possesses the properties of ID and 2D liquid depending on the type of 
experimental probe. Some possibilities of observation of ID — > 2D crossover in transport 
measurements were discussed in Ref. [To] . 

In transport measurements, the geometric factors regulate the crossover from anisotropic 
to isotropic resistivity of QCB: one may study the dc response for a field applied either 
parallel to one of the constituent arrays or in arbitrary direction. One may also study 
spatially nonuniform response by means of two probes inserted at different points of QCB 
and regulate the length scale, i.e., the distance between the two probes in comparison 
with the periods of the crossbar superlattice. These methods give information about the 
nearest vicinity of LL fixed point at (q, (J,T) — > 0. 

Several crossover effects such as appearance of non-zero transverse space correlators 
and periodic energy transfer between arrays ("Rabi oscillations") were discussed in the 
previous Chapter. Unlike transport measurements, the methods of infrared spectroscopy 
provide an effective tool for investigating the excitation spectrum in a rather wide (q, u) 
area well beyond the sliding phase region 61J. We will show that the IR spectroscopy 
allows scanning of the 2D Brillouin zone in various directions and thereby elucidates 
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dimensional crossover in the high symmetry points of the BZ. 

The direct manifestation of dimensional crossover is through the response to an exter- 
nal ac electromagnetic field [60, 61, 62J. To estimate this response one should note that the 
two main parameters characterizing the plasmon spectrum in QCB are the Fermi velocity 
v of electrons in a wire and the QCB period a (we assume both periods to be equal). These 
parameters define both the typical QCB plasmon wave numbers q = |q| ~ Q = 2n /a and 
the typical plasmon frequencies u ~ ojq = vQ. Choosing according to Refs. (HI HE] 
v ~ 0.8 • 10 6 m/sec and a ~ 20 nm, one finds that characteristic plasmon frequencies lie 
in the far infrared region uo ~ 10 14 sec" 1 , while characteristic wave vectors are estimated 
as q ~ 10 6 cm" 1 . 

Here we study high frequency properties of the simplest double square QCB (general- 
ization to more complicated geometries is straightforward). We start from QCB interact- 
ing with an external infrared radiation. The plasmon velocity v is much smaller than the 
light velocity c and the light wave vector k is three orders of magnitude smaller than the 
characteristic plasmon wave vector Q corresponding to the same frequency. Therefore, 
infrared radiation incident directly on an isolated array, can not excite plasmons at all 
(it could excite plasmon with u ^ 0). However in QCB geometry, each array serves as a 
diffraction lattice for its partner, giving rise to Umklapp processes of wave vectors nQ, n 
integer. As a result, excitation of plasmons in the BZ center q = with frequencies nvQ 
occurs. 

To excite QCB plasmons with q ^ one may use an additional diffraction lattice (DL) 
with period A > a coplanar to the QCB. Here the diffraction field contains space har- 
monics with wave vectors 2ir Mj A, M integer, that enables one to scan plasmon spectrum 
within the BZ. Dimensional crossover manifests itself in the appearance of additional 
absorption lines when the wave vector of the diffraction field is oriented along specific 
directions. In the general case one observes the single absorption lines forming two sets 
of equidistant series. Instead of that, in the main resonance direction (QCB diagonal) 
an equidistant series of split doublets can be observed. In the case of higher resonance 
direction, absorption lines form an alternating series of singlets and split doublets demon- 
strating new type of dimensional crossover related to the frequency change with direction 
fixed. 

The structure of the present Chapter is as follows: In Sections 14.21 and 14.31 we study 
interaction of QCB with an external field. In Section 14.21 we consider the case when 
the incident infrared radiation falls directly on the QCB. In Section T4.3I we study possi- 
ble scanning of QCB spectrum with the help of an external DL. In the Conclusions we 
summarize the results obtained. 
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4.2 Long Wave Absorption 



In the case of a dielectric substrate transparent in the infrared region, one can treat QCB 
as an isolated grid (without substrate) interacting directly with the incident radiation. 
Consider the simplest geometry where an external wave falls normally onto QCB plane, 
and its electrical field 

E = .Eb^i cos (kr — ut) 

is parallel to the lower (first) array (see Fig. I4.1l for details). In this geometry the field E is 
longitudinal for array 1 and transverse for the array 2. The eigenfrequencies of transverse 
modes in array 2 substantially exceed the IR frequency of the incident wave and even 
the standard LL ultraviolet cutoff frequency. Thus, the incident wave can be treated as 
a static polarization field for this array, and the factor cos tut can be omitted. Then, 
the polarization waves in array 2 form a longitudinal diffraction field for array 1 with 
quasi wave vectors nQ (n integer). Further, the characteristic order of magnitude Q of a 
QCB plasmon wave vector is much larger than the wave vector k of the incident light, so 
one can assume the latter to be equal to zero from the very beginning. Then, the light 
wavelength is much larger than a nanotube diameter and the geometrical shadow effect 
can be neglected. As a result the total field which affects array 1 consists of an external 
field and a diffraction field produced by a static charge induced in array 2. 




Figure 4.1: The incident field orientation with respect to QCB. The axes X\ and %i are 
directed along the corresponding arrays, and d is the inter-array vertical distance (along 
the X3 axis). 

To calculate the diffraction field, we consider first the field E° produced by the quan- 
tum wire of array 2 which is located at x% = X3 = and labelled by n\ = 0. The large 
distance between the wire under consideration and its neighbor partners from the same 
array allows us to neglect the influence of the charges induced on them. The static po- 
tential on the surface of the wire includes external potential of an incident field and the 
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potential $° of the charge induced on the wire. On the other hand, this static potential 
should be equal to a constant which we choose to be zero. In cylindrical coordinates 
r, X2, Xi — r cos 23 = r sin this condition reads 

$ o (R o ,ti,x 2 ) = E o r o coB0. (4.1) 

Outside the wire, the induced potential $° satisfies the Laplace equation A$° = 0. Solving 
this equation with boundary condition (|4.1j) we obtain the static part of the induced 
potential 

$°(r,tf,x 2 ) = ^costf 
r 

and the corresponding static part of the induced field along the x\ direction, 

,y> 2 ( ,-yj2 ,-y»2 \ 

'0 l x 3 x l) 



E%(x u x 3 ) = -Eq 



{xl + x\f 



The first component of the diffraction field is the sum of the fields induced by all wires of 
the upper array, 

Etixr, t) = cosut ~ n ia , -d) = -Eq cosujt 7^ 7 ^ ~ • (4-2) 

This field is a periodic function of x\ with a period a. Therefore, its Fourier expansion 
contains only wave vectors k\ n = nQ (n is the order of diffraction). This means that only 
frequencies uj n = nvQ can be excited. In this case it is more convenient to expand the 
field over Bloch eigenf unctions of an "empty" wire jSH]- These functions are labelled by 
quasimomentum q±, \qi\ < Q/2, and the band number s. The expansion includes only 
qi = components and has the form 

E^Xi't) = CQ8UJty)E[ a /2]U a (xi), 
s 

where u s (x) is the q\ = + Bloch amplitude u sqi (x) (|A.2j) within the s-th band and 

2 

E n = -E ^j-nQde- nQd . (4.3) 
ad 

The excited eigenfrequency u n = 0J[ s /2] belongs simultaneously to the top of the lower 
even band with number s = 2n and to the bottom of the upper odd band with number 
s = 2n + 1 (this is the result of E(x) = E(— x) parity). The incident field cannot excite 
plasmons at all and we do not take it into account. 

Turning to the (q, s) representation with the help of the expansion 



L 
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and similarly for 6 2 and 7r li2 , one easily sees that only the q = components are involved 
in the interaction of plasmons with the incident radiation. 

Consider an initial frequency to close to w„. In a resonant approximation, only four 
equations of motion for the "coordinate" operators 9 S with s = 2n, 2n + 1 are relevant 

01,2n + ^l,2n + (#2,2n ~ 02,2n+l) = L fn COS ut, 

#1,271+1 + Un@l,2n+1 ~ (6*2,271 ~ #2,2n+l) = Lf n COSUjt, 

o 2 ,2n + uj 2 n e 2M + (f)uj 2 n (e^ 2n -e^ 2n+1 ) = o, (4.5) 

#2,2n+l + ^n^2,2n+l — <f>^n (#l,2n — $l,2n+l) = 0, 



where 



V2vge 

Jn t /— LJ 'n 



and we assume that Cs(li,2) (j3.12j) is equal to 1 for the first few bands. 

The homogeneous part of this system defines four eigenfrequencies: u) gs = u u ^ g = ui n , 
^g/u^u ~ kVi(l =t 4>)- The corresponding eigenvectors are symmetrized combinations of the 
four operators which enter Eq. (j4.5|) . They have a fixed parity with respect to permutation 
of arrays (the first index) and neighboring bands (the second index). Only two modes 
(even with respect to band index) 

@g/u,g — r (^l,2n + $1,271+1 =t $2,2n i #2,2n+l) > 

interact with an external field. Therefore only the unperturbed frequency u n = u gg = u ug 
will be absorbed. The two equations of motion for the operators 9 gg ^ ug have the same 
form 

6 a + 2-y8 a + tu 2 8 a = Lf n cos tot, 

where a = gg,ug. Employing standard procedure in the vicinity of the resonance \u — 
uj n \ <C 0J n immediately yields the relative absorption of the Lorentz type 

Io g hc\ad) {u-uj n f + ^ [UQde J ' ^ 



where 



cL 2 2 
-El 



is the intensity of light that falls on the QCB. 

Due to the exponential term in the r.h.s of Eq. (|4.3|) . E n decreases fast with n and only 
the first few terms contribute to absorption. The characteristic dimensionless scale of the 
induced field ^/(ad) for typical values of QCB parameters equals 0.004. We tabulate 
below the lowest dimensionless Fourier components of the induced field. 
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The results show that one can hope to probe at least the first five spectral lines corre- 
sponding to oj n with n — 1, 2, . . . , 5. 

The width of the absorption line ()4.6)1 is governed by an attenuation coefficient 7. We 
introduce it phenomenologically, but one can (at least qualitatively) estimate its value. 
The attenuation is caused by decay of plasmon into phonons. The one phonon decay 
of the plasmon with wave number k and frequency u = v\k\ into a single phonon with 
the same u and k occurs at a single point in ID and does not yield finite attenuation 
at all. Multi-phonon decay is weak because of the small anharmonic coupling within the 
wire. As a result, the form of the absorption lines should be determined mainly by the 
instrumental line-width. 



4.3 Scanning of the QCB Spectrum within the BZ 

Within a geometry considered in the previous subsection, one can probe plasmon spectrum 
only at the BZ center. To study plasmons with nonzero wave vectors one should add to 
the system an external diffraction lattice (DL), namely, a periodic array of metallic stripes 
parallel to the Y axis (see Fig. 14. 2)1 . The DL plane Z = is parallel to the QCB planes 
Z = —D for the upper second array and Z = — (D + d) for the lower first array (the Z 
axis is parallel to the 23 axis). The distance D between DL and second array is of the 
same order as the inter-array distance d = 2 nm. The angle between the DL wires and 
the second array is if (0 < ip < tt/2). To get a wave number K of a diffraction field much 
smaller than Q one needs a DL with a period A much larger than the QCB period a. In 
the following numerical estimations we choose A pa 200 nm. 























































































































jj 
















a 


a 




\ 




— > 



















































X 



x 



Figure 4.2: QCB and DL. The X, Y axes are oriented along the DL stripes and the wave 
vector K of the diffraction field respectively. The DL (QCB) period is A (a). 
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Consider an incident field with electric vector E = E e x cos(kr — uot) oriented along 
the X axis (perpendicular to the DL wires). The radius Rq of a DL wire is assumed to 
be not much larger than the nanotube radius ro- In this case, light scattering on the DL 
is similar to that considered in subsection 14.21 Then the diffraction field is concentrated 
along the X direction and has the form (compare with Eq. (|4.2j) ) 

?2,"72 IV AT A^2\ 

N (Z* + (X-NA)*f 
The Fourier transform of the diffraction field is 



E X (X, Z, t) = -E cos^t £ R }lf ~l X 



E x {K,Z) = -E^\KZ\e-\ KZ \ 

where K(M) = Ke x — K (sin cp, cos cp) = (K^ K 2 ), with K — 2ttM/A, M being a positive 
integer. This means that all the points K lie on the same ray oriented along the positive 
direction of the X axis. The vector K(M) for a fixed M can be uniquely represented as 
a sum of quasimomenta lying in the first BZ and two reciprocal lattice vectors 

K(M) = q(M) + mi(M) + m 2 (M). 

The field components 

E 1K = E x (K,D + d)sincp, E 2 k = E X (K, D) cosy? (4.7) 

parallel to the quantum wires, can excite plasmons and contribute to the absorption 
process. 

The Hamiltonian describing the interaction of QCB with an external field reads, 
HL 



He = — [/i.k (Vk + <k) + / 2 ,k (Ovc + 9{ K ) ] , (4.8) 
v 9 K 

where 



fj,K = j^^ Ejjt, m = mi + m 2 . (4.9) 

In this Section we are interested not in the form of the absorption line but only in the 
resonant frequencies. Therefore, we do not introduce any phenomenological attenuation. 
The equations of motion for boson fields have the form 

#l,q+mi + ^g 1+mi Q01,q+mi + <P ^ ^q+m^2,q+m 2 = L ^ /l,K^K,q+m, 

ni2 M,m 2 

#2,q+m 2 + ^g 2 + m2 Q#2,q+m 2 + 4> / ] ^q+m^l,q+mi = L ^ /2,K^K,q+m, (4-10) 

mi M,mi 
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where is the dimensionless coupling strength (|3.8|) . $k is given by Eq. ()3.9|) . 

Only the first few terms in the sum over in the r.h.s. of Eq. (j4.10|) really excite the 
QCB plasmons. Indeed, the diffraction field f)4.4j) is proportional to the same dimension- 
less function of the type te~ f (t = \KZj\) as in the previous subsection (see Eq. ()4.3j) ). This 
function has its maximum at t — 1 and differs significantly from zero for 0.2 < t < 2.7. 
For a = 20 nm, D = 2 nm, it is of order unity within the interval 0.18Q < \K\ < 2.13Q 
for the first array {Z\ = D + d), and within the interval 0.36Q < \K\ < 4.26Q for the 
second array (Z\ = D). This means that one can excite the modes of the four lower bands 
(K < 2Q) of the first array and the modes of eighth lower bands (K < 4Q) of the second 
array. 

According to Eqs. (j4.7j) the field -EjK(M) is coupled with plasmons of wave vectors 
q + rrij = q(M) + m.j(M) within the j-th array. The nature of the excited plasmons 
as well as their frequencies depend on the direction of the vector K(M). For simplicity 
we restrict ourselves by acute angles < <p < n/2 describing orientation of both the DL 
and the vector K(M). There are four kinds of dimensional crossover events depending 
on the specific directions in the BZ. Each type of crossover is characterized by its own 
set of absorption lines. The first one takes place in a common case when K(M) for 
any M never reaches neither a resonant direction nor the BZ boundary. The second 
case corresponds to the bisectorial direction (p = tc/4 where the main resonant condition 
ou(Ki) = u{K-i) is fulfilled. The third set of directions is determined by another resonant 
condition io(K\) = uinQ^K^). Finally, the fourth set is formed by directions intersecting 
with the BZ boundaries for some values of M. In what follows we consider these four cases 
separately. 

1. In the general case, the points K(M) for all M are far from the BZ diagonals and 
boundaries. Therefore each of them corresponds to a couple of plasmons mostly propa- 
gating along the j-th array, j = 1,2, with unperturbed frequencies ojk {m) — vKj(M). 
The inter-array interaction slightly renormalizes the eigenfrequencies 



Thus, increasing the frequency of an incident light one observes a set of single absorption 
lines that consists of two almost equidistant subsets with frequencies corresponding to 
excitation of plasmons in the first or second arrays. The frequency spacing between 
adjacent lines within each subset are 




Au% = vAKi = 27rvsmip/A, 



Acj 2 = vAK 2 = 2irvcosif/A, 
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and their ratio depends only on the DL orientation ip 

— — = tan. (p. 
Au 2 

2. In the resonant case ip = 7r/4, the relation K\(M) = K 2 (M) is satisfied for all 
M. Therefore modes propagating along the two arrays are always degenerate. Inter- 
array interaction lifts the degeneracy. Indeed, in the resonant approximation, the coupled 
equations of motion for the field operators read 

9 1K + UJk^ik + 0^^2K = flK, @2K + ^K^2K + (p^K^lK = flK- 

After symmetrization 9 gtU = (9ik ± #2k)/v^2, they have the same form 

q + U 2 a 6 a = f a , 

where u g / u = ioks/^- ± <f> are the renormalized frequencies, f g / u = (/ik ± iWVv^, and 
a = g,u. The amplitudes f 9yU are of the same order of magnitude because the distances 
D and d are different but have the same order of magnitude. As a result, increasing the 
frequency of an incident light one observes an equidistant set of absorption doublets with 
distance ir^/2v/A between adjacent doublets. 

3. Consider now the directions ip determined by the equation 

7r\ nA 



sin {ip ± 



where n and Mq are mutually prime integers. For this direction, two components of the 
first Mo — 1 points K(M) do not satisfy any resonant condition while the Mo-th one does 

K 1 (M )±K 2 (M )=nQ. (4.11) 

With increasing M this situation is reproduced periodically so that all the points K.(pM ) 
with p integer satisfy a similar condition with pn standing instead of n, while all interme- 
diate points are out of resonance. 

In zero order approximation with respect to the inter-array interaction we expect to 
observe two set of absorption lines with frequencies pujj = vKj(pM ), j = 1,2, corre- 
sponding to excitation of plasmons within the prrij(Mo)-th band of the j-th array. The 
ratio of the frequencies u>j is defined by the DL orientation 

— = tan ip. 

However, due to the resonance condition (|4.11j) . a plasmon in the first array with wave 
vector Ki(pMo) and frequency uj\ = vK\(pMo) is coupled with a plasmon in the second 
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array with the same frequency and wave vector K' 2 = ^(npQ — Ki(pM )) (inter- array 
degeneracy). Similarly, a plasmon in the second array with wave vector K 2 (pM ) and 
frequency 002 = vK^ijpMo) is coupled with a plasmon in the first array with the same 
frequency and wave vector K[ = npQ =p K2(pMo). This degeneracy of two modes corre- 
sponding to the same band but to different arrays is lifted by the inter-array interaction. 
As a result one has two sets of doublets instead of two sets of single lines. 

Thus, for such orientation of the DL, increasing the frequency of an incident wave one 
should observe two equidistant sets of single absorption lines with two sets of equidistant 
doublets built in these series 



In the case n = 1 the lower doublet lies in the first energy band, whereas the upper one 
lies in the second band. For A/ a = 10 (that corresponds to the realistic values of the 
parameters a = 20 nm and A = 200 nm) the lowest doublet (p = 1) will be observed for 
example for integer M = 8, at the angle <p(8) ~ 17°, around frequencies co>i(8) = 0.76vQ, 
cu 2 (8) = 0.2AvQ. 

4. A similar behavior will be manifested in the case when the points K p jv/ he at one 
of the BZ boundaries, i.e., they satisfy the relation 



with some specific values j, Mj and n. Such situation is realized at specific angles that 
depend on the integers j,n,Mj. In the vicinity of the points K.(pMj) two frequencies 
corresponding to the unperturbed modes of the j-th array from the np-th and (np + l)-th 
bands coincide. This is the case of inter-band degeneracy that is also lifted by inter-array 
interaction. Due to the square symmetry (invariance with respect to Xj — > —Xj inversion), 
only one of the two components with frequency u = v\Kj(pMj)\ may be excited by a 
diffraction field. Therefore, this case is not distinct from case 1 considered above and two 
sets of equidistant single lines can be observed. 

We emphasize that by studying absorption of light by QCB one can expose, beyond 
the studied above |HH] dimensional crossover with respect to an angle (direction), also the 
occurrence of a new type of crossover with an external frequency as a control parameter. 
This occurs for special directions of type 3 where, with increasing frequency, the set of 
single lines is periodically intermitted by doublets. 
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4.4 Conclusions 



In conclusion, we investigated the possibility of spectroscopic studies of the excitation 
spectrum of quantum crossbars, which possesses unique property of dimensional crossover 
both in spatial coordinates and in (q, a;) coordinates. It follows from our studies that 
the plasmon excitations in QCB may be involved in resonance diffraction of incident 
electromagnetic waves and in optical absorption in the IR part of spectrum. 

In the case of direct interaction of external electric field with QCB, infrared absorption 
strongly depends on the direction of the wave vector q. One can observe dimensional 
crossover from ID — > 2D behavior of QCB by scanning an incident angle. The crossover 
manifests itself in the appearance of a set of absorption doublets instead of the set of 
single lines. At special directions, one can observe new type of crossover where doublets 
replace the single lines with changing frequency at a fixed q direction. 

Dimensional crossover in QCB plays a significant role in all the above phenomena. 
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Chapter 5 



Landau Damping in a 2D Electron 
Gas with Imposed Quantum Grid 

5.1 Introduction 

In this Chapter, QCB interaction with semiconductor substrate is studied. Any surface 
wave excited in the substrate is coupled with QCB-plasmon modes due to the substrate- 
QCB interaction. This interaction might be strong enough compared with the frequency 
spacing of surface waves and QCB plasmons because surface plasmon waves exist in the 
same frequency and wave vector plasmons in QCB (see subsection l5.2.2l for details). 

Therefore by exciting the substrate plasmons one can probe the QCB characteristics 
jHIl E2 • Indeed, substrate-QCB interaction substantially changes the conventional 
picture of substrate dielectric losses. Due to such interaction, new regions of Landau 
damping appear. The existence of these regions themselves, as well as their structure and 
the density of losses are sensitive both to the QCB period a and the direction of the wave 
vector k of the initial wave. Thus, dielectric losses in QCB-substrate system serve as a 
good tool for studying QCB spectral properties. 

The structure of this Chapter is as follows. In Section 15. 2| we briefly describe double 
square QCB interacting with the dielectric substrate and introduce the necessary defi- 
nitions. Dielectric properties of the system considered are studied in Section f5.3l where 
Dyson-type equations for the polarization operator are obtained and analyzed. The de- 
tailed description of new regions of Landau damping is presented in Section 15.41 In the 
Conclusion Section we summarize the results obtained. 
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5.2 Quantum Crossbars on Semiconductor Substrate 



Let us consider a square QCB on a semiconductor substrate (see Fig. 15.1)1 . We choose 
coordinate system so that 1) the axes Xj and the corresponding basic unit vectors ej are 
oriented along the j-th array (j = 1,2); 2) the x 3 axis is perpendicular to the QCB plane; 
3) the X3 coordinate is zero for the second array, — d for the first one, and — (d + D) for 
the substrate. The basic vectors of the reciprocal superlattice for a square QCB are Qei.2, 
Q = 27t/a so that an arbitrary reciprocal superlattice vector m is a sum m = mi + ni2, 
where m,, = rajQej (mj integer). An arbitrary vector k = fe^ei + k 2 e 2 of reciprocal space 
can be written as q + m where q belongs to the first BZ |g 12 | < Q/2. 




Figure 5.1: QCB on a substrate. e M (/x = 1,2,3) are basic vectors of the coordinate 
system. The vector ei (e2) is oriented along the first (second) array. The inter-array 
distance is d and the distance between the substrate and the first (lower) array is D. 



5.2.1 Substrate Characteristics 

The substrate is described by the Hamiltonian 

H S = H K + H C , (5.1) 

where 

qm 

is a kinetic energy of the substrate electrons with effective mass m and quadratic dispersion 
law (we omit the irrelevant spin variables), and 

1 tt t 1 v— v t TT 2vre 2 

H C = g ^q+mPq+mPq+m, Pk = y ^ C lc' C k+k', ^k = -T~ , k = q + HI, 

qm k' 

is Coulomb interaction within the substrate. 

Dielectric properties of the substrate per se are described by its dielectric function 
e s (k,w), 

' l + U k Il s {k,u). (5.2) 



e s (k,w) 



56 



Within the RPA approach, the polarization operator II s (k, u) is approximated by the 
Lindhard expression 



with 



(iyk^))- 1 = (iyk^))- 1 -^, 



no(k,a;) - fa-iew-e^+iO 



(5.3) 



(5.4) 



at T = 0. 

Active branches of substrate excitations are the surface density fluctuations which 
consist of 2D electron-hole pair continuum and surface plasmon mode with dispersion 
law 



oo s (k) = vpk\ / 1 + 



/i 2 



2krg me" 
The RPA spectrum of surface excitations is shown in Fig. 15.21 



0.25 



Ikl. 




0.05 0.1 0.15 



Figure 5.2: Dispersion of the substrate plasmons (upper line) and quasi-continuum spec- 
trum of electron-hole excitations, (dashed area). Frequency is measured in u = vf/tb 
units. 

In the case of GaAs, the substrate parameters are m = 0.068m , m Q is the free 
electron mass, vf = 8.2 ■ 10 6 cm/sec, vb = 0.78 nm, ujq = vf/tb = 1-05 x 10 14 sec -1 . For 
k < k* w 0.104r^ 1 , the plasmon frequency lies above the continuum spectrum of electron- 
hole pairs and the substrate plasmons are stable. Besides, one may easily satisfy the 
resonance condition for the collective plasmon mode near the stability threshold k ~ k* 
and QCB excitations with frequency u ~ u* = 2.6-10 13 sec -1 . For large enough k > k* the 
plasmon dispersion curve lies within the quasi-continuum spectrum and plasmons become 
unstable with respect to decay into electron-hole pairs (Landau damping of the substrate 
plasmons). Dielectric losses of an isolated substrate are described by an imaginary part 
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£ye s (k, u) of its dielectric function ()5.2|) . This imaginary part is non-zero within the dashed 
region in Fig. 15.21 due to the appearance of an imaginary part 



sn (k,c) "' 1 



2nh 2 n 2 



d{K 2 - u 2 + )^k 2 -u\ -$(k 2 - v 2 _)J~K 



i_ v 2_ 



k o i 00 

K-- — , V± = V ± K /2, V 



kp vpkF 
of the bare polarization operator n (k, oj) ()5.4|) . 

5.2.2 Interaction 

Interaction between QCB and the substrate is a capacitive coupling of charge fluctuations 
in the substrate with collective modes in the quantum wires. Assuming the distance D 
between the first array and the substrate to be much smaller than the distance d between 
arrays, one can keep only the interaction between the substrate and the first array. 
The interaction between the substrate density fluctuation at the point r = (xi,^) and 
the density fluctuations located in the vicinity of the point x\ which belongs to the ri2-th 
wire of the first array is described by its amplitude W{x\ — x[,x 2 — n 2 a), where 

w(r) = V2eH( Xl /r ) 

The function C( a; i/ r o) i n the numerator describes the screening of Coulomb interaction 
within a wire or nanotube j3H]- I n the momentum representation the interaction Hamil- 
tonian between substrate and array has the form: 



V ^ mq 



where 



W k = ik lx P^- [ drWMe 1 ^ (5.6) 
V ha J 

is proportional to the Fourier component of the interaction amplitude W(r). 

Finally, the Hamiltonian of QCB interacting with a semiconductor substrate is the 
sum of Hamiltonians ()3.7|h (|5.1|) and (|5.5j) . 

H = H QCB + H s + H sl . (5.7) 
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5.3 Dielectric Function 



High frequency properties of the system at zero temperature are determined by zeroes of 
its dielectric function 

' l + U k U(k,u). (5.8) 



e(k,w) 
Here 



U(k,tu) = '- / ilh UCl 



is the polarization of the substrate interacting with QCB, Pk(t) = 

e iHt/h pke -im/n ig the 

density of the substrate electrons in the Heisenberg representation, and averaging is per- 
formed over the ground state of the Hamiltonian (|5.7j) . 

The Umklapp processes stimulated by the interaction between the substrate and the 
first array ()5.5|) as well as the interaction between arrays ()2.27|) . produce modes with wave 
vectors q+m with various inverse lattice vectors m. This necessarily leads to appearance 
of non-diagonal polarization operators 



n(q + m,q + m» = -^ J dte^( [p q+m (t)^ +m , (0)] >. 

o 

In what follows we always consider a fixed frequency u and a fixed wave vector q from 
the BZ. So the variables q and u are omitted below for simplicity. In the framework of 
RPA approach, n(m, m') satisfies the Dyson-type equation 

n(m,m') = U a {m)S mM +U a {m)W m E 1 {m li m'). (5.9) 

The first term n s (m) in the right hand side is the substrate polarization ()5.3|) of the 
isolated substrate itself, W m = W^ +m is a bare vertex (|5.fi|) which describes substrate - 
(first) array interaction, and 



E j( m j, m ') = ~jdt e ^([e^ +mj (t),pl +mf (0)}) (5.10) 



is the correlation function of the j-th array mode and the substrate plasmon. 

The Dyson equation ()5.9j) should be completed by two equations for the correlation 
functions (j = 1,2) 

Hi(mi,m') = J D?(m 1 )^H/ m ^(m,m / ) + J D?(m 1 )^<^> m H 2 (m 2 ,m , ), (5.11) 



: 2 (m 2 ,m') 



D 2 (m 2 )^<^> m S 1 (m 1 ,m , ). (5.12) 



mi 
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Here D®{rrij) (J = 1, 2) is the bare correlation function of the j-th array modes 

oo 

1 



£>°(m 7 ) = -— / cfte itJ * 



9 jt(l+m .(t),9 j)(i+ (0) 



o cu 2 — f 2 (gj + rrij)' 



and another bare vertex $ m describes the separable inter-array interaction ()3.9|) . 

Solving the system of equations ()5.9j) . (|5.11jl and (|5.12j) one obtains the diagonal 
element Il(m) = n(m, m) of the polarization operator 

[n(m)]- 1 = [Il^m)]- 1 - |W m | 2 L>(m). (5.13) 

The second term on the right-hand side of this equation describes renormalization of the 
substrate polarization operator n s (m) by interaction between the substrate and QCB. 
The factor D(m) is a renormalized correlation function of modes of the first array 

[LKm)]- 1 = [D^m^]' 1 - (w(m) + V ( mi )) . (5.14) 

The first term w(m) describes the effective interaction between the first array and the 
substrate 

w(m) = F(mi) - |^ m | 2 n s (m), F( mi ) = ^ \W m \ 2 U s (ui). (5.15) 

rn-2 

The second one f{rni) is the effective interaction between arrays 



/ . 2 



rn 2 



*mi, uj mj = v\qj + rrijQl, (5.16) 



renormalized by Coulomb interaction of array modes with the substrate plasmons, 

= V" ^ (5 17) 

mi ^WK))- 1 -^)' 

Equations (|5.13j) - (|5.17j) together with definition ()5.8|) solve the problem of elucidation 
of the dielectric properties of the combined system QCB-substrate. 

The spectrum of collective excitations in QCB-substrate system is determined by zeros 
of the dielectric function e(q, u) = 0. The key question here is the robustness of the QCB 
spectrum against interaction with 2D substrate excitations. Detailed analysis shows that 
in the long wave limit q Q the interaction just renormalizes the bare dispersion laws 
of the arrays, conserving its LL linearity. This result verifies stability of QCB plasmons 
with respect to substrate-QCB interaction. 

The QCB-substrate interaction also results in occurrence of some special lines in the 
BZ. These lines correspond to resonant interaction of the substrate with the first or the 
second array. The resonance condition w s (k) = ^j(k) is fulfilled along the line LJIN for 
j = 1 and along the line KBM for j = 2 in Fig. 15 .5| below. 
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5.4 Landau Damping 



As was mentioned in subsection 15.2.11 above, dielectric losses of an isolated substrate 
are related to the Landau damping due to decay of substrate plasmons with momentum 
k > k* « O.KMr^ 1 into electron- hole pairs. The substrate-QCB interaction remarkably 
modifies the conventional picture of substrate plasmon dielectric losses. Due to QCB- 
substrate interaction, new domains of Landau damping appear in addition to the dashed 
region in Fig. 15.21 Indeed, outside the initial instability region where 3e s (k, uj) = 0, 
nonzero imaginary part 9fe(k, uj) ()5.8|) exists if the imaginary part of the bare polarization 
operator 3n (k + m, to) differs from zero at least for one of the reciprocal lattice vectors 
m. The main contribution to 9te(k, uj) is related to the renormalization term w(m) 
in Eq. (j5.15j) due to Umklapp processes along the x 2 axis (summation over m2 in the 
expression for the function F{m\) is implied). It is proportional to the fourth power 
of QCB - substrate interaction W 4 . The Umklapp processes along both directions 
contribute also to the renormalization term <p(mi) in Eq. (|5.1fi|) . However, they contain 
an additional small parameter 4 related to inter-array interaction within QCB. These 
terms are not taken into account. Thus, the possible Umklapp vectors have the form 
m 2 = m 2 Qe 2 , m 2 = ±1, ±2, . . . and in what follows we will label them by an integer 
number m 2 . 

a/r B 

100 
80 
60 
40 
20 


0.01 0.03 0.05 0.07 0.09 k Fr B 

Figure 5.3: Phase diagram describing appearance and structure of new regions of Landau 
damping. Lines 1 — 6 separate different types of new damping regions. Points a — f 
correspond to the structures displayed below in Figs. I5.4H5.91 

The structure of the new Landau damping regions and their existence themselves is 
governed by interplay between the Fermi momentum of the substrate kp and the QCB 
period a. The first of these parameters defines the width of the two-particle excitation 
band (dashed region in Fig. 15 .2j) while the second determines the minimal reciprocal 
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vector Q. In the case of sufficiently thick QCB superlattice (small a) and sufficiently low 
electron density within the substrate (small kp), Umklapp processes are always ineffective 
because they change an initial plasmon wave- vector into the outer part of the instability 
region. This means that only plasmons with momenta \k\ > k* decay into electron- hole 
pairs. 

Increasing the QCB period or the Fermi momentum turn the Umklapp processes 
effective, and additional Landau damping regions appear within the circle \k\ < k*. The 
first factors involved are the smallest Umklapp vectors ±1, then new damping regions 
appear corresponding to the Umklapp vectors ±2 and so on. As a result one gets a 
rich variety of possible damping scenarios. We describe them with the help of a "phase 
diagram" in the a — kp plane displayed in Fig. 15 .HI (actually dimensionless coordinates 
a/rs and kprs are used). Here the set of curves labelled by numbers 1 — 6 separate 
the regions of parameters corresponding to the different Umklapp vectors and different 
structures of the new damping regions. There is no additional Landau damping regions 
below the first line. Above the sixth line Landau damping takes place within the whole 
circle k < k*. Above lines with numbers 2n — 1, the Umklapp vector ±n becomes effective. 
The corresponding additional damping region has the form of a tail touching the initial 
Landau damping region \k\ > k*. This tail turns to the additional damping band well 
separated from the initial one above the line number In in < 3) within some sector of 
directions in fc-space (in what follows, these directions will be labelled by corresponding 
arcs of the circle). 



k 2 r B 




Figure 5.4: New Landau damping region PUW for QCB with period a = 20 nm (point 
a in Fig. I5.3J1 corresponds to the Umklapp vector —I. Other details of this figure are 
explained in the text. 

Possible structures of new damping regions corresponding to some representative 
points a — / in the a — kp plane (see Fig. 15. 3 j) are displayed in detail in Figs. 15.41 
15.91 All these figures correspond to the GaAs value of kptB ~ 0.038. Generally speaking 
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Figure 5.5: New damping region STUW for a = 30 nm (point b in Fig. 15.3)1 corresponds 
to the Umklapp vector —1 and describes Landau damping tail (within the arc TU) or 
separate landau band (within the arc TP). 



we should display new damping region within the whole circle \k\ < k* in the plane ki,k 2 - 
The circle center T is placed at the origin (we did not put the letter T in Figs. l5.4H5.9j) . 
But this region is always symmetric with respect to reflection k% — > — k\ and with respect 
to the combined reflection k 2 — > — k 2 , m 2 — > —m 2 . This enables us to describe the damping 
scenarios only within the quarter k^ 2 > (complete picture of the new damping region 
can be easily obtained from the displayed one with the help of the reflection symmetries 
mentioned above). 

Damping of the substrate plasmon occurs inside the arc PR in Figs. 15 .4115. 91 when at 
least one of the points in the phase space with coordinates (k + m 2 ,uj s (k)) lies within 
the quasi-continuum spectrum of the electron-hole excitations, whereas the "mother" 
point (k,u s (k)) lies above the continuum (above the dashed area in Fig. 15. 2j) . As was 
mentioned above, for small enough QCB period a < a± = 17.3 nm, the basic reciprocal 
lattice vector Qe 2 is too large, the points (k + m 2 , w s (k)) lie outside the quasi-continuum 
for all 777-2 and additional Landau damping region does not exist. It appears only for 
a > at (Q < 2k* + 2k F ). For a x < a < a 2 = 23.6 nm (2k* + 2k F > Q > 2k*), this is the 
region PUW (see Fig. I5.4J) corresponding to the 777.2 = — 1 (in all Figs. 15 .4115 .91 the regions 
related to this Umklapp vector are always hatched by the hatching tilted to left). As a 
result, the damping tails touching the initial Landau damping region appear in certain 
directions of the k plane. 

For a 2 < a < 03 = 34.6 nm (2k* > Q > k* + kp), the new damping region is related 
to the same Umklapp vector —1, but now it has a strip-like structure bounded by the 
line STUW in Fig. 15.51 Note that the damping is absent within the region PTS. As a 
result, it is possible to divide the angular region < (p < it/ 2 into three sectors. Within 
the first one PT, < (p < ^3, a new damping region is separated from the initial one. 
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The second sector TU, (f% < <p < <^6> corresponds to a new damping tail. Finally, within 
the third sector UR , ip s < ip < tt/2, new damping region does not exist at all. 




Figure 5.6: In the case a = 40 nm (point c in Fig. I5.3|) . new Landau damping regions 
PVD, STRUW, and UHR correspond to the Umklapp vectors —2, —1, and +1,. 



k 2 r B 




Figure 5.7: New regions of Landau damping for a = 50 nm (point d in Fig. 15. 3J) . Regions 
ACVD, STR, and FHR correspond to the Umklapp vectors —2, —1, and +1 

For larger QCB period, CI3 < a < = 47.1 nm (k* + kp > Q > k*), the new damping 
regions have a more complicated structure. In fact the damping area consists of three 
parts (see Fig. 15. 6J) . The first one, STRUW corresponds to the Umklapp vector — 1. 
Note that it is shifted to the bottom with respect to the previous case a = 30 nm. This 
part overlaps with the second part HRU . The latter corresponds to the Umklapp vector 
+ 1 and is hatched in Figs. I5.6H5.9I by horizontal hatching. The second region PDV 
corresponds to the Umklapp vector —2 (in Figs. l5.6H5.9l such a regions are always hatched 
by vertical hatching). As a result in the direction close enough to the &2 axis, one gets a 
new damping tail with m 2 = —2) and well separated new damping band with Umklapp 
vector —1. 

Further increase of QCB period < a < = 52.3 nm, k* > Q > 2(k* + kp)/3 
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(k* pa 2.8kp for GaAs) leads to further extension of new damping regions. The region 
ADVC corresponding to the Umklapp vector —2 is partially separated from the initial 
damping region. It overlaps with the region FHRT [m^ — +1) which in its turn overlaps 
with the region STRT (m,2 = — 1)- Actually the two latter regions do not include an 
extremely small vicinity of the origin T which is not shown in Fig. 15.71 Visible coincidence 
of the points V and T in Fig. 15.71 is an artefact of the accuracy of the figure. 




Figure 5.8: New regions of Landau damping for a = 60 nm (point e in Fig. 15. 3|) . Besides 
the regions ACVD, STR, and FHR corresponding, as in the case a = 50nm, to the 
Umklapp vectors —2, — 1, and +1, new Umklapp vector —3 appears (region PEG). 




Figure 5.9: The case a = 70 nm (point / in Fig. 15. 3 J) . New regions of Landau damping 
PEG, ACVD, STR, whole sector TPR, and VOR correspond to the Umklapp vectors 
-3, -2, -1, +1, and +2. 

Within the next interval of QCB periods < a < = 65.2 nm (2(k* + fcp)/3 > Q > 
2k f)) new damping region GPE corresponding to the Umklapp vector —3 appears. This 
region is hatched in Figs. 15.81 and 15.91 bv the hatching tilted to the right. Beside that, the 
regions ADVC (m 2 = —2), FHRT (+1), and STRT (—1) are present. As in the previous 
figure, visible coincidence of the points E and C in Fig. 15.81 is an artefact of the accuracy 
of the figure. 
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Finally for a > ciq (Q < 2kp), Landau damping emerges in the whole circle \k\ < k* 
(Fig. 15. 9|) This occurs due to processes with Umklapp vector +1. The corresponding 
damping region covers the whole quarter. Therefore we did not hatch it at all, and used 
the same horizontal hatching for the new region VOR corresponding to the Umklapp 
vector +2. The region DACRV is related to the Umklapp vector —2. We emphasize that 
at the same time the vertex V of this region is the vertex of the region VOR. This is not 
an accidental approximate coincidence as in the two previous figures. The regions GPE 
and TSTR are related to the Umklapp vectors —3 and —1 respectively. 
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Figure 5.10: Damping tail for <p w 16°. Precursor of the resonant peak and the resonant 
peak are resolved quite well. 

Thus, the general structure of the additional damping regions is described in Figs. 
15 .4115. 91 However, there is an additional structure of these regions. This fine structure is 
related to possible resonance interaction between the substrate plasmons and the QCB 
plasmons of the first or second array. The resonance condition for the first (second) array 
is written as u s (k) = Ui(k) (u s (k) = w 2 (k)). 

hi e 




0.0285 0.0286 0.0287 0.0288 0.0289 



Figure 5.11: Damping tail for tp « 33°. Precursor of the resonant peak and the resonant 
peak are resolved quite well. 

Consider such fine structure of new Landau damping region in details for QCB with 
period a = 30 nm. Here the resonance conditions are satisfied along the lines LJIN 
and KBM (see Fig. 15.5)1 . These lines intersect with the damping region boundaries 
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at the points J, I, B and MB which define four rays 01, OJ, OB, and OU and four 
corresponding angles <p\ ~ 16°11', <p 2 ~ 33°19', <^4 ~ 69°, ip 5 ~ 76°. The resonance 
interaction takes place within two sectors (p% < (p < ip 2 and (p^ < ip < (p§. 

For each (p < (pi the new damping region is a well separated damping band. The 
damping amplitude is small because of the small factor of order W 4 mentioned above. 
When (p — > — 0, small peak appears near the "blue" boundary of this damping region. 
This peak is a precursor of the resonance between the substrate plasmon and the first 
array QCB plasmon (Fig. 15.10)1 . The same happens from the opposite side of the sector 
(tpi, (p 2 when ip —> ip 2 + (Fig. 15.11)) . 

Within the sector ip x < <p < ip 2 , the damping band contains a well pronounced peak 
corresponding to resonant interaction between the substrate plasmon and the first array 
plasmons (see Figs. 15.101 and 15.11)) . The peak amplitude is of order of the damping 
amplitude within the initial damping region. It has a Lorentz form placed on a wide 
and low pedestal. The peak is especially sensitive to the strength of the QCB-substrate 
interaction which is governed by the distance D between QCB and substrate. 

To study this D dependence, let us consider the imaginary part Se(k, uS) of the di- 
electric function within the considered sector ipx < (p < (p 2 . In the vicinity of the plasmon 
frequency uj ~ uj s (k) this imaginary part is written as 

cwv ^ |W k | 2 -$sw(k,u a (k)) 

U k ( W 2 _ ^2)2 + (g w ( k>Ci;s (fc))) 2 ' 

with u>(k, oj s (k)) being of order \W\ 2 . So the resonance peak indeed has the Lorentz like 
shape with height of order unity, 
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whereas its half-width 

r = ^(k,t;|A; 1 |)~|^ k _Q e2 | 2 , 

is of order W 2 . The peak is displayed in Fig. 15.121 for different values of the distance D 
between the substrate and the nearest (first) array. It is seen that the amplitude changes 
slowly with increasing distance D while its width squeezes sharply, W 2 ~ 1/D 2 . 

There is no resonance interaction within the sector ip 2 < (p < ip^ but further increase 
of the angle v?4 < ip < ip$ leads to re-appearance of the resonant peak within the damp- 
ing tail. In this case one deals with a resonance between the substrate plasmon and the 
QCB plasmon in the second array. Existence of this resonance is caused by inter- array 
interaction that brings additional small parameter to the imaginary part of the dielec- 
tric function. As a result, the width of the peak is much smaller in the second sector 
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Figure 5.12: Damping tail for tp = 20° for different distances D between QCB and sub- 
strate. The curves 1, 2, 3, and 4 correspond to D — 1 nm, 1.5 nm, 2 nm, and 2.5 nm, 
respectively. With increasing D, the resonant peak narrows and slowly increases, whereas 
the area under the curve decreases. 

while surprisingly, the peak amplitude has the same order of magnitude as in the case of 
resonance with the first array (closest to the substrate). 

Existence of the additional QCB bands (tails) of Landau damping and appearance of 
the resonant peaks within the bands (tails) is a clear manifestation of interplay between 
real 2D surface plasmons and quasi-2D QCB plasmons. 



In conclusion, the possibility of spectroscopic studies of the excitation spectrum of quan- 
tum crossbars interacting with semiconductor substrate is investigated. A capacitive 
contact between QCB and substrate does not destroy the LL character of the long wave 
excitations. However, the interaction between the surface plasmons and plasmon-like ex- 
citations of QCB essentially influences the dielectric properties of the substrate. The QCB 
may be treated as a diffraction grid for the substrate surface, and Umklapp diffraction 
processes radically change the plasmon dielectric losses. Due to QCB-substrate interac- 
tion, additional Landau damping regions of the substrate plasmons appear. The structure 
of these regions and the spectral density of dielectric losses are strongly sensitive to the 
QCB period. The surface plasmons are more fragile against interaction with superlattice 
of quantum wires than the LL plasmons against interaction with 2D electron gas in a 
substrate. 



5.5 Conclusions 
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Chapter 6 

Ultraviolet Probing of Quantum 
Crossbars 

6.1 Introduction 

The IR based methods mentioned above are not very convenient from two points of view. 
First, as it was mentioned, one needs an additional diffraction lattice to tune the light 
wave vector and that of the QCB plasmon. Second, they probe QCB spectrum only in 
some discrete points. The alternative method of studying QCB spectrum by ultraviolet 
(UV) light scattering is the subject of the present Chapter. The advantages of this method 
are evident. It does not require any additional diffraction lattice. It probes QCB spectrum 
in a continuous region of wave vectors. Finally, its selection rules differ from those for IR 
absorption. This gives rise to the observation of additional spectral lines not visible in IR 
experiments (S51EE]. 

In this Chapter we formulate the principles of UV spectroscopy for QCB and study 
the main characteristics of scattering spectra. The Chapter is organized as follows. In 
Section 16.21 we discuss light scattering on QCB and present basic equations describing 
this process. The main results of the Chapter are contained in subsection 16.21 where 
we classify the basic types of the scattering indicatrices (angular diagrams of differential 
cross section) corresponding to various detector orientations. The results obtained are 
summarized in the Conclusion Section. Technical details are concentrated in Appendix 
IE1 This Appendix is devoted derivations of an effective QCB-light interaction and the 
basic formula for differential cross section of light scattering. 
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6.2 Light Scattering on QCB 



The simplest process contributing to the Raman-like light scattering is an annihilation 
of an incident photon and creation of a scattered photon together with a QCB plasmon 
(Fig. 16.1b ). In terms of initial electrons, this is in fact the second order process. Since 
the energies of incident and scattered photons significantly exceed the electron excitation 
energy in nanotube, one may consider the emission/ absorption process as an instantaneous 
act (see Appendix |E]) . This process may be treated as the inelastic photon scattering 
accompanied by emission or absorption of a plasmon (Fig. 16.1b ) . 



Figure 6.1: (a) The scattering process geometry (all notations are explained in the text), 
(b) Second order diagram describing light scattering on QCB. Solid lines correspond to 
fermions, whereas dashed lines are related to photons. The vertices are described by the 
interaction Hamiltonian (jfj.lj) . (c) Effective photon-plasmon vertex. QCB excitation is 
denoted by wavy line. 

Let Kj (Kf) and Qi = cKi = c|Kj|, (Qf = cKf) be momentum and frequency of the 
incident (scattered) photon. For simplicity we restrict ourselves by the case of normal in- 
cidence Kj = (0, 0, K,j). Direction of the scattered wave vector is characterized by the unit 
vector n((p,i)), where (p and $ are polar and azimuthal angles of the spherical coordinate 
system with polar axis oriented in the e3 direction. The momentum of the excited QCB 
plasmon is k, the projection of the scattered photon momentum onto QCB plane is — k. 
This process is displayed in Fig. 16.1b . If the QCB plasmon frequency u;(k) coincides with 
the photon frequency loss to = Qi — Qf, then a detector oriented in a scattered direction n 
will register a sharp well pronounced peak at the frequency loss uo = uj(k). The frequency 
loss is much smaller than the incident and scattered photon frequencies. Therefore in 
what follows we use the same notation K for both Ki and Kf where it is possible. Scan- 
ning the frequency loss u, (or, that is the same, modulus of the scattered wave vector 
Kf) for a fixed detector orientation n, two or more (up to six for a square QCB) such 
peaks can be observed. The number of peaks and their location strongly depend on the 
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azimuthal angle § in the QCB plane. Scanning this angle, one can change the number of 
observable peaks. This is yet another manifestation of dimensional crossover mentioned 
above. 

An arbitrary vector k = (ki, k 2 ) of reciprocal space can be written as k = q+m, where 
q = (qi, q 2 ) belongs to the first BZ. QCB eigenstates are classified by quasi-momentum q 
and 2D band number. However, the specific QCB geometry makes its spectral properties 
rather unusual. Consider an isolated array 1. As was mentioned in Section E21 within the 
(xi, x 2 ) plane, its excitations are characterized by a pair of 2D coordinates (xi, nza), i.e. a 
continuous longitudinal coordinate X\ parallel to e x direction, and its discrete transverse 
partner n 2 a parallel to e 2 . As a result, the longitudinal component k\ = q\ + rriiQ of 
the excitation momentum changes on the entire axis —00 < k\ < 00 while its transverse 
momentum q 2 is restricted to the region \q 2 \ < Q/2. Thus an eigenstate (plasmon) of 
the first array is characterized by the vector k x = q + ni! = (ki,q 2 ) and the frequency 
o;(ki) = v\ki\ which depends only on the longitudinal component k\ of the momentum 
ki. Similar description of the second array is obtained by replacing 1^2. 

The exact differential cross section of the scattering in a given direction, considered 
as a function of the scattered frequency, has almost equidistant sets of peaks. However, 
because of weak inter-array interaction ()3.6|) . only few peaks are well pronounced, whereas 
other peaks are very low and will be omitted below. Therefore in the approximation 
adopted here the scattering cross section is characterized by a number of peaks, their 
positions and intensities. Calculation of these quantities as functions of the frequency 
loss to and azimuthal angle for various fixed values of the polar angle (p, is the goal of 
our study. 

To perform quantitative analysis, we should derive an expression for the scattering 
cross section. All the details of our calculations are contained in Appendix |E| Here we 
present only the main steps of the derivation. We start with the Hamiltonian 



which describes interaction between a single nanotube oriented, e.g., along the ei direction 
(the first array nanotube) and an external electromagnetic field. The field is described by 
its vector potential in the Landau gauge 



The indices a = A, B enumerate sublattices in a honeycomb carbon sheet, (r, 7) are polar 
coordinates in the (x 2 ,x 3 ) plane, ip a (xi,j) = ip a (xi, To, 7) and ip^ are slowly varying 
electron field operators at the nanotube surface r = ro, and the vector of Pauli matrices 




(6.1) 



A = A&x + A 2 e 2 + A 3 e 3 . 



(6.2) 
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a is cr = &\O x + e 7 cr y , e 7 = — e 2 sin 7 + e 3 cos 7. The light wavelength is much longer than 
the nanotube radius, so the vector potential can be taken at its axis, A(xi, 0, 0). 

Such form of a nanotube- light interaction leads to the following expression for an 
effective QCB-light interaction Hamiltonian 

Hm = -^JT \) ^ / dxid Xl 9i(xi,n2a)Al(xi,n 2 a, 0) + 

+ 4t|~ (~) ^ f dx2d X2 9 2 ( y nia,x 2 )Al(n 1 a,X2,d). (6.3) 

Here Aj(r) = A(r) + (v^2 — l) Aj(r)ej (j = 1,2) are two effective vector potentials 
affecting two arrays and Aj(r) are two (of three) cartesian components of the full vector 
potential ()6.2j) . The main subject of our interest is the scaled scattering differential cross- 
section cr(u;,n) defined by equation 

*, = ALV(„,n)^ A = 3 JL-(Q 2 (H£)\ (6 . 4) 

where L 2 is the QCB area and do = sin (pdipdi9. A standard procedure applied to the 
Hamiltonian ()6.3|) . leads to the following form of the scaled cross section 

n ) = 47 E K p l h l°) P l 2 ^ - ( 6 - 5 ) 

p 

Here = Qi — Qf is the frequency loss and 

h = " E -jt\ P ^f^ + 4- kj ) (6-6) 

is the interaction Hamiltonian reduced to the subspace of QCB states and summation 
is performed over all one-plasmon states \P). Pj-x ft Xi(v,^), j = 1,2, are polarization 
matrices. In the basis (||, _L) they are given by 



/ zsmv 1 cos v \ , / 7r \ , 

Pi{<P,&) = -[ . ) , P 2 (^) = Pi(-+tf,<A (6.7) 

\ 1% cos it cos (p sin it cos (p J \z / 

Equations ()6.5|) - ()6.7|) serve as a basis for the subsequent analysis. 



6.3 Scattering Cross Section 
6.3.1 Cross Section: Basic Types 

According to Eqs. ()6.5j) . ()6.6j) . in order to contribute to the cross-section (16 .5|) for a fixed 
detector orientation at n, an excited QCB plasmon \P) must contain at least one of two 
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single-array states |1,— ki) or |2, — k 2 ). Analysis shows that there are five basic types 
of excited QCB plasmons depending on the location of the point k. Here we describe 
these types of plasmons and the corresponding structure of the differential scattering 
cross section. In this description beside the polar angle (p of the scattered photon wave 
vector K.f, we use in a sense a mixed representation. It is based on the excited plasmon 
momentum k and the azimuthal angle $ of the transverse component of the scattered 
photon wave vector K/. 




Figure 6.2: Part of inverse space. The small square kj/Q < 0.5 is the quarter of the first 
QCB BZ. High symmetry lines are parallel to the coordinate axes. Resonant lines are 
parallel to BZ diagonals. The arcs coming over the point Bi, B 2 , B 3 , B4, B 5 , correspond 
to different wave numbers of excited plasmons: \k/Q\ = 0.3; 0.5; 0.6; 0.7; 0.8. 



i. The general case: The point k lies away from both the high symmetry 
lines and the resonant lines. This case is illustrated by the point A in Fig. 16.21 

There are two QCB plasmons \P), mostly propagating along the first or second array 
with the frequencies ()B.2j) . which contribute to the scattering. The differential cross- 
section is a sum of two peaks centered at these frequencies. After averaging over initial 
polarizations and summing over final polarizations, it has the form a(u, n) = F 1 (n)<5(c<j — 
uJi)+F 2 (n)5(uj— u> 2 ), where the functions F!(n)=F(y9, $) and F 2 (n)=F(<p, r d+n/2) describe 
universal angle dependence of the peak amplitudes, 

r 3 l 1 

F((p, $) = \ sxmp cos i?| 1 — - sin 2 cos 2 $ + - cos 2 <p . 

The functions F± j2 are related to the corresponding array plasmons. Each one of them 
vanishes when the scattered photon is perpendicular to the corresponding array. However, 
these functions describe the scaled cross section. The absolute value amplitude of each 
peak has an additional factor A (see Eq. ()6.4J0 . Strong electron-electron interaction 
in a nanotube corresponds to small values of the Luttinger parameter g and therefore 
suppresses the scattering cross-section. 
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ii. Inter-band resonance in one of the arrays: the point k lies on a high 
symmetry line of only one array. This case is illustrated by the points 6*2,3 (1-st 
array) and -D2.3 (2-d array) in Fig. 16.21 Consider for example point C2 where k\ = Q/2, 
^2 7^ nQ/2. Here three QCB plasmons contribute to the scattering. The first one is |2, k%) 
plasmon propagating along the second array with quasimomentum fc 2 and frequency 002, 
Eq. (jB.2|) . The two others are even or odd superpositions of the 1-st array states (Eq. 
fFTTQj) with j = 1) of the two first zones with eigenfrequencies (jB.lljl . Due to weakness 
of the inter-array interaction, three peaks of the scattering cross section form a singlet 002 
and doublet U\ g , uj\ u . After averaging over initial and final polarizations, the cross section 
has the form u(u, n) = 2 _1 J?. Fj(n)[6(u — u>j g ) + S(u — u)j U ))]- 

iii. Inter-band resonance in both arrays: The point k is a crossing point 
of two high symmetry lines away from all resonant lines. This case is illustrated 
by the points C 4 and D4 in Fig. 16.21 Consider for example point C 4 . Here ki = Q/2, 
k 2 = 2Q/2, and four QCB plasmons contribute to the scattering. The first pair consists 
of even and odd superpositions of the 1-st array states of the first and the second bands. 
These states and their frequencies are described by Eqs. (|B.14|) . (|B.15|) with j = 1. The 
second pair consists of the same superpositions of the 2-d array states from the second and 
third bands and is described by the same equations with j = 2. As a result, four peaks, 
which form two doublets (|B.15|) . j = 1,2, can be observed. After averaging over initial 
and final polarizations the cross section is ct(uj, n) = 2 _1 Fj(n)[S(uj — ujj g )+S(uj — ujj U ))]. 

iv. Inter-array resonance: The point k lies only on one of the resonant 
lines away from the high symmetry lines. This case is illustrated by the points 
i?i_3 5, E, and F in Fig. 16.21 Here the QCB plasmons which contribute to the scattering 
are two even and two odd superpositions of the first and second array states (|B.12|) whose 
eigenfrequencies form two doublets (jB.13|) . As in the previous case, the scattering cross 
section contains four peaks which form two doublets. After averaging over initial and final 
polarizations the cross section is a(u, n) = 2 _1 V. Fj(n)[5(uj — u) g (\c.j)) + 5{u — uj u (k.j))]. 

Thus, the inter-array splitting is proportional to the main small parameter (jrS.8)l of the 
theory, <p w 0.007 jSHl- For the set of parameters described in the beginning of subsection 
13.51 the inter-band splitting defined by Eq. (jB.ll|) . is five times smaller because it contains 
an additional factor (fra/rQ. 

v. Inter-array and inter-band resonance: The point k lies at the intersection 
of two resonant lines. There is only one such point B4 in Fig. 16.21 In the general case, 
where the parameter n 7^ for both crossing resonant lines (the point B 4 is not the 
case), the QCB plasmons involved in Raman scattering form two quartets. The first 
quartet consists of four symmetrized combinations (jB.lfij) of the single-array states. QCB 
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eigenstates for the second quartet are obtained from these equations by replacing 1 «-> 2. 
The corresponding eigen-frequencies are described by Eqs. (jB.17|) . (|B.18[) . The scattering 
cross section in this case contains six peaks, and two of them are two-fold degenerate 

a(u,n) = -^F i (n)[2%-a; S)U/s (k^) + %-a; M (k j )) + %-a; us (k J -)]. 

i 

The point B± in Fig. 16.21 lies on the main resonance line with n = 0. Here \ki\ = \k 2 \, 
the frequencies of both quartets coincide, and the scattering cross section contains one 
four-fold degenerate peak and a symmetric pair of its two-fold degenerate satellites. 

This classification of all types of excited plasmons enables us to describe completely 
the UV scattering on QCB. 

6.3.2 Scattering Indicatrices 

In this part we describe the results of the scattering process with the help of a family of 
scattering indicatrices. 

To explain the indicatrix structure, we start with some preliminary arguments. Con- 
sider the case where the detector is tuned on the frequency Q f and is oriented in direction 
n. These parameters determine a point k by a unique way. There are two ways of scanning 
QCB plasmons. The first way is to scan the polar angle The corresponding points k 
in Fig. 16.21 form a circle with radius k = (flf/c) sin (p. The second way is to tune the 
detector frequency as Qf. We are interested in frequency loss of order of the plasmon 
frequencies within the two - three lowest bands. This loss is of order of Qf(v/c) <C Qf. 
Therefore in this case the point k remains on its place with a very good accuracy. 

Each scattering indicatrix corresponds to a circular arc in Fig. 16.21 and the structure 
of the indicatrix is completely determined by the arc radius k. The indicatrix represents 
a set of curves displayed in polar coordinates with polar angle, which coincides with the 
polar angle $ used above, and the (dimensionless) radius u/(vQ), where u is the frequency 
loss. Each point of the indicatrix corresponds to an excitation of a QCB plasmon and 
therefore to a sharp peak in the scattering cross section. The number of peaks depends 
on the polar angle. Scanning the polar angle $ results in changing the number of peaks. 
This is one more example of dimensional crossover in QCB (see Chapter 0] for similar 
effects in IR spectroscopy). 

We start with the case of the smallest radius K sin if = — 0.3Q (arc AB\ in Fig. 16. 2|) . 
Here all points beside the point Bi are points of general type (i). Each one of them, e.g., 
point A, corresponds to excitation of two plasmons in the two arrays and therefore leads 
to two separate peaks in the scattering spectrum (see Fig. 16 .H| left panel). The peaks 
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corresponding to the point A in Fig. 16.21 lie on the ray defined by the angle 9 a- The point 
Bi corresponds to an inter-array resonance (iv) and in this direction a split doublet can 
be observed. 




Figure 6.3: Left panel: Positions of the scattering peaks for \k/Q\ = 0.3. The doublet in 
the resonance direction (point B\ in Fig. 16. 2|) is well pronounced. Right panel: Positions 
of the scattering peaks for \k/Q\ = 0.5. Two doublets appear at the BZ boundaries (points 
d and Dt in Fig. IfHjl . 

In the next case Ksimp = —0.5Q (arc C\B 2 D\ in Fig. I6.2J) . as in the previous one, 
in the general directions one can observe two single lines which form a split doublet in 
the resonant direction B 2 . However at the final points C\ and D\ the arc touches the 
high symmetry lines. Here the low frequency single line vanishes (there is no scattering at 
k = 0) while the high frequency line transforms into a doublet because of an inter-band 
resonance (ii) in one of the arrays (Fig. 16 .3| right panel). 

Further increase of the arc radius K simp = — O.QQ (arc C 2 B%D 2 in Fig. 16. 2j) leads to 
appearance of two points D 2 and C 2 where the arc intersects with high symmetry lines 
(BZ boundaries). Each of these points generates an inter-band resonance doublet, which 
coexists with the low frequency single peak (Fig J6.4[ left panel). 

In the case if sin = —\^2Q/2 the corresponding arc includes the BZ corner B±. This 
is the point of a double inter-array and inter-band resonance (v). Moreover, here, the 
two quartets described above coincide. Therefore, there are three lines in Fig. 16.41 right 
panel. The low-frequency line as its high-frequency partner is two-fold degenerate while 
the central line is four-fold degenerate. We emphasize that each quartet manifests itself 
in three lines, contrary to the IR absorption, where selection rules make two of them 
invisible. 

The last case K sin if = — 0.8Q demonstrates one more possibility related to inter- 
band resonance simultaneously in two arrays (iii). Each point E and F generates (in 
the corresponding direction) two doublets describing the inter-band splitting in different 
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Figure 6.4: Left panel: Positions of the scattering peaks for \k/Q\ = 0.6. Two doublets 
at the BZ boundaries (points C2 and D2) in Fig. I6.2)l are shifted from the high symmetry 
directions. Right panel: Positions of the scattering peaks for \k/Q\ = y/2/2. Resonance 
triplet corresponding to point B 4 in Fig. I6.2I (two of four frequencies remain degenerate 
in our approximation). In IR experiments only one of the triplet components is seen. 



arrays (see Fig. I6.5|) . 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 <»/vQ 



Figure 6.5: Positions of the scattering peaks for \k/Q\ = 0.8. Two pairs of doublets appear 
corresponding to excitation of two pairs of plasmons in the two arrays (points E and F 
in Fig. Q. 



6.4 Conclusions 

In conclusion, we studied inelastic UV Raman scattering on QCB. We derived an effective 
Hamiltonian for QCB-light interaction which is expressed via the same Bose variables that 
the QCB itself. With the help of this Hamiltonian we calculated the differential scattering 
cross section as a function of detector orientation and scattered frequency. Scanning these 
parameters, one can observe a set of sharp peaks in the scattering spectrum. The number 
of peaks and their positions strongly depend on the direction of the scattered wave vector. 
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This results in a dimensional crossover. It manifests itself in the splitting of the peak 
frequencies and therefore in appearance of multiplets (mostly doublets) instead of single 
lines in the scattering spectrum. 

The sizes of peak splitting are determined by the nature of interaction which lifts 
the corresponding degeneracy. In the case of initial inter-array degeneracy, the splitting is 
proportional to the dimensionless interaction strength cf) « 0.007. The inter-band splitting 
is proportional to the square of this parameter. For chosen QCB parameters it is smaller 
than the inter-array interaction strength in spite of an additional large multiplier a/ro. In 
all cases, the splitting increases with increasing the interaction in QCB crosses. The peak 
amplitudes are proportional to the Luttinger parameter g in a single nanotube. Therefore 
strong electron-electron interaction suppresses the peaks. 

The effectiveness of UV scattering is related to the possibility of changing continuously 
the excited plasmon frequency. Due to other selection rules, some lines which are invisible 
in the IR absorption spectrum, become observable in the UV scattering. UV scattering 
spectroscopy enables one to restore parameters describing both interaction in QCB crosses 
and electron-electron interaction in a single QCB constituents. 

Our studies of optical properties of QCB ( this Chapter and Chapters^ and{Ej) show that 
these nanoobjects possess a unique combination of optical spectra. Firstly, they are active 
in IR and UV frequency range. Secondly, they may be observed in various kinds of optical 
processes, namely direct and indirect absorption, diffraction, energy loss transmission, and 
Raman-like spectroscopy. 
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Appendix A 
Empty Super-Chain 



Here we construct eigenfunctions, spectrum, and quasi-particle operators for an "empty 
super-chain" - quantum wire in an infinitely weak periodic potential with period a. Exci- 
tations in an initial wire are described as plane waves L~ 1//2 exp(ifcr) with wave number 
k = 2nn/L, with integer n, and dispersion law u{k) = v\k\ (the array number is tem- 
porarily omitted). The following orthogonality relations are valid 

f L/2 \^ 

J -L/2 k 

where 8l stands for periodic delta-function 

5 L (x — x') = 5(x — x' — nL). 

n 

An "empty super-chain" is characterized by a space period a and a corresponding re- 
ciprocal lattice wave number Q = 2n/a. Each excitation in such a super-chain is described 
by its quasi-wavenumber q and a band number s (s = 1,2, . . .) that are related to the 
corresponding wave number k by the following relation, 



k = q + Q{-\) 



s-l 



S 

L2J 



sign q. 



Here, square brackets denote an integral part of a number. The corresponding wave 
function ip s>q {x) has the Bloch-type structure, 

iM*) = -^e^u s , q (x), (A.l) 

and satisfies the orthogonality relations 

"L/2 



^*s, q (x)^ s ', q '(x)dx = 5 S , S , V 5 q+mQiq/ , V i)* (x)i> SA (x) = 5 L (x - x'), 

-L/2 V 



s,q 



79 



Within the first BZ, — Q/2 < q < Q/2, Bloch amplitude and dispersion law ui s have the 
following form 



u s , q ( x ) = ex P {iQx(-l) s 1 I sign q j , w a (g) = ^ 



+ (A.2) 



Taking into account that both Bloch amplitude u s>q (x) and dispersion law u> s (q) are peri- 
odic functions of q with period Q, one obtains general equations for the Bloch amplitude 

sin£ n 



u 



t ■ cos [(2s - l)f n ] exp y-^tn^j 



4^n = Q(x-na), 



2s -1 ^ 2(-l) 



27r(2n + l)q 
cos - 



and dispersion law cu s (q) 

u s (q) 

vQ " 4 ^^7r 2 (2n+l) 2 ~" Q 

The relations between quasiparticle operators for a free wire, for momentum k ^ nQ/2 
with n integer, and those for an empty super-chain, C S)(? , for quasimomentum q from the 
first BZ, -Q/2 <q< Q/2, look as 



c fc = C^sign*;, 



1 + 



2\k[ 

Q 



c. 



s,q 



:-!)%, k = q + (-l)-Q 



z/ = s + l + 



2g 



where curely brackets denote a fractional part of a number. For obtaining these relations 
we used the following expression 



f L/2 \^ 

/ *l>k( x )*l>s, q (x)dx = 5 sAq) signk \ 5 q+mQ:k , s(q) = 1 + 

J -L/2 



2M 

Q 



for the transition amplitude (k\s, q). In case when k = nQ/2 with n integer, hybridization 
of the neighboring bands should be taken into account. This modifies the above relations 
by the following way 

CnQ/2 = 9(n) [a n C n>qn + (3 n C n+hq J + 6(-n) [/3* n C-„, gn - a*_ n C '- n +i,q n ] , 

<ln = Q ^ j ^ l - ^^j > C s,q s = &* s C s Q/2 + f3 s C- s Q/ 2 , C S +l,g s = (3*C s Q/ 2 ~ a s C^ sQ / 2 , 

where a, (5 are hybridization coefficients. Corresponding relations between wave functions 
follow immediately from these formulas. 

Alternatively, each mode with the quasi-wavenumber q in the energy band s (reduced 
BZ description) can be described by the wave number k = q + mQ (extended BZ descrip- 
tion), where m = (— l) s sign(g) [s/2\. Within this scheme, Bloch amplitude and dispersion 
law uj{k) have the following simple form Uk(x) = e lkx , u(k) = v\k\. 
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To write down any of these formulas for a specific array, one should add the array 
index j to the wave function ip, Bloch amplitude u, coordinate x, quasimomentum q, and 
to the periods a and Q of the super-chain in real and reciprocal space. 
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Appendix B 

Spectrum and Correlators of Square 
QCB 

B.l Square QCB Spectrum 

Here we obtain analytical expressions for dispersion laws and wave functions of QCB. For 
quasimomenta far from the BZ boundaries, the energy spectrum of the first band can be 
calculated explicitly. Assuming that oj 2 <C u 2 (qj), s = 2, 3, 4, . . ., we omit u 2 in all terms 
in the r.h.s. of Eq. (|3.15J) except the first one, s = 1. As a result the secular equation 
(j3.15j) reads 

^(R, C?feV i 
where oo 2 (q) = w\(q). The solutions of this equation have the form: 



2 gi(q) +^ 2 2 (q) , 

<^/lq - g ± 



^?(q) - ^ 2 2 (q) 



1/2 

(B.l) 



Here q = 2 Ci (?i)Ci feV 2 ((?i V 2 (<?2), ^ = +, — is the branch number, Uj(q) is determined 



as 



^(q)=^te) 1 _ £(Fo _ y2fa))(Fo _ y2M) , ^fe) = T C,fe) (B.2) 

for j = 1. Expression for u5f (q) can be obtained by permutation 1 <-> 2. The expression 
inside the parentheses on the r.h.s. of Eq. ()B.2|) describes the contributions to F from 
higher bands. Therefore u;|(q) is the j-th array frequency renormalized by the interaction 
with higher bands. In principle, contribution of higher bands may turn the interaction 
to be strong. However, for the specific case of carbon nanotubes, one stays far from the 
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critical value e c (see estimates at the end of subsection 13. 4|) . Therefore the interaction 
with higher bands is weak in almost all the BZ except its boundaries. 

The resonance line equation is u(qi) = oj(q2)- Out of this line the branch number is in 
fact the array number and the renormalized frequencies are frequencies of a boson prop- 
agating along one of the arrays slightly modified by interactions with the complementary 
array. In case when u(qi) > u;(g 2 ), one obtains 

<i q « uj 2 ( gi ) (1 - eF 0lf 2 ( qi )) , u 2 _ M « u 2 (q 2 ) (1 - eF oV 2 (q 2 )) . (B.3) 

In the opposite case one should replace indices + <-> — . 

Consider the frequency correction in the latter equation in more details. The correction 
term can be approximately estimated as u 2 (qj)S(qj) with 

S(q) = eF oV 2 {q) = e^C!(q) J #C 2 (0- (B.4) 

Due to the short-range character of the interaction, the matrix elements Ci(l) ~ 1 vary 
slowly with the quasimomentum q < Q. Therefore, the r.h.s. in Eq. (|B.4|) can be roughly 
estimated as S(q) ~ eR^/a = O.lRo/a <C 1. One should also remember that the energy 
spectrum of nanotube remains one-dimensional only for frequencies smaller than some 
uj m . Therefore, an external cutoff arises at s ~ ak m where k m ~ ou m /v. As a results one 
gets an estimate S(q) ~ ek m R\ja. Hence, one could hope to gain additional power of 
the small interaction radius. However, for nanotubes, k m is of the order of 1/Rq (see 
Refs. |4TH I18j) and both estimates coincide. For quasimomenta close to the BZ center, 
the coefficient S(q) can be calculated exactly. For exponential form of £(£) oc exp(— 1£|), 
one obtains, S(0) = 0.14R /a. Thus, the correction term in Eq. (jB.3|) is really small. 

The eigenstates of the system are described by renormalized field operators. Within 
the first band they have the form 

oo 

#llq = "lq (Wq011q ~ ^q#21q) ~ (01sq M q#2sq + fasqVqOlsq) , (B.5) 

s=2 
oo 

#21q = « 2 q (fq^llq + Mq#21q) ~ (01 S qVq#2sq + 02sqMq#l S q) • (B.6) 

s=2 

Here the coefficients w q and v q describe mixing between the modes with different array 



indices, within the first band, 



'A* + < + A q 



1/2 




and 



^q = ~ , , ^q = ~ , , (B.7) 



A q = (^ 2 (g 2 ) - uj 2 (qi))/2, 0i q = v / ^Ci(9i)Ci(?2)^i(gi)wi(g 2 ), 
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£ = 4>a/R . The parameters 0j sq , s = 2,3,..., in Eqs. f)B.5|) . ()B.6|) correspond to 
inter-band mixing, 

and the coefficients «i q take into account corrections from the higher bands 



s=2 



Expressions for (/> 2sq and a 2q can be obtained by permutation 1 «-> 2. 

Equations (|B.1|) . ()B.5|) and ()B.6|) solve the problem of QCB energy spectrum within 
the first band and away from the BZ boundaries. However, due to smallness of the 
interaction, the general expressions for the eigenstates of QCB for arbitrary energy band 
can be obtained. For quasimomenta far from the high symmetry lines kj = ±Q/2, the 
energy spectrum can be calculated explicitly. Let us consider the s-th energy band and 
assume that 

\u 2 - (Qj ) I < tf^l [qj), s'^s. 

Moreover, we consider the wave vector lying away from the resonant lines k\ = ±/c 2 + mQ 
with m integer. Then, in zeroth approximation, the QCB eigenstates coincide with "empty 
QCB" eigenstates |1, si, gi) and |2,s 2 ,<Z2), i-e., belong to a given array and energy band. 
(Here we represent the wave vector k as k = q + m, where the quasi-momentum q belongs 
to the first BZ, whereas m = mi + m 2 is reciprocal super-lattice vector, m,, = rrijQej, 
j = 1,2, rrij is integer. The band numbers Sj are Sj = 1 + [2|fcj|/Q].) The inter-array 
interaction mixes the modes propagating along both arrays. The eigenstates of QCB are 
described by renormalized field operators 

ft n, ft i ■ ^ Cs'feVs'feKsfaVsfah , R , 

01,s,q - ®l,s,q?l,s,<i + <P > . 27 \ 2? \ y 2,s',q, (o-8j 



where 



oo 
s'=l 



Ca'(g2)^(g2)Cs(gl)^ B (gl 

o#(g 2 ) -o;2(g x ) 



The corresponding eigenfrequency is a slightly modified frequency of bosons propagating 
along the first arrays, 



Expressions for 2 ,s,q> a 2,s,q5 an d ^2,s,q can be obtained by permutation 1 <-> 2. 
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As was mentioned in subsection 13.51 description of QCB plasmons in terms of array 
plasmons fails at the points lying at two specific groups of lines in reciprocal space. 
Nevertheless, QCB plasmons even at these lines can be considered as finite combinations 
of array plasmons generated by a more complicated way. Below we consider various types 
of the interference of array plasmons in QCB. 

The first group of specific lines is formed by the high symmetry lines kj = PjQ/2 with 
Pj an integer (the corresponding quasi-momentum has a component qj = for pj even 
or Q/2 for pj odd, the reciprocal lattice vector is m,, = [|pj|/2]). In an "empty lattice" 
approximation, these lines are degeneracy lines which separate Sj-th and (sj + l)-th bands 
of the j-th array (sj = \pj\). Inter-array interaction mixes the array states, lifts the (inter- 
band) degeneracy and splits the corresponding frequencies. As a result, the QCB plasmons 
related to a point k lying at a high symmetry line (e.g. all C, D points in Fig. 16.2)1 . are 
built from the array plasmons associated not only with the points k 12 but also with the 
symmetric (with respect to coordinate axes of a reciprocal space) points ky = (—ki,q 2 ) 
and k 2 = {q%, —k 2 )- (Indeed, if the condition kj = PjQ/2 is fulfilled, — kj = kj — PjQ, and 
then both these wave-numbers correspond to the same quasi-wave-number qj in the first 
BZ.) 

Consider for definiteness the case j — 1 and assume first that the ratio 2k 2 /Q is non- 
integer. In this case, the point k = q + m generates three QCB plasmons. The first of 
them is plasmon 02,s 2 ,q QB.8j) propagating along the second array with quasimomentum 
q 2 and frequency cu 2jS2iq ()B.9|) . The two others 



01,ai,5,q = — ^ (01,si,q + 01,ai+l,q) , (B.10) 
<W,q " ^(^ 1 ,q-^ 1 +l,qJ+V20^ ^(fc) " 

are even or odd superpositions of the 1-st array states from the two zones s\ and si + 1 
with eigenfrequencies 



~,2 



U) 



X> (g 2 )^ s 2 . (q2)(i {qiW Sl (qi 



l,si,S,q 



w M llU , q = w sl (9i)-^ L — ; ,2 (n \ ,,2 (n ^ • ( Rn ) 

t~. u s ,{q 2 ) - to Si {q 1 ) 



The case j = 2 is described similarly after change 1 «-> 2. 

The second set is formed by the resonant lines defined by the equation k% + rk 2 = nQ, 
where r = ±1 and n are two integer parameters determining the line. The corresponding 
quasi-momentum satisfy the equality q\ + rq 2 = 0, then it lies on the diagonal of the first 
BZ (see Fig. 13. 4jl . In an "empty lattice" approximation, these lines are also degeneracy 
lines. However this degeneracy has more complicated nature. Here there is a dual point 
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k = (— rk 2 , —rki) which satisfies the equalities — rk 2 = k\ — nQ and — rk\ = —k 2 + rnQ, 
and therefore this point corresponds to the same quasimomentum q as the point k. As a 
result, the array plasmons associated with the points k and k are involved in the resonance. 
Inter-array interaction mixes degenerate modes, lifts the degeneracy, and splits degenerate 
frequencies. As a result, the QCB plasmons related to a point k at the resonant line (e.g. 
all B, F, E points in Fig. 16. 2|) . are built from the array plasmons associated not only with 
the points ki 2 but also with the symmetric (with respect to one of the bisector lines of the 
coordinate system of a reciprocal space) points ki = (—rki, —rq 2 ) and k 2 = (—rqi, —rk 2 ). 
Moreover, if the point k lies at an intersection of two resonant lines (the point W in Fig. 
I3.4J) . the QCB plasmons include also array states associated with the kj points. 

In the general case, the point k = q + m lies only on one of the resonant lines q± = ±q 2 
away from the high symmetry lines. Here the point k generates two pairs of even and odd 
superpositions of the first and the second array states (two doublets). The first of them, 
9g/u, ai w lies in the energy band s x = 1 + [2|&i|/Q], 

a ft , C 8 '(g2)Mg2)Csi(giHi(gi) fl / r10 n 

0<?/it,Si,q = — ^ (^l,si,q ± ^2,si,q) , j = 1, 2. 

The corresponding eigenfrequencies form two doublets 



The eigenfunctions and eigenfrequencies of another doublet are obtained similarly by 
replacing the indices 1 «-> 2. 

The crossing points of the resonant lines can also be divided into two groups. The first 
group is formed by the crossing points of only two high symmetry lines (points X\ and X 2 
in Fig. 13. 3J) . Here, each of the generated vectors kx and k 2 corresponds to a pair of array 
eigenstates belonging to two adjacent ID bands. Thus an inter- band mixing is significant 
in both arrays. Consider for definiteness the point X\. In this case the wave vector k (kj) 
of the extended BZ can be represented as k = q + m (kj = q + m-,), with q = (Q/2, 0). 
As a result, in "empty lattice" approximation, we have two sets of degeneracy lines. The 
first one of them corresponds to odd and next even bands of the first array (q\ = Q/2), 
whereas the second one describes degeneracy of even and next odd bands of the second 
array (q 2 = 0). The corresponding QCB plasmons are even or odd combinations of the 
j-th array plasmons (j = 1, 2). For j = 1, we have, 

S -ft a -ft i o-tf C2n(g2V2n(g2)Ci(glVsi(gl) a (T\ ^ A\ 

fl,si,p,q — fl,si,g,q) fl,si,u,q — fl,si,u,q T £<p > ~ 7 r Yl \ t'2,2n,u,q* l - tJ -- L4 J 

~{ W 2nW2j - U (qi) 



86 



where 



The corresponding eigenfrequencies are 

" (Li^PLi^Xl (?i) 



cJ 2 



i,si>s,q 



w lj , li14iq = a; ai (gi) -40 ^ ; ,2 /„ n , ,2 (n ^ • ( R15 ) 



n=l 



The case j = 2 is described similarly after change 1 •*->• 2 and 2n — > 2n + 1 in Eqs. (jB.14|) 
and (iRToTl . 

The second group consists of the crosses of two resonant lines (points T and W in Fig. 
13.4)1 . These points are also the crosses of two high symmetry lines (this is always true). 
Here the QCB plasmons, generated by the point k = q + m with q = (q, q) = (0, 0) (r 
point) or {Qq/2, Qq/2) (W point), form two quartets. The first quartet is really generated 
by the point ki. It consists of four symmetrized combinations of single-array states 

a _n n _n , \- C' Vs' feRsx (?1 V*i fah , R , ^ 

where v — g,u; the band numbers Si and s' are even for q = and odd for q = Q/2\ 

®s,g/u,g,q = ~ (^j,s,q ^ $l,s+l,qj j @s,g/u,u,q = ~ (@j,s,q ^ ^j,a+l,qj • 

J'=l,2 J=l,2 

QCB eigenstates of the second quartet are generated by the point k 2 . They are obtained 
from the equations (jB.16|) by replacing 1 <-> 2. Even array eigenstates are degenerate with 
the frequencies 

^, M / M = Sfe) 2 : ( B - 17 ) 

while the odd array eigenstates are split 

/ uo n 2/ \ 4 ,2 \- C(g2)^ 8 2 /(g2)C 2 1 (glX(gl) 



Equations ()B.10|) - ()B.18|) exhaust all cases of QCB plasmons, generated by the point 
k 2 lying on a specific line, via the interference of array plasmons. 

B.2 AC Conductivity 

For interacting wires, where (j S {qj) ^ 0, the correlator (|3.19)l may be easily calculated 
after diagonalization of the Hamiltonian ()3.7|) by means of the transformations considered 
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in Appendix lB.il For the first band, for example, these are the transformations (jB.5|) and 
(IB. 61). As a result, one has: 



■W*),4lq(0) 

■M*),4i q (o) 



-2ivg (uq^ +> i q sin(^ +i i q t) + v^uj- M sxn(uj- M t)) , 
^ivgu^Vq (u;_ lq sin(u;_ i q *) - w +)lq sin(u; +)lq t)) . 



where ii q and v q are defined in Eqs. ()B.7|) . Then, for the optical absorption a' one obtains 



Trvg 



U 2 5 ((jj - U+,lq) + V^5 (u - &-,lq) 



TTVgUqVq 



8{UJ- £-,lq) -5(0J- UJ +M ) 



(B.19) 
(B.20) 



1 and v q ~ 



if q for 



For quasimomentum q away from the resonant coupling line, 
A q > (fq ~ 1 and ui ~ </>f for A q < 0). Then the longitudinal optical absorption 
()B.19J) (i.e. the absorption within a given set of wires) has its main peak at the frequency 
LJ +M ~ v \Ql\ f° r > ( or ^-.lq ~ v \Qi\ f° r ^q < 0)) corresponding to the first 
band of the pertinent array, and an additional weak peak at the frequency cj_ )lq ~ v\q 2 \, 
corresponding to the first band of a complementary array. It contains also a set of weak 
peaks at frequencies u>2 }S q ~ [s/2]vQ (s = 2,3,...) corresponding to the contribution 
from higher bands of the complementary array (in Eq. (|B.19|) these peaks are omitted). 
At the same time, a second observable becomes relevant, namely, the transverse optical 
absorption (|B.20|) . It is proportional to the (small) interaction strength and has two peaks 
at frequencies cJ+,i q and W-,i q in the first bands of both sets of wires. 

If the quasimomentum q belongs to the resonant coupling line A q = 0, then ui = 
fq = 1/2. In this case the longitudinal optical absorption ()B.19|) has a split double peak 
at frequencies w+,i q and w_,i q , instead of a single main peak. The transverse optical 
absorption ()B.20)1 . similarly to the non- resonant case ()B.20)1 . has a split double peak at 
frequencies ct>+,i q and W-,i q , but its amplitude is now of the order of unity. For |q| — > 
Eq. flB.19|) reduces to that for an array of noninteracting wires (|3.2U|) . and the transverse 
optical conductivity (|B.20|) vanishes. 

The imaginary part of the ac conductivity cr"-,(q, a>) is calculated within the same 
approach. Its longitudinal component equals 



2vg 



OJ 



0J\ 



,lq ~ u ^-,iq 

Beside the standard pole at zero frequency, the imaginary part has poles at the resonance 
frequencies ^+,i q , ^-,i q , and an additional series of high band satellites (omitted here). 
For quasimomenta far from the resonant lines, only the first pole is well pronounced 
while amplitude of the second one as well as amplitudes of all other satellites is small. 



+ 



0J Z 
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At the resonant lines, the amplitudes of both poles mentioned above are equal. The 
corresponding expression for <7 22 (q, a;) can be obtained by replacement 1^2. 

The transverse component of the imaginary part of the ac conductivity has the form: 



'12 



2vg 



LO 



LO 



,lq 



LO 



+,iq 



CO' 



00 



u ~ ^+,iq- 

It always contains two poles and vanishes for noninteracting wires. For quasimomenta 
far from the resonance lines the transverse component is small while at these lines its 
amplitude is of the order of unity. 



col 
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Appendix C 
Tilted QCB 



C.l Geometry, Notions and Hamiltonian 

A tilted QCB is a 2D grid, formed by two periodically crossed arrays of ID quantum 
wires or carbon nanotubes. Like in a square QCB, arrays are labelled by indices j — 1,2 
and wires within the first (second) array are labelled by an integer index n-i [n\)- The 
arrays are oriented along the unit vectors ei^ with an angle ip between them. The periods 
of a crossbars along these directions are a\ and and the corresponding basic vectors 
3X6 Si j == CljGj. In experimentally realizable setups, QCB is a cross-structure of suspended 
single-wall carbon nanotubes lying in two parallel planes separated by an inter-plane 
distance d. Nevertheless, some generic properties of QCB may be described under the 
assumption that QCB is a genuine 2D system. We choose coordinate system so that the 
axes Xj and corresponding basic unit vectors are oriented along the j-th array. The basic 
vectors of the reciprocal superlattice for a square QCB are Qi,2gi,2 ; Qj — 2n/aj so that an 
arbitrary reciprocal superlattice vector m is a sum m = mi + 1112, where nx, = mjQjgj, 
(rrij integer). Here gi^ are the vectors of the reciprocal superlattice satisfying the standard 
orthogonality relations (e^ • gj) = <5y. The first BZ is | 2 1 < Qi.,2/2. 

A single wire of j-th array is characterized by its radius rj, length Lj, and LL in- 
teraction parameter gj. The minimal nanotube radius is 0.35 nm[Hl], maximal nanotube 
length is L = 1 mm, and the LL parameter is estimated as 0.3 - In typical experimental 
setup the characteristic lengths mentioned above have the following values 

dw2 nm, Lj w 0.1 mm, 

so that the inequalities 

r 1)2 < d < a 12 < Li j2 
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are satisfied. 

The QCB Hamiltonian 

H = H 1 + H 2 + H 12 . (C.l) 
consists of three terms. The first two of them describes LL in the first and second arrays 

Li/2 



H2 -In/2 



L 2 /2 

H 2 = J dx ^92^l (niQ!, x 2 ) + — (d X2 9 2 (n 1 a 1 ,x 2 )) 2 



ni -L2/2 



x 



The inter-array interaction is described by the last term in Eq. (|C.lJ) 

2e^ ^ \ r ( Xi - niaA ( n 2 a 2 - x 2 \ 

Hl2 = T 22 J dx ^ {—^-~ J c {~^2~~ ) 

ni,na 

xd Xl 9 1 (x 1 ,n 2 a2)d X2 9 2 (nia 1 , x 2 ). (C.2) 

It results from a short-range contact capacitive coupling in the crosses of the bars. 
The dimensionless envelope function (introduced phenomenologically) C(£j) describes re- 
distribution of a charge in a tube j induced by the interaction with tube i. This function 
is of order unity for |£| ~ 1 and vanishes outside this region so that the dimensionless 
integral 



J C(Oe lkr ^ ~ 1 



is of order unity for all \k\ smaller than a certain ultraviolet cutoff. 

The QCB Hamiltonian (jC.lj) is a quadratic form in terms of the field operators, so it 
can be diagonalized exactly. Such a procedure is rather cumbersome. However, due to 
the separability of the interaction (jC.2|) the spectrum can be described analytically. 



C.2 Spectrum 

Now we consider the spectrum of a generic double QCB. The resonance condition (|3.18j) is 
fulfilled not at the BZ diagonal but at the resonant polygonal line. Its part ODE, lying in 
the first quarter of the BZ, is displayed in Fig. Kill (all figures of this Section correspond 
to the specific values v 2 Q 2 = 1, V1Q1 = 1.4). This results in qualitative modifications of 
the spectrum that are related first of all to the appearance of two points D and E of the 
three-fold degeneracy for a titled QCB (Fig. IC.1J) instead of a single point W of four-fold 
degeneracy for a square QCB (Fig. 13.4)1 . 
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Figure C.l: BZ of a tilted QCB. 



We start with the resonant line ODE (Fig. 13. 3J) . The dispersion curves at its OD 
part and the symmetry properties of the corresponding eigenstates are similar to those at 
the OC resonant line for the square QCB (Fig. 13 .4jl . The only difference is that instead 
of the four-fold degeneracy at the BZ corner C of the square QCB, there is a three-fold 
degeneracy at the point D lying at the BZ boundary. A completely new situation takes 
place at the DE line, where two other modes (1, 1) and (2,2), corresponding to different 
arrays and different bands, are degenerate. The interaction lifts this degeneracy and the 
two middle lines in Fig. 1(121 describe even (g) and odd (it) combinations of these modes. 
The even mode corresponds to the lowest frequency and the odd mode corresponds to the 
higher one. At the point E one meets another type of a three-fold degeneracy described 
in more detail in the next paragraph. 

Dispersion curves corresponding to quasi momenta lying at the BZ boundary q\ = 
Qi/2, < q 2 < Q 2 /2 (FC line in Fig. £HJ and q 2 = Q 2 /2, < q l < Q l /2 (CF' line 
in Fig. IC.lj) . are displayed in Fig. IC.3I The lowest and highest curves in the FE part 
of the latter figure, describe two waves propagating along the second array. They are 
nearly linear, and deviations from linearity are observed only near the point E where 
the interaction has a resonant character. Two modes propagating along the first array, 
in zero approximation, are degenerate with an unperturbed frequency uj = 0.7. The 
interaction lifts the degeneracy. The lowest of the two middle curves corresponds to (1, u) 
boson, and the upper of one describes (l,g) boson. Note that (l,g) boson conserves its 
unperturbed frequency to = 0.7. The latter fact is related to the symmetry Cj(0 = Cj( — 
of the separable interaction ()3.3|) . At the point E, the two modes propagating along 
the first array and the mode propagating along the second array in the second band are 
degenerate. The interaction lifts the degeneracy, and, as a result, the (l,u) and (2,2) 
waves are strongly mixed and the eigen-modes are their even (highest frequency) and odd 
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(lowest frequency) combinations, and the (l,g) mode (middle level). 
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Figure C.2: The energy spectrum of a tilted QCB (solid lines) and noninteracting arrays 
(dashed lines) for quasimomenta on the resonant line of the BZ (line ODE in Fig. IC.1|) . 
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Figure C.3: The energy spectrum of a tilted QCB (solid lines) and noninteracting arrays 
(dashed lines) for quasimomenta at the BZ boundary (line FCF' in the Fig. IC.lj) . 



There are two separate degeneracies within each array at the corner C of a titled QCB 
BZ. Both of them are related to inter-band mixing conserving array index. The spectral 
behavior along the CF' boundary of the BZ is similar to that considered above but in 
the vicinity of the point D of three-fold degeneracy. Here, two modes propagating along 
the second array in the separable potential approximation (|3.1()j) remain degenerate. This 
degeneracy is lifted only if deviation from separability is accounted for. 

The diagonal OC of a tilted QCB BZ represents a new type of generic line, that crosses 
a resonant line (Fig. IC.4J) . Here the spectrum mostly conserves its initial systematics, 
i.e. belongs to a given array, and mostly depends on a given quasimomentum component. 
However, at the crossing point B, the modes (1,1) and (2,2), corresponding to both 
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Figure C.4: The energy spectrum of a titled QCB (solid lines) and noninteracting arrays 
(dashed lines) for quasimomenta on the BZ diagonal (line OC in Fig. IC.lj) . 

f>-^ q2/Q 2 




Figure C.5: Lines of equal frequency for a tilted QCB (solid lines) and noninteracting 
arrays (dashed lines). Lines 1,2,3 correspond to frequencies uj\ = 0.1, u 2 = 0.25, uj^ = 
0.45. 

different arrays and bands, become degenerate (two middle dashed lines in Fig. IC.4J) . 
Interaction between the wires lifts the degeneracy. The eigenstates of QCB have a definite 
parity with respect to transposition of these two modes. The lowest and upper of the two 
middle lines correspond to even (g) and odd (u) mode, respectively. 

Like in a square QCB, bosons with quasimomenta close to the resonant lines are 
strongly mixed bare ID bosons. These excitations are essentially two-dimensional, and 
therefore lines of equal energy in the vicinity of the resonant lines are modified by the 2D 
interaction (see Figs. IC.5l and lC.6J) . Deviations from ID behavior occur only in this small 
part of the BZ. For uj < 0.5^2^2 the lines of equal energy within BZ consist of closed line 
around the BZ center and four open lines (within the extended bands scheme these lines 
are certainly closed) around the BZ corners (lines 1, 2, 3 in Fig. IC.5|) . At the line OD in 
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Figure C.6: Lines of equal frequency for a tilted QCB (solid lines) and noninteracting 
arrays (dashed lines). Lines 4,5 correspond to frequencies uj^ = 0.55, CJ5 = 0.65. 

BZ, the modes of QCB are strongly coupled bare bosons propagating along both arrays 
in the first band. 

For 0.5t>2Q2 < uj < 0.5i>iQi (lines 4, 5 in Fig. IC.6J) the topology of lines of equal 
energy is modified. In this case, lines of equal energy within the BZ consist of four open 
lines. The splitting of lines at the direction DE corresponds to strong coupling of modes 
propagating along the first array in the first band with those propagating along the second 
array in the second band. 
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Appendix D 
Triple QCB 



D.l Notions and Hamiltonian 



Triple quantum bars is a 2D periodic grid with m — 3, formed by three periodically 
crossed arrays j = 1, 2, 3 of ID quantum wires. These arrays are placed on three planes 
parallel to the XY plane and separated by an inter-plane distances d. The upper and the 
lower arrays correspond to j = 1,2, while the middle array has number j = 3. All wires 
in all arrays are identical. They have the same length L, Fermi velocity v and Luttinger 
parameter g. The arrays are oriented along the 2D unit vectors 

--m »=«»>' «=«- 

The periods of QCB along these directions are equal, ctj = a, so we deal with a regu- 
lar triangular lattice. In what follows we choose a 12 = ae 12 as the basic vectors of a 
superlattice (see Fig. ID. If) . 

The wires within the j-th array are enumerated with the integers rij. Define 2D coor- 
dinates along the n 3 -th wire rj as rj = xjej + njae 3 for upper and lower arrays (j = 1,2) 
and r 3 = X3e 3 + n 3 aei for the middle array. Here Xj are ID continuous coordinates along 
the wire. The system of three non-interacting arrays is described by the Hamiltonian 

Hq = Hi + H 2 + H 3 , 



where 



Ho 



hv 



hv 



hv 



j dxi 



E 

712 

E 



dx- t 



dx-i 



gnl(xiei + niae 3 ) + - {d^OAxxei + niae 3 ))' 
9 



gTcl(x 2 e 2 + n 2 ae 3 ) + - (d X2 6 2 (x 2 e 2 + n 2 ae 3 )Y 

1 

5f7r 3 (a; 3 e 3 + n 3 aei) + - (d X3 9 3 (x 3 e 3 + n 3 ae!))^ 



(D.l) 
(D.2) 
(D.3) 
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and 7r j and d x .9j are canonically conjugate fields describing LL within the j-th array. 

3\/l 




Figure D.l: Triple QCB. 




Figure D.2: Elementary cell BIJL of the reciprocal lattice and the BZ hexagon of the 
triple QCB. 

Interaction between the excitations in different wires of adjacent arrays j,j' is con- 
centrated near the crossing points with coordinates TijSLj + rij'&ji. It is actually Coulomb 
interaction screened on a distance r along each wire which is described by Hamiltonian 

Hint = H13 + ^23; 



where 



H 



J"3 



Vq J dxjdx 3 § 



Xj — n 3 a x 3 — n,ja 
— e,' — e 3 I x 



r 



r 



x 



d Xj 9 j (x j e j + n j ae 3 )d X3 9 3 (n 3 ae j + x 3 e 3 ). 



(D.4) 



Here the effective coupling strength Vq is defined by Eq. (13.5)1 . the dimensionless interaction 
$ is separable 

H&j + £ 3 e 3 ) = C(&)C(&), j = 1, 2, (D.5) 
and is a dimensionless charge fluctuation in the j-th wire (see Eq. (jrS.rij) ). 
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Such interaction imposes a super-periodicity on the energy spectrum of initially one 
dimensional quantum wires, and the eigenstates of this superlattice are characterized 
by a 2D quasimomentum q = gigi + g2g2 = (<?i,<?2)- Here gi 5 2 are the unit vectors 
of the reciprocal superlattice satisfying the standard orthogonality relations (e« ■ gj) = 
5ij, j = 1,2. The corresponding basic vectors of the reciprocal superlattice have the form 
Q(m,igi + m2g2), where Q = 2n/a and m 12 are integers. In Fig. ID.2l the elementary cell 
BIJL of the reciprocal lattice is displayed together with the hexagon of the Wigner-Seitz 
cell that we choose as the BZ of the triple QCB. 

To study the energy spectrum and the eigenstates of the total Hamiltonian 

H = H + H int , (D.6) 
we define the Fourier components of the field operators 

0i(ziei+mae 3 ) = -±= J^^e^+^Us^), (D.7) 



'NL 

s,q 



e 2 ( X2 e 2 + n 2 ae 3 ) = -jL= ^ 9 2s(l e^ + ' 1 ^u Stq2 (x 2 ), (D.8) 

V iV Lj 

s,q 

^(xses + nsaex) = £ 9 3s ^ x ^ n ^u sm {x 3 ). (D.9) 



'NL 

s,q 



Here 



Q = <?igi + <?2g2, qa = Qt-qi, 

and N = L/a is the dimensionless length of a wire. In the q representation, the Hamil- 
tonians Hj (Eqs. (ID. 1}) - (ID. 3^ and H j3 (Eq. (ID .4)) ) can be written as 

H o = — Z^ 7r Sq 7r ^ + ^ZX^iq^' J = 1,2,3, 

s,q s,q 

#j3 = X] Ml^fis'ilj^siqs^s'iqj) {QtsqQjs'q + h.c.) , j = 1, 2, 

s,s',q 

where 

w.(g)=«([|]g + (-l)- 1 M), Q = ^, 

Thus the total Hamiltonian (|D.6j) describes a system of coupled harmonic oscillators, and 
can be diagonalized exactly like in the case of double QCB. 



D.2 Spectrum 

Separability of the interaction (jD.5|) allows one to derive analytical equations for the 
spectrum of the total Hamiltonian ()D.6)) . Here we describe the behavior of the spectrum 
and the corresponding states along some specific lines of the reciprocal space. 
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To diagonalize the Hamiltonian (jD . 6|) . we write down equations of motion 

[uj 2 s {qj) - uj 2 ] e jsC[ + \/e(j) s (q j )u s (q j )—^2 <M<?3V s '(<?3)#3s' q = 0, (D.10) 

a s' 



a 

j,s' 



Here j = 1,2, and e is defined by Eq. (|3.17|) . The solutions of the set of equations (|D.10j) 
- (jD.ll|) have the form: 



Substituting this equation into Eqs. ljD.lOj) and (jD.llj) . we have three equations for the 
constants Ay. 

A x + A 3 ^F qz {uo 2 ) = 0, A 2 + A 3 ^F q3 (u 2 ) =0, A 3 + ^ A,-V^ 9 > 2 ) = 0, 

i=i,2 

where 

Dispersion relations can be obtained from the solvability condition for this set of equations 

eF q3 (uj 2 ) {F qi (u 2 ) + F q2 (u 2 )) = 1. 

The function F qa {oj 2 ) has a set of poles at uj 2 = uj 2 (q), s = 1, 2, 3, . . . . For cj 2 < ou 2 (q), i.e. 
within the interval [0,c<j 2 (g)], F qs (uj 2 ) is positive increasing function. Its minimal value F 
on the interval is reached at uj 2 = and does not depend on the quasi-momentum q 



F «(°) = 7 £#(*) = / ^(0 = ^ 



If the parameter e = rj 2 is smaller than the critical value 

1 



then all the solutions uj 2 of the characteristic equation are positive. When e increases, the 
lowest QCB mode softens and its square frequency vanishes in the whole BZ at e — e c . 
For the exponential interaction model = exp(— 1£|), one obtains e c ~ 1. 

The high symmetry of the triple QCB leads to a number of lines where inter-array 
or inter-band resonant interaction occurs: all lines in Fig. ID. 21 possess some resonant 
properties. These lines may be classified as follows: 

On the Bragg lines where one of the three array wave-numbers qj is a multiple integer 
of Q/2, there is a strong intra-band mixing of modes of the j-th array. In Fig. ID.2[ these 
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lines are the boundaries of the elementary cell of the reciprocal lattice I J LB, axes qi and 
q 2 , lines OB and EH. In particular, along the lines OA (q 2 = 0) and OB (q 3 = 0) two 
modes corresponding to the second and third bands and to the second (OA) or third (OB) 
array are mixed. Along the line AB (qi = Q/2) the same mixing occurs between (1, 1) 
and (1, 2) modes. Moreover, the resonant mixing of different arrays within the same band 
occurs along the medians OA, OB, etc. There are two types of such resonance. The first 
one (e.g., OA line) is the resonance between neighboring arrays (q\ = —(73) and therefore 
it is of the main order with respect to interaction. The second one (e.g., OB line) is the 
resonance between remote arrays (qi = q 2 ) and it is one order smaller. 

The second family consists of resonant lines formed by the BZ hexagon boundaries 
and diagonals. Thus, the diagonal OC realizes a first order resonance between the first 
and the third arrays q\ = (73, and the BZ boundaries HD and AN correspond to the same 
resonance up to an umklapp process (qi = q 3 — Q and q\ = q 3 + Q respectively). Along 
the diagonal OD and the BZ boundary NC a second order resonance takes place with 
resonance conditions q 2 = —qi and q 2 = —qi + Q respectively. 

In the reciprocal space of the triple QCB there are four different types of crossing 
points. Two of them include the bases of BZ medians (e.g., points A, B, E and so on). 
Here one deals with the four-fold degeneracy of the modes corresponding to the first order 
resonance between the neighboring arrays (e.g., point A, ujx jS = CJ3 )S >, s,s' = 1,2), or to 
the second order resonance between remote arrays (like point B, lo^ s = uj 2)S i, s, s' — 1,2). 
One more family consists of crossing points of the BZ diagonals and the lines connecting 
the bases of its medians (points M, F, G and so on). Here one deals with three types of 
two-fold degeneracy simultaneously. For example, at the point M two separate pairs of 
modes corresponding to neighboring arrays (2, 1), (3, 1), and (2, 2), (3, 2), are degenerate, 
as well as two modes corresponding to the first array, (1, 1), (1, 2). Finally the BZ hexagon 
vertices form the most interesting group of points where the three-fold degeneracy between 
modes corresponding to all three arrays takes place. The typical example of such a point 
is the vertex C where the resonance condition qi = —q 2 + Q = q 3 = Q/3 is satisfied. 

Almost all these peculiarities of the triple QCB spectrum can be illustrated in Fig. 
ID. 31 where the dispersion curves along the closed line OABO are displayed. We emphasize 
once more that in the infrared limit u, q — > triple QCB like double QCB preserves the 
characteristic LL properties of the initial arrays. 
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O A M B O 

Figure D.3: Dispersion curves at the OAMBO polygon of BZ. 

D.3 Observables 

The structure of the energy spectrum analyzed above strongly influences optical and 
transport properties of the triple QCB. As in the case of double QCB (part 13. 6|) . one 
expects to observe four peaks of the optical absorption near the points A, B, E, H of the 
four-fold degeneracy. Then, specific features of space correlators like those considered in 
ID. 21 can be observed. But the most pronounced manifestation of a triangular symmetry 
of the triple QCB are its Rabi oscillations. 

Consider the vicinity of the point C of the BZ of three-fold degeneracy, mixing all 
three arrays, 

qi = <?3 = -<?2 + Q = ^u(<5/3) = cj 2 i(2<2/3) = u 31 (Q/3) = lu . 

Equations of motion at this point in the resonance approximation read 

9 3 + v/^o 2 (9, + 6 2 ) = 0, 
where 9j = Jq . General solution of this system looks as 
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where one of the eigenfrequencies coincides with ujq, while the other two are 



u± = oo y 1 ± V2(j) 2 , 
and 9o t ± are the corresponding amplitudes. 



dt 2 



9 3 + v^ 2 ^ = 0, 



dt< 
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Choosing initial conditions 0i(O) = i0„, 0\(O) = w o o , 0a(O) = 3 (O) = 0, 2 (O) 
6*3(0) = 0, we obtain for the field amplitudes at the coordinate origin 



0i(O,O;t) 
2 (0, 0;t) 
03(0,0; t) 



0o 
4 

0o 
4 

0o 



— sin(u; + rj H sm(u;_t) 

— sin(u; + rj H sm(u;_cj 



0n 



sin(a> t), 



2v/2 



w . / ,\ w . / ,x 

— sin(w + t) sin(w_tj 



In the limiting case e < 1 these formulas lead to the following time dependence of the 
field operators in the coordinate origin in real space 



0i(O,O;t) = O sia(u t) cos 2 



2 (O,O;t) = O cos(a;ot) sin 2 



UJ t 



Wot 



(D.12) 



3 (O,O;t) = O sin(o;ot) cos 



2e^ 



u t . 



The field operators of all three arrays demonstrate fast oscillations with the resonant 
frequency ujq modulated by a slow frequency. It is the same for the two remote arrays, 
and doubled for the intermediate array. These beatings are synchronized in a sense that 
zero intensity on the intermediate array always coincides with the same intensity on one 
of the remote arrays. At these moments all the energy is concentrated solely within one 
of the remote arrays. These peculiar Rabi oscillations are displayed in Fig. ID. 41 
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Appendix E 

Derivation of QCB-Light Interaction 
Hamiltonian 



Light scattering on QCB is described by equations (jfi.lj) - (|fi.7j) . In this Appendix we 
briefly explain the main steps which lead to such a description. 

1. Nanotube-light interaction, Eq. (J6.1|) . We start with a consideration of a 
single nanotube of the first array, interacting with an external electromagnetic field. The 
characteristic time of all nanotube energies, including Coulomb interaction, is of order of 
an inverse plasmon frequency. The scattering process occurs during much shorter time 
interval, which is of the order of an inverse photon frequency. Hence, Coulomb interaction 
is irrelevant for the scattering processes. This enables us to restrict ourselves by a kinetic 
part of the nanotube Hamiltonian. Within the k — p approximation, this part is given by 

m 

h k in = v F [ -^^(ii, r 0l 7)fik ■ ff Q a /^(ii, r , 7). (E.l) 

J 27T 

In the presence of a magnetic field, one should add eA/c to the electron momentum oper- 
ator Kk. An additional part of the Hamiltonian (jE.l|) exactly coincides with the nanotube- 
light interaction Hamiltonian ((SHJ- 

The Hamiltonian (|E.1|) is diagonalized by a two step canonical transformation. The 
first one is Fourier transformation 

M%ul) = 4= E c <*,k,me ikxi+im \ (E.2) 

k,m 

where the orbital moment m is restricted by the condition |m| < m = [7rr /ao] due to 
the finite number of honeycomb cells along the nanotube perimeter (czoa/3 is the lattice 
constant). The second rotation is defined as 

1 

71 



Cp,k,m — (pAkm + P e ^ m CBfcm) j P — ij 
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where 

, kro . , m 

cos^™ = — . sirii 



sj (kr ) 2 + m 2 a/ (kro) 

As a result, the Hamiltonian takes the form 



m 2 



H kin = ^p( k ^ m ) C lkm C pkm (E.3) 

p,k,m 



with eigenfrequencies 




u) p (k,m) = pvf\ k 2 H g-. (E.4) 



2. QCB-light interaction, Eq. (J6.3|) . Substituting into Eq. (|6.1|) the scalar product 
A ■ <j in the form 

(A-«x) = 

where A 1 * 1 = A\ ± zA 7 , we can write the nanotube-light interaction as 




c J 2tc V 

We are interested in an effective QCB-light interaction Hamiltonian obtained in second 
order of perturbation theory which describes transitions between initial states \i) and final 
ones |/). Initial states are one-photon states of the electromagnetic field and the electron 
ground state of the nanotube whereas final states consist of one photon and an electron 
above the Fermi level. It will be seen later that just these states form a one-plasmon array 
state. The energy of the incident photon E = hck is much higher than the excitation 
energies of the nanotube. Therefore, absorption of the incident photon and radiation 
of the scattered photon occur practically without retardation. All this results in the 
interaction Hamiltonian 

'-> 2 /7)„\2 



'int 



\v) (v\M^i, 7')4(zi, 7) \v) (v\A+(x u 0, 0)Mxu 7) + A <- b) , (E.5) 



which corresponds to the diagram shown in Fig. 16.1b (there is no photons in an interme- 
diate state). 

Consider now the matrix element (v\if) a (x' 1 , 7 / ) , 0J,(xx, j)\v) which enters this Hamilto- 
nian. Due to our choice of initial and final states, only diagonal elements with respect to 
both virtual states \v) and sublattice indices a survive. In the c a ^, m representation (|E.2jl . 
they have the form 

k,k' m,m' 
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Internal matrix elements on the r.h.s. of the latter equation are 

(here the symmetry property n (E_(k, m)) + n (E + (k, m)) = 1 is used). Therefore 

where 



mo 



2,71 ' ' 

tn=-rao 



Thus, the interaction (|E.5|) takes the form 
^ = ^(-) 2 / y g(7 ~ VM'^i, 0: 0M + (*i, 0, 0) rO^afo, 7)- 

C C J 7T q. 

The field dependent factor here is 

A"(a;i, 0, 0)A + (x!, 0, 0) = 0, 0) + A 2 ( Xl , 0, 0) sin( 7A - 7 ') sin( 7A - 7 ), (E.6) 

and (Ax,A,j) are cylindrical components of the vector potential A ()6.2|) . Taking into 
account the angular dependence of the field (|E.6[) . we can omit it in the electron operators. 
Indeed, according to Eq. (|E.4J) (see also Fig. 12. 4j) in an energy-momentum region where 
we work, the m — 1 spectral band with nonzero orbital moment is separated from that 
with m = by an energy of order of h~Vp/r , which is much higher than the QCB plasmon 
energy. Keeping only the zero moment field operators which form the electron density 
operator 

y^^lix^ajxi) = p(xi) = V2d Xl 9(xi), 

a 

and integrating over 7 , 7 ', we obtain an interaction Hamiltonian in form 

h in t = (t)V ^^OAi^i-O-O). (E.7) 

where 

Ai = A + (V2- l)Axe x . 

Straightforward generalization of this expression to the QCB case leads exactly to the 
Hamiltonian Ijfi.Hj) . 

3. Polarization matrix, Eq. (|6.7|) . To study the scattering process, we should modify 
the last expression for the interaction Hamiltonian. To proceed further, we define Fourier 
transforms 9j,Qj of the bosonic fields 

6 1 (x 1 ,n 2 a) = -L=^^_ ki e- ifc ^-^ 2a , 
V iVL , 

6 2 (n ia ,x 2 ) = -4E^-^" %m " ite . (E.8) 
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Here N = L/a is the number of QCB cells in both directions. The electromagnetic 
field also can be expanded in a sum of harmonics with a wave vector k and polarization 

A = ||, -L, 

A(r) = 5>KA^KAe lKr . (E.9) 

KA 

The polarization vectors 

iK 2 e 1 iKxe 2 

nK " ~ ^/W+M'7W+M } 

Ky/K* + K\ Ky/K* + Kl K 

are normalized, |iika| = 1; an d satisfy the orthogonality conditions, iika'K = nK| nx± = 
0. The field operators Aka satisfy the condition A^ x = A_ka, so that A^(r) = A(r). 
Substituting equations IjE.Sjl . and (|E.9J) into the Hamiltonian ()6.3|) . we obtain 

V2NLe 2 fv F \ 2 /K' K\ jb .* A „ 

K,K',k j,A,A' v 7 

Here 

P ^ A '' A ( T? J = K,K',A' • Ki,K,A) (E.ll) 



is the polarization matrix, k = q + m, and 



*%,k,a = n K A + (\/2 -II (n K A • e^J e,. 

In the case of normal incidence, Eqs. (jE.lljl . (|E.12|) result in form (|fj.7jl of polarization 
matrix. 

4. Scattering Hamiltonian, Eq. (16.61) . In the next step, we express the Fourier 
transforms of the Bose fields 6 via creation (a') and annihilation (a) operators of the 
array plasmons 

The electromagnetic field amplitudes Ak,a should also be expressed via photon creation 
(c*) and annihilation (c) operators 



A K> x(t) = J^^ + ct 



2V# V~~' A ' ~" K '\ 

Substituting these expansions into Eq. (jE.lOJ) we obtain the final form of the effective 
interaction. In the case of normal incidence, this interaction is written as 



Hlnt= ^(£|) 2 EE; ^ (k04vCK , A , (E12) 



K,K' j,X,X' 
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where hy^K.') is the Hamiltonian ()6.6|) . where A/, Aj, K are replaced by A', A, K'. 
5. Scattering cross section, Eqs. (J6.4|) - (|6.5|) . Standard procedure based on the 
Fermi golden rule leads to the following expression of the differential scattering cross 
section per unit QCB square 

1 da 1 (Vkf\ | , | rrri / £; - e f 



L 2 dtudo 7r 



(M (-3) 



(here bar denotes averaging with respect to polarization of both incident light and scat- 
tered quanta). Choose an initial ket-state |i) such that it contains an incident photon 
with momentum Kj, frequency f2j = cKi, and polarization Aj, and does not contain any 
QCB plasmon. This state can be written as i) = |Kj)^ (^) |0) p , where |0) p is the plasmon 
vacuum, |Kj)j = c K . A .|0);, and \0)i is the photon vacuum. A final bra-state (/| contains a 
scattered photon with momentum Kj, frequency Qf = cKf, and polarization A/. It con- 
tains also a QCB plasmon P with the frequency up (its quantum numbers will be specified 
below). The final state is written as (/| = (P| (££)(K/|, where (K/| = (^ Kf \ f \tyi^j ■ 
A matrix element of the interaction which enters Eq. (jE.13|) . is 

(f\H int \i) = (P\H mt \0) p , (E.14) 

where 



T7 /t^ \ itr \tr \ \ i^JgNL (ev F \ 2 sr-^ kj . 
Hint = (iV/, A f \H int \ K i5 Aj) = \7T-F7) / , /jT~ r ^" k-H-u x 



V \2cKJ ^ ,/\k~\ K f' K i- k 



XP ^(^'f)( Q ^ +a ^' (K15) 

where k|j = Kjj — (Kj / • e 3 )e 3 . In the case of normal incidence Eqs. (|E.13|) - (|E.15|) 
are equivalent to Eqs. ()6.3|) - ([6.6)1 . 
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